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EXECUTIVE  SUMMARY 


1.9  Introduction 


In  the  early  I960' a,  the  United  States  Army  supported  an  extensive  exper¬ 
imental  and  theoretical  research  program  [47]  on  radiowave  propagation  in  the 
environment  of  a  tropical,  thickly  vegetated  jungle.  The  experimental  data  ac¬ 
quired  during  this  program  were  later  used  by  Tamir  [93,107,108]  to  validate  a 
theoretical  propagation  model  which  shows  that,  for  frequencies  less  than  200 
NHs,  the  principal  mechanism  responsible  for  long-distance  propagation  is  the 
so-called  lateral  wave.  According  to  Tamir' s  theory,  a  lateral  wave  propagates 
upward  from  the  transmitting  antenna  through  the  vegetative  canopy  to  the  tree 
tops,  along  the  air-canopy  interface,  and  downward  from  the  tree  tops  through 
the  canopy  to  the  receiving  antenna. 

Recently,  the  Army  has  become  interested  in  the  development  of  several 
spread-spectrum  radio  systems  (e.g.  Packet  Radio,,  PLRS,  JTIDS)  operating  in 
the  frequency  band  200-2000  MHx .  From  the  perspective  of  radio  wave  prop¬ 
agation,  these  newer  systems  differ  from  earlier  ones  in  three  important  re¬ 
spects:  (1)  their  higher  operational  frequencies,  (2)  their  broader  spectral 
occupancy,  and/or  (3)  their  pulsed  (digital)  mode  of  operation.  The  Tamir 
model,  however,  in  valid  primarily  below  200  MHx  and,  further,  considers  only 
the  propagation  of  an  unmodulated,  time-harmonic  signal  in  a  continuous  iso¬ 
tropic  medium. 

This  report  describes  a  stochastic  radiowave  propagation  model  useful  for 
assessing  the  effects  of  forests  end  other  vegetation  upon  digital  spread-spec¬ 
trum  radio  communication  systems  operating  in  the  200  -  2000  MHx  band.  Accord¬ 
ing  to  this  model,  the  forest  is  represented  as  a  time -invariant  ensemble  of  lossy, 
randomly  positioned  and  oriented  discrete  canonical  scatterers.  Tree  trunks  are 
modelled  as  infinite ly-long,  circular,  dielectric  cylinders;  branches  as  finitely- 
long,  circular  dielectric  cylinders;  and  leaves  as  flat,  circular,  dielectric 
discs.  The  orientation  of  these  elements  must  be  specified  (statistically)  be¬ 
cause  their  scattering  is  directional.  Thus  the  model  is  anisotropic,  with  dif¬ 
ferent  properties  in  different  directions.  The  model  is  developed  in  sufficient 
detail  to  be  useful  in  the  prediction  of  attenuation  and  pulse  distortion  of  mean 
scattered  radiowaves.  Areas  for  future  study  are  identified  below  in  Sections  3, 

4,  5,  7,  and  8. 


2.0  Radiowave  Scattering  in  Discrete  Random  Media 

The  forest  is  viewed  as  a  random  ensemble  of  tree  trunks,  branches,  and 
leaves  all  having  prescribed  location  and  orientation  statistics.  Because  of 
the  inherent  randomness  associated  with  this  medium  of  discrete  scatterers,  the 
behavior  of  propagating  radiowaves  within  the  forest  cannot  be  described  by  tradi¬ 
tional,  deterministic  electromagnetic  models  but  rather  only  by  modem  stochastic 
models.  Such  models  provide  tht  basis  for  determining  the  two  most  important  char¬ 
acteristics  of  the  propagating  electromagnetic  wave:  the  mean  field  component, 
and  the  space -frequency  correlation  function  which  characterizes  the  random  (or 
fluctuating)  field  component.  See  Section  2.2.4. 

The  stochastic  electromagnetic  model  employed  by  Cyber Com  as  the  basis  for 
this  study  of  radiowave  propagation  through  the  forest  is  described  in  Section  2 
of  thin  jeport.  The  model  was  first  developed  by  Foldy  [53]  and  later  extended 
by  Lax  [55a] ,  Twersky  [64] ,  Lang  [57] ,  and  others.  The  forest  is  considered  to 
be  a  discrete  medium  representable  as  a  time -invariant  ensemble  of  randomly 


positioned  and  oriented  discrete  canonical  scatterers.  Tree  trunks  are  modeled 
a-  in finitely-long,  circular,  dielectric  cylinders!  branches  aa  finitely-long, 
circular,  dielectric  cylinders!  and  leaves  as  flat,  circular  dielectric  discs. 

The  electromagnetic  waves  scattered  by  the  canonical  scatterers  (the  cyl¬ 
inders  and  discs)  are  related  to  the  fields  induced  within  them  by  the  incident 
fields.  The  fields  induced  within  a  single  scatterer  are  related  to  the  inci¬ 
dent  electromagnetic  wave  through  the  transition  operator.  This  latter  rela¬ 
tion  arises  as  a  direct  consequence  of  the  linearity  implicit  in  maxwell* s 
equations.  Scatterers  are  also  characterized  by  the  amplitude  of  the  scattered 
field  when  the  scatterer  is  illuminated  by  a  plane-wave.  Both  the  transition 
operator  and  the  scattering  amplitude  are  employed  in  the  development  of  the 
theory.  Cyber Coo  has  utilized  the  Foldy-Lax  model  in  considering  radiowave 
scattering  within  an  unbounded  forest  comprised  of  either  two-dimensional  scat¬ 
terers  (infinitely-long,  circular,  cylinders  representing  tree  trunks)  or 
three-dimensional  scatterers,  (finitely-long,  circular,  cylinders  representing 
branches,  and  flat,  circular,  discs  representing  leaves).  For  these  media. 

Cyber Com  used  the  Twer sky  model  to  derive  mathematical  expressions  describing 
the  behavior  of  the  mean  scattered  field  and  the  space-frequency  correlation 
function.  These  expressions,  which  represent  the  principal  results  of  Section 
2,  provide  the  basis  for  subsequent  studies  addressed  in  later  sections  of  this 
report  relating  to  the  effective  permittivity  (susceptibility)  of  the  forest,  the 
specific  attenuation  of  radiowaves  propagating  through  the  forest,  the  estimation 
of  contributions  by  the  lateral  wave  propagating  above  an  anisotropic  forest,  the 
dispersion  of  broadband  radiowave  pulses,  and  the  assessment  of  the  relative 
strength  of  the  mean  and  random  components  of  the  scattered  radiowave. 


3.0  Dyadic  Scattering  Amplitude 

The  stochastic  electromagnetic  model  employed  by  cyber Com  as  the  basis  for 
this  study  of  radiowave  propagation  through  forests  is  predicated  upon  the  asser¬ 
tion  that  the  far-field  scattering  behavior  of  the  individual  canonical  scatterers 
(the  cylinders  and  discs)  can  be  characterized  by  a  (dyadic)  transition  operator 
and/or  scattering  amplitude,  m  Chapter  3  of  this  report,  mathematical  expressions 
are  derived  relating  the  scattering  amplitudes  of  the  three  canonical  scatterers 
to  the  dimensions,  orientation,  and  permittivity  of  the  scatterer  and  to  the  fre¬ 
quency  and  polarization  of  the  propagating  radiowave.  Approximate  expressions 
valid  only  in  the  low-frequency  (Bayleigh)  regime  have  been  obtained  for  all  three 
canonical  scatterers.  In  addition,  exact  expressions  valid  for  all  frequency  bands 
have  been  obtained  for  infinite-length,  circular,  dielectric  cylinders  (tree  trunks) 
and  approximate  expressions  suitable  for  the  resonant  regime  have  been  obtained  for 
discs  (leaves).  These  expressions  represent  the  principal  results  of  Section  3. 

A  resonant  model  for  finite-length,  circular,  dielectric  cylinders  (branches)  has 
yet  to  be  developed. 


4.0  Coherent  Forest  Scattering 

A  physically-appealing  representation  for  the  mean  field  cooponent  can  be 
obtained  by  postulating  that  the  mean  (or  coherent)  field  component  satisfies 
Maxwell's  equations  "in  the  mean"  and  that  the  ensemble  of  discrete  scatterers 
can  be  replaced  by  an  equivalent  continuous  medium  described  by  an  effective 
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dyadic  permittivity  ^  Zn  general,  g  ha  a  baan  found  to  dapand  upon  tha  direction 
of  radiowave  propagation  through  tha  foraat)  such  media  ara  terms?  spatially  dis- 
parsiva.  Saa  Section  4*1.  Because  £  ia  a  dyadic,  tha  foraat  is  found,  in  ganaral, 
to  ba  anisotropic. 

A  constitutive  parameter  of  tha  equivalent  continuous  medium  closely  allied 
to  tha  affective  dyadic  permittivity  £  ia  tha  affective  dyadic  susceptibility  g, 
so  defined  that  g  •  £  +  x  wh*r*  2,  ia  the  unit  dyadic.  CyberCom  has  found  that 
X  ia  directly  proportional  to  the  fractional  volume  of  tha  forest  occupied  by  tha 
scattarers.  As  a  consequence,  it  is  a  parameter  which  is  conveniently  scaled  to 
scattexer  density  and,  therefore,  is  preferred  over  g  for  the  characterization 
of  the  aquivalant  medium.  Mathematical  expressions  for  the  affective  dyadic  sus¬ 
ceptibility  of  unbounded  forests  of  tree  trunks,  branches,  and  leaves  have  been 
obtained  in  Section  4  under  the  hypothesis  that  the  orientation  distribution  of 
the  canonical  scattarers  is  asimuthally  uniform  about  tha  vertical.  Numerical 
computations  based  upon  typical  forest  parameters  show  reasonably  good  agreement 
with  values  inferred  from  measurements. 

CyberCom  has  solved  the  wave  equation  for  the  mean  scattered  field  prop¬ 
agating  through  an  unbounded  forest  consisting  solely  of  tree  trunks,  branches, 
or  leaves.  Hie  solution  reveals  that,  to  a  first  approximation,  the  horisontally- 
and  vertically  polarized  waves  propagate  independently  and  without  any  depolarisa¬ 
tion.  Further,  these  waves  propagate  with  velocities  dej-:.  lent  upon  the  inclina¬ 
tion  of  their  wave-normal  to  the  forest  floor.  These  effr :ts  are  a  consequence 
of  the  anisotropy  and  spatially-dispersive  character  of  the  equivalent  forest 
medium. 

m  general,  the  wave  propagation  constant  has  both  real  and  imaginary  com¬ 
ponents.  The  real  part  is  expressed  in  radians  par  mater;  the  imaginary  pert, 
also  called  the  specific  attenuation,  is  ejqpresssd  in  nspers  par  meter  or,  alter¬ 
natively,  in  decibels  per  meter.  Mathematical  expressions  for  the  specific  at¬ 
tenuation  of  radiowaves  propagating  through  unbounded  forests  of  tzea  trunks, 
branches,  and  leaves  have  been  obtained  in  Section  4  using  the  dispersion  re¬ 
lation  and  the  effective  dyadic  susceptibility.  Numerical  computations  based 
upon  typical  forest  parameters  show  reasonably  good  agreement  with  measured  values. 

Experimental  verification  of  tha  electromagnatic  modal  and,  ultimately,  pre¬ 
diction  of  radio  system  performance  requires  the  identification  of  omasureeble 
quantitative  parameters  to  characterize  the  forest.  Soam  of  these  parameters 
have  already  been  identified  by  CyberCom  during  the  course  of  this  study.  Some 
of  these  parameters  (e.g.,  the  size  and  relative  permittivity  of  individual  tree 
trunks,  branches  or  leaves). are  termed  "microscopic";  others  (e.g.,  number  of 
trees  per  hectare  or  branches  or  leaves  per  cubic  meter,  foliage  orientation  dis¬ 
tribution)  are  termed  "macroscopic".  Zn  the  low-frequency  (Rayleigh)  regime, 
CyberCom  has  been  able  to  derive  fairly  simple  engineering  expressions  relating 
specific  attenuation  to  salient  parameters  of  the  forest.  One  compoeite  parameter, 
especially  important  at  low  frequencies  is  tha  fractional  volian  occupied  by  the 
vegetation.  Prior  to  this  study,  no  such  expressions  were  available  for  predicting 
the  relationship  between  specific  attenuation  (dB/e)  and  ealiant  forest  parameters. 
More  work  is  required  to  identify  real-world  forest  types  and  to  derive  the  asso¬ 
ciated  susceptibilities  and  specific  attenuations. 
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5.0  Anisotropic  Forest-Slab  Model 

The  forest  model  described  in  previous  sections  of  this  report  (an  unbounded 
ensemble  of  randomly-positioned  tree  trunks,  branches  and  leaves  having  prescribed 
location  and  orientation  statistics)  can  be  refined  by  assuming  the  trees  to  be 
spread  uniformly  in  height  above  a  smooth  forest  floor  and  bounded  above  by  air. 
Earlier  efforts  [93,104,107]  have  shown,  however,  that  the  presence  of  the  ground 
complicates  the  model  significantly.  These  complications  can  be  avoided,  however, 
by  allowing  the  ground  plana  to  recede  from  the  air-forest  interface  so  that  the 
model  adopted  by  Cyber Com  provides  a  very  good  approximation  to  the  basic  slab 
model  if  both  the  transmitting  and  receiving  antennas  are  not  located  too  close 
to  the  ground  [107].  Further  work  should  include  the  effects  of  the  ground  plane. 

The  lower  half-space  representing  the  forest  is  characterised  (at  least  so  far 
as  the  mean  scattered  fields  are  concerned)  by  effective  dyadic  permittivity  £  and 
so  may  be  considered,  in  general,  to  be  electrically  anisotropic  and  spatially  dis¬ 
persive.  The  relation  between  the  effective  dyadic  permittivity  £  and  the  biophys¬ 
ical  parameters  of  the  forest  has  been  discussed  previously.  A  transmitting  antenna 
representable  as  a  vertical  tHertsian)  dipole  having  a  time-harmonic  current  moment 
of  angular  frequency  u  is  assumed  to  be  imaersed  a  distance  d  below  the  air-forest 
interface. 

The  electromagnetic  boundary  value  problem  suggested  by  the  anisotropic  half¬ 
space  model  has  been  solved  by  CyberCom  using  the  classical  approach  first  described 
by  Soamerfeld  [ill]  and  later  extended  by  Brekhovskikh  [87],  In  effecting  this 
solution,  Cyber Com  considered  the  forest  to  be  uniaxially  anisotropic  (a  consequence 
of  the  prescribed  azimuthal  uniformity  in  the  s cot te re reorientation  probability  dis¬ 
tributions)  .  Earlier,  Sachs  and  Wyatt  [lOl],  and  Tamir  and  others  [93,107]  had  con¬ 
sidered  similar,  but  isotropic,  slab-type  models  based  upon  postulated  effective 
permittivities  which  could  not  be  related  directly  to  the  biophysical  parasmters  of 
the  forest.  A  principal  conclusion  derived  from  the  earlier  efforts  was  the  exper¬ 
imentally-confirmed  conjecture  that  radiowave  propagation  over  long  distances  be¬ 
tween  antennas  within  forests  is  dominated  by  a  so-called  lateral  wave.  According 
to  those  models,  the  lateral  wave  propagates  from  the  transmitting  antenna  up 
through  the  forest  at  the  critical  angle  to  the  air-forest  interface,  through  the 
air  along  the  air-forest  interface,  and  down  through  the  forest  at  the  critical 
angle  from  the  air-forest  interface  to  the  receiver.  Because  a  substantial  frac¬ 
tion  of  the  transmission  path  of  the  lateral  wave  can  lie  in  the  dissipationless 
air,  the  transmission  loss  associated  with  the  lateral  wave  can  be  significantly 
lower  than  that  associated  with  the  direct  wave  through  the  forest.  The  former 
was  characterized  by  Tamir  as  inversely  proportional  to  the  square  of  the  distance. 

The  anisotropic  half-space  model  considered  by  CyberOom  also  shows  that  the 
propagation  of  the  mean  field  scattered  from  randomly-positioned  tree  trunks, 
branches,  and  leaves  is  dominated  at  low  frequsrcies  and  large  distances  by  the 
lateral  wave.  Preliminary  studies  reported  in  Section  7  suggest  that  in  the  fre¬ 
quency  band  200  -  2000  MHz  the  non -coherent  random  field  is  larger  that  the  co¬ 
herent  mean  field,  causing  large  spatial  fluctuations. 
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6.0  Forest  Pulse  Response 

The  anisotropic  forest-slab  model  developed  by  Cyber Coo  in  Section  5  of  this 
report  is  a  time-harmonic  model  in  the  sense  that  the  signal  radiated  by  the  an¬ 
tenna  is  a  sinusoidal  waveform  of  angular  frequency  w.  However,  because  the 
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equivalent  continuous  medium  characterized  by  the  effective  dyadic  permittivity 
i  it  a  linear  medium,  Pourier-transf^.rr  techniques  can  be  employed  to  generalize 
the  model  so  that  it  can  accommodate  such  arbitrarily-modulated  waveforms  as  the 
pulse  transmissions  employed  for  spread-spectrum  digital  systems.  Due  to  the 
complex  frequency  dependence  exhibited  by  the  effective  dyadic  permittivity  of 
the  forest,  the  model  employs  numerical  techniques  based  upon  the  fast  Fourier 
transform.  Using  this  model,  in  Section  6  of  this  report,  the  forest  pulse  re¬ 
sponse,  defined  here  as  the  vertically -polarized  component  of  the  mean  scattered 
electric  field,  is  found  for  a  transmitted  5.8  nanosecond  rectangular  pulse 
having  a  carrier  frequency  of  600  Megahertz.  The  model  can  be  used  to  assess 
medium-induced  pulse  distortion  and  inter-symbol  interference  and  can  be  easily 
extended  to  accommodate  other  field  components  of  the  mean  wave,  arbitrary  an¬ 
tenna  types,  anu  even  Doppler  effects  induced  by  terminal  antenna  movement. 


7.0  Non-Coherent  Forest  Scattering 

As  a  radiowave  propagates  through  the  forest,  power  associated  with  the  mean 
(or  coherent)  field  is  transformed  to  the  random  (or  non-coherent)  scattered 
field,  in  Section  7,  this  phenomenon  is  examined  and  the  mechanism  related  to 
specific  biophysical  parameters  of  the  forest.  To  ease  the  mathematical  burden, 
the  forest  has  been  represented  as  an  unbounded  medium  of  infinitely-long,  par¬ 
allel  tree  trunks  with  transmitting  antenna  represented  as  a  line-source  par¬ 
allel  to  the  trees.  Attention  has  been  focused  upon  the  Rayleigh-scattered  field 
intensity  of  a  radiowave  propagating  normal  to  the  trunks.  The  results  show 
that  the  intensity  of  the  non-coherent  scattered  field  increases  relative  to  that 
of  the  coherent  (mean)  field  with  increasing  distance  from  the  source  and  with 
increasing  frequency.  Is  is  desirable  to  study  the  non-coherev.t  field  further 
and  to  attempt  to  include  the  effects  of  leaves. 


8.0  Conclusions 

In  summary,  the  homogenous,  isotropic,  refracting  slab  of  a  forest  has 
been  replaced  by  Cyber Com  with  an  inhomogenous ,  anisotropic,  scattering  ensemble 
of  trunks,  branches,  and  leaves.  In  consequence: 

1.  The  lateral-wave  contribution  has  been  found  even  above  700  Megahertz. 

2.  Preliminary  validation  of  the  model  has  shown  rough  agreement  with 
experiments. 

3.  Results  have  been  obtained  for  narrow  pulse  transmission  at  600  MHz. 

4.  A  preliminary  study  ha3  emphasized  the  importance  of  the  incoherent 
component  of  the  transmitted  field. 

As  the  Cyber Com  approach  is  ambitious  in  both  scope  and  depth,  the  following 
remain  to  be  done: 

a.  The  techniques  *  '.ready  developed  must  be  exercised  to  determine  the 
relative  importance  of  forest  components  and  the  effects  of  varying 
critical  parameters.  For  example,  if  the  contribution  of  branches 

is  major,  the  difficult  characterization  of  the  scattering  properties 
should  be  advanced  in  the  resonant  region. 

b.  The  effects  of  antenna  directivity,  the  ground,  and  terminal  move¬ 
ment  should  be  incorporated  into  the  model. 

c.  The  difficult  transport  theory  for  the  important  non-coherent  scat¬ 
tered  wave  should  be  developed  as  far  as  practical. 

d.  Forest  studies  must  be  pursued  to  quantify  important  biophysical  para¬ 
meters  in  areas  of  interest. 


1.0  Introduction 


This  report  describes  a  stochastic  radiowave  propagation 
model  useful  for  assessing  the  effects  of  forests  and  other  vege¬ 
tation  upon  radio  communication  systems  operating  in  the  ?00  -  2000 
Megahertz  frequency  band.  In  this  introductory  section,  the  back¬ 
ground  leading  to  the  requirement  for  such  a  model  is  presented, 
with  a  summary  of  the  approach  employed  for  its  development  and  an 
outline  of  this  report. 

1.1  Background 

The  United  States  Army  has  had  a  long  and  continuing  interest 
in  radiowave  propagation  through  forest,  jungle,  or  otherwise 
vegetated  environments.  In  the  early  1960's,  the  Army  supported  an 
extensive  experimental  and  theoretical  research  program  [47]  on 
radiowave  propagation  in  the  environment  of  a  tropical,  thickly 
vegetated  jungle.  The  experimental  data  acquired  during  this  pro¬ 
gram  were  later  used  by  Tamir  [93,  107,  108]  to  validate  a  theoret¬ 
ical  propagation  model  which  shows  that,  for  frequencies  less  than 
200  MHz,  the  principal  mechanism  responsible  for  long-distance 
propagation  is  the  so-called  lateral  wave.  According  to  this 
theory,  a  lateral  wave  propagates  upward  from  the  transmitting 
antenna,  through  the  vegetative  canopy  to  the  tree  tops,  along  the 
air-canopy  interface,  and  downward  from  the  tree  tops  through  the 
canopy  to  the  receiving  antenna.  Prior  to  the  Army's  research 
program,  it  was  generally  believed,  primarily  on  the  basis  of  mea¬ 
surements  made  in  England  by  Saxton  and  Lane  [ 36 ] ,  that  the  excess 
(specific)  attenuation  contributed  by  the  foliage  per  meter  of  path 
length  was  independent  of  distance  and  unacceptably  high  for  the 
tactical  deployment  of  VHF  radio  sets  on  long  paths  in  forests. 

Some  of  the  experimental  data  used  by  Tamir  and  the  Army  in  refuting 
this  erroneous  contention  are  shown  in  Figure  1-1  along  with  the 
experimental  data  of  Saxton  and  Lane. 

Recently,  the  Army  has  become  interested  in  the  development  of 
several  spread-spectrum  radio  systems  (e.g.,  Packet  Radio,  PLRS, 
JTIDS)  [163].  From  the  perspective  of  radiowave  propagation,  these 


newer  systems  differ  from  earlier  ones  in  three  important  respects; 
(1)  their  higher  operational  frequencies ,  (2)  their  broader  spec¬ 
tral  occupancy,  and  (3)  their  pulsed  (digital)  mode  of  transmission. 
The  intelligent  deployment  of  these  systems  in  vegetative  environ¬ 
ments  requires  a  radiowave  propagation  model  capable  of  describing 

(1)  the  attenuating  characteristics  of  whatever  propagation  mech¬ 
anism  is  likely  to  prove  dominant  at  these  higher  frequencies,  and 

(2)  the  time-variant,  dispersive  characteristics  of  the  vegetation 
as  it  affects  wideband,  pulsed,  transmission  modes. 

To  support  this  requirement  for  enhanced  propagation  models, 
the  Army  has  again  embarked  upon  an  ambitious  research  and  develop¬ 
ment  program,  [164]  but  now  directed  toward  the  higher  frequency 
bands  (200  -  2000  MHz)  and  digital  spread-spectrum  modulation  tech¬ 
niques  typifying  the  modern  electronic  battlefield.  Experimental 
aspects  of  this  program  have  included  wideband  propagation  measure¬ 
ments  in  the  presence  of  forests  [29] ,  the  derivation  of  empirical 
path-loss  propagation  models  [45],  and,  most  recently,  the  develop¬ 
ment  of  sophisticated,  microprocessor-controlled,  spread-spectrum, 
mobile  radio  measurement  equipment.  Complementary  theoretical 
aspects  of  the  program  are  described  in  this  report. 

Prior  to  the  publication  of  this  report,  there  had  been  no 
theoretically-based  radiowave  propagation  model  for  vegetative  en¬ 
vironments  suitable  for  frequencies  in  the  band  200  -  2000  MHz  and 
capable  of  describing  not  only  transmission  path  attenuation  but 
also  pulse  distortion.  The  Tamir  model  is  valid  primarily  below 
200  MHz  and,  further,  considers  only  the  propagation  of  an  unmodu¬ 
lated,  time-harmonic  wave;  no  attempt  is  made  to  examine  the  effect 
of  frequency-dispersive  multipath  on  pulsed  signals.  Recently, 

Brown  and  Curry  [25]  have  developed  a  UHF  vegetative  model  for  air¬ 
borne,  synthetic  aperture  radars  operating  near  grazing  incidence. 
Their  model,  however,  neglects  the  effects  of  leaves  and,  further, 
does  not  consider  the  lateral  wave  found  to  be  so  important  by 
Tamir. 
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1.2  Technical  Approach 

The  foreat  is  considered  to  be  a  discrete  random  medium  repre¬ 
sentable  as  a  time- invariant  ensemble  of  randomly  position  J  and 
oriented  discrete  canonical  scatterers.  Tree  trunks  are  modelled 
as  infinitely-long,  circular*  dielectric  cylinders;  branches  as 
finitely-long,  circular  dielectric  cylinders;  and  leaves  as  flat* 
circular  dielectric  discs.  The  electromagnetic  wave  propagating 
within  this  medium  is  representable  as  the  sum  of  two  components : 
a  mean  field  component  and  a  residual  random  (or  fluctuating) 
field  component.  The  scatterers  are  assumed  to  be  characterised  by 
either  their  transition  operator  or  their  scattering  amplitudes. 
Using  a  discrete  scattering  model  originally  proposed  by  Foldy  [53] 
and  later  extended  by  Tweraky*  (64]  Lax,  (55a]  Lang  (57]  and  others* 
equations  are  derived  which  describe  the  behavior  of  the  mean  field 
and  the  space- frequency  correlation  function  of  the  random  field. 

In  this  report*  attention  is  concentrated  on  the  behavior  of 
the  mean  field*  It  is  shown  that  with  regard  to  this  component* 
the  ensemble  of  discrete  scatterers  can  be  represented  by  an  equiv¬ 
alent  continuous  medium  characterised  by  an  effective  dyadic  per¬ 
mittivity  which  can*  in  turn*  be  related  directly  to  the  sise* 
shape,  orientation,  number  density  and  permittivity  of  the  scat¬ 
terers  themselves.  This  equivalent  nedium  is  then  used  to  determine 
the  wave  propagating  within  an  unbounded  lores t ,  and  to  define  an 
anisotropic  forest  slab  model  analogous  to  the  isotropic  slab  model 
introduced  by  Sachs  and  Wyatt  (101]  and  studied  by  Tamir  [107]. 

The  electromagnetic  fields  scattered  from  the  tree  trunks, 
branches*  and  leaves  are  assumed  to  be  linearly  related  to  the 
mutually-induced  currents  excited  within  them  by  the  propagating 
radio  wave.  As  a  consequence*  the  forest  transmission  channel  be¬ 
tween  the  transmitter  and  receiver  can  be  considered  linear*  and 
Fourier  spectral  techniques  employed  to  extend  the  time-harmonic 
forest  slab  model  to  encompass  arbitrarily-modulated  transmissions. 


•Chapter  7  shows  ths  importance  of  the  random  or  noncoherent  field. 
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1.3  Scope 

This  report  uses  a  stochastic  model  for  describing  the  be¬ 
havior  of  radiowaves  propagating  through  forests  or  other  vegetated 
regions.  The  discrete  scattering  theory  supporting  the  model  is 
developed  in  Sections  2  and  3.  In  Section  4,  the  forest  is  repre¬ 
sented  (at  least  so  far  as  the  mean  scattered  field  is  concerned) 
by  an  equivalent  continuous  medium,  and  characterized  by  an  effec¬ 
tive  dyadic  permittivity  which  is  directly  related  to  the  biophy¬ 
sical  parameters  of  the  forest.  Also  in  Section  4,  expressions  for 
the  specific  attenuation  are  developed, and  evaluated  for  typical 
biophysical  parameters  of  a  deciduous,  hardwood  forest.  In  Section 
5,  an  anisotropic  forest  slab  model  is  introduced  and  the  relative 
contributions  of  the  different  propagation  modes  (direct  wave,  re¬ 
flected  wave,  and  lateral  wave)  assessed.  The  anisotropic  forest 
slab  model  is  applied  in  Section  6  to  accommodate  arbitrarily- 
modulated  waveforms;  numerical  results  are  presented  for  the  case 
of  a  broadband  r-f  pulse.  Section  7  considers  the  relative  impor¬ 
tance  of  the  random  (fluctuating)  field. 

Although  the  model  described  herein  can  accommodate  transmit¬ 
ting  and  receiving  antennas  of  arbitrary  directivity  and  polariza¬ 
tion,  this  initial  study  emphasizes  vertically-polarized,  electri- 
cally-short,  linear  antennas.  Further,  although  Doppler  effects 
due  either  to  wind-induced  tree  motion  and/or  terminal  (antenna) 
motion  can  also  be  accommodated  by  the  model,  they  are  disregarded 
here.  Finally,  although  the  validity  of  the  model  is  supported 
tentatively  by  comparison  with  measurements,  this  will  be  pursued 
further  as  the  model  is  developed  and  additional  experimental  data 
become  available. 
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2.0  Radiowave  Propagation  in  Discrete  Random  Media 

The  basic  methodology  for  describing  radiowave  propagation 
through  the  forest  is  formulated  in  this  section.  Because  of  the 
complexity  of  vegetation,  stochastic  or  random  methods  are  used 
rather  than  deterministic  techniques.  The  forest  is  viewed  as  a 
random  ensemble  of  trunks,  branches  and  leaves.  Equations  are 
then  developed  for  the  mean  and  correlation  of  the  electromagnetic 
fields.  From  these  averages,  physical  quantities  such  as  attenua¬ 
tion,  propagation  delay  and  pulse  dispersion  are  obtained. 

2.1  Representation  of  Scattered  Electromagnetic  Fields 

The  stochastic  approach  takes  the  viewpoint  that  the  field 
quantities  are  composed  of  a  mean  component  and  a  random  (or  fluc¬ 
tuating)  component.  If  ip  represents  an  electromagnetic  field 
quantity  such  as  E  or  H,  it  is  in  general  a  function  of  the  loca¬ 
tion  and  orientation  of  the  vegetative  components.  Since  these 
are  randomly  located  and  oriented,  the  field  quantity  (i  is  a  random 
variable.  As  such,  it  can  be  broken  up  into  a  mean  component,  <^># 
and  a  fluctuating  component,  $,  i.e. 

«  <^>  -r  $  (2-1-1) 

Here  the  brackets,  <  >,  have  been  used  to  denote  an  ensemble 
average.  Taking  the  average  of  Equation  (2-1-1)  shows  that  <$>»0, 
a  reasonable  result  since  the  mean  has  already  been  extracted 
from  \p. 

The  correlation  function  of  ip  (see  Section  2.2.4)  can  be  de¬ 
termined  by  multiplying  ip  and  its  complex  conjugate  i|>*  and  averag¬ 
ing  over  the  ensemble.  By  using  Equation  (2-1-1)  and  the  fact  that 
the  average  fluctuations  are  zero,  the  correlation  function  of  <|> 
can  be  written  as 

m  <\p><\p*>  +  <^*> 

Here  the  correlation  is  broken  up  into  two  components:  a  mean  part, 
<i|>>  <^*>,  and  a  fluctuating  component  <W*>.  Wh®n  absorption  is 
important  than  scattering,  the  fluctuating  component  is  small. 

This  is  typically  the  case,  when  the  wavelength  is  large  compared 


a 
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to  the  size  of  the  scat terers *  In  this  regime,  the 

correlation  can  be  approximated  by  J; 

*  *  '•T; 

•>  <\^><^*>  (2-1-3) 

It  is  this  approximation  that  is  inherent  in  early  attempts  at 

forest  modeling  [92,  93,  100,  101,  106,  107,  108].  The  forest  was  j 

replaced  by  a  dielectric  slab  and,  in  effect,  mean  fields  were 

computed. 


The  development  of  equations  for  the  mean  and  correlation  can 
be  approached  in  two  different  ways.  They  are  the  continuous  and 
discrete  modeling  procedures.  The  continuous  approach  uses  the 
mean  and  correlation  of  the  effective  forest  permittivity  as  input 
information  to  the  model.  The  discrete  approach,  on  the  other  hand, 
uses  the  scattering  amplitudes  of  individual  scatterers,  as  well  as 
the  position  and  orientation  statistics  of  scatterers.  Examples  of 
modeling  by  the  continuous  approach  are  provided  by  Keller  156] , 
Besieris  and  Kohler  [49]  and  Tsang  and  Kong  [62],  while  the  discrete 
approach  has  been  used  by  Twersky  [63,  64],  Ishimaru  [54,  55]  and 
Lang  [57,  58].  The  continuous  modeling  approach  tends  to  be  some¬ 
what  simpler  than  the  discrete  method,  nils  is  due  to  the  need  for 
scattering  amplitudes  of  individual  scatterers  in  the  discrete 
method.  The  discrete  technique,  however,  provides  a  closer  connec¬ 
tion  with  reality.  Individual  scatterer  sizes  and  dielectric  con¬ 
stants  can  be  measured;  the  effective  dielectric  permittivity  and 
its  correlation,  as  required  by  the  continuous  approach,  are  diffi¬ 
cult  quantities  to  determine. 

For  the  discrete  case,  an  approximate  equation  for  the  mean 
field  was  obtained  by  Foldy  [53]  for  dipole  scatterers  and  later 
by  Lax  [55a]  for  resonant-size  scatterers.  The  mean  equation  is 
valid  when  the  fractional  volume  of  vegetation  is  small.  This  cri¬ 
terion  appears  to  be  satisfied  within  most  forests  (see  Section  4.3). 
Twersky  [63]  subsequently  derived  an  approximate  equation  for  the 
correlation  which  again  is  valid  when  the  fractional  volume  is 
small.  It  is  these  basic  equations  that  will  be  used  to  obtain  the 
macroscopic  effects  of  vegetation  on  radiowave  propagation. 
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2.2  Two-Dimensional  Scattering  Media 


Scattering  by  an  aggregate  of  parallel  dielectric  cylinders 
is  considered  in  this  section.  Because  of  the  planar  symmetry  of 
the  problem,  scalar  rather  than  vector  equations  can  be  used.  This 
substantially  simplifies  the  analysis  and  thus  makes  the  derivation 
of  the  mean  and  correlation  equations  more  transparent.  In  addi¬ 
tion,  the  problem  provides  a  good  model  for  a  forest  consisting 
wholly  of  trunks  or  or.e  in  which  trunks  have  the  dominant  effect 
on  the  channel  properties. 

2.2.1  Model  Formulation 

Consider  a  collection  of  N  identical  parallel  dielectric 
cylinders  each  having  complex  relative  permittivity  ^  and  cross- 
sectional  surface  area  S^.  These  scatterers  are  shown  in  Figure 
2-1.  In  this  figure  a  cross-sectional  view  of  the  forested  region 
is  shown.  The  forest  is  totally  contained  within  the  area  S. 

Analysis  of  the  problem  begins  by  considering  Maxwell's  time- 
harmonic  equations  (with  an  expijwt)  time  dependerce  assumed) 


VxE  «  -j<i)y0H 
7xH  =  jwe0e(xt)E  + 


(2-2-1) 


where  y0  and  e©  are  the  permeability  and  permittivity,  respectively, 
of  free  space,  t'  i.s  the  angular  frequency,  and 


t  (x  '  a  1  r  l  xs(x.) 

i-1  3 


(2-2-2) 


is  the  relative  permittivity  of  the  scatterers  as  a  function  of 
transverse  vector  position  x^.  The  susceptibility  function  Xj (xt) 
represents  the  susceptibility  of  the  scatterer  as  a  function 
of  position.  It  is  given  by 


Kj(5t> 


(  Xi  t  inside  scatterer 


I. 


outside  scatterer 


(2-2-3) 
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Here  is  the  susceptibility  of  a  scatterer,  which  is  related  to 
ep  x*  “  H  “  l- 

To  reduce  the  vector  field  problem  to  a  scalar  field  problem, 
the  current  Jg  is  assumed  to  have  the  following  form 

Jg  -  J2(xt)*»  *  (2-2-4) 

It  has  been  tacitly  assumed  in  writing  Equation  (2-2-4)  that  the 
gcatterers  are  parallel  to  the  z-axis.  The  electromagnetic  fields 
are  decomposed  into  tr&acverse  and  z-directed  components  so  that 


E  *  E.  +  E„z • 

—  -v  I- 

H  -  H.  +  H  Z# 

—  -t  z — 


(2-2-5) 


Substitution  of  Equations  (2-2-4)  and  (2-2-5)  into  Equation  (2-2-1) 
reveals  that  E_  and  H*.  are  the  only  field  components  excited,  and 
that  they  satisfy 

IV|  +  k*e(xfc)  JEx(xt)  -  (2-2-6) 


Here  is  the  transverse  del  operator  and  k 2  ■  w2y0e0.  The  problem 

has  been  scalar ized.  As  a  result,  the  complete  field  behavior  can 

be  obtained  from  the  scalar  field  E_. 

z 

It  is  convenient  at  this  point  to  introduce  an  operator  nota¬ 
tion.  This  notation  will  highlight  the  important  aspects  of  the 
development  while  suppressing  unimportant  details.  It  will  also 
make  the  parallelism  that  exists  between  the  scalar  and  vector 
models  more  apparent.  The  operator  notation  is  introduced  by  de¬ 
fining  the  following  quantities: 

L  -  -(V2  +  k2)  Vj  *  k*Xj(St>  (2-2-8) 


f  -  Ex  g  -  -ju>u«Jr  (2-2-9) 

By  using  this  notation  in  Equation  (2-2-6)  the  wave  equation  can  be 
written  as 


(2-2-10) 


*«'  is  a  unit  vector  in  ths  s  direction. 
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For  later  use,  the  field  Y  can  be  decomposed  into  incident  and 
scattered  parts.  This  is  written  as 


*1  ♦ », 


(2-2-11) 


where  Y i  is  the  incident  wave  and  Y^  is  the  scattered  wave.  The 
incident  field  Y^  is  the  field  that  would  exist  if  the  scatterers 
were  not  present,  thus 

LYa  -  g  (2-2-12) 

2.2.2  Single-Scatterer  Characterisation 

Before  the  N-particle  problem  can  be  addressed  the  single 
scatterer  must  be  characterised  electromagnet ically .  For  this  pur¬ 
pose,  the  transition  operator  and  scattering  amplitude  are  intro¬ 
duced  in  this  section. 

Consider  a  single  scatterer  located  at  the  origin  as  shown  in 
Figure  2-2.  The  electric  field  satisfies  the  wave  equation  given 
by  Equation  (2-2-10)  with  N*l.  This  equation  is 


where 


with 


(L  -  V)* 


V  -  k*X(xt) 


-jwy0J, 


(2-2-13) 


(2-2-14) 


(2-2-15) 


X(xt) 


inside  S. 


outside  S. 


(2-2-16) 


Lower-case  notation  has  been  employed  for  the  field  as  a  reminder 
that  only  one  scatterer  is  being  considered. 

The  transition  operator  will  now  be  introduced.  This  operator 
relates  the  induced  sources  within  the  scatterer  to  the  incident 
field  upon  the  scatterer.  To  develop  this  relation,  the  field  is 
separated  into  incident  and  scattered  components  so  that 
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(2-2-17) 


♦  =  ♦j,  +  ♦, 


»  *  g 


Substitution  of  Equation  (2-2-17)  into  Equation  (2-2-13)  gives  the 
equation  for  the  scattered  field,  viz. 


geq  -  v* 


(2-2-18) 


The  equivalent  sources  that  generate  the  scattered  field,  g  ,  are 
given  by  Vip.  The  transit  ion  operator  T  can  now  be  defined  as 

geq  -  T*a  (2-2-19) 

From  the  above  equation,  it  is  seen  that  once  the  incident  field  is 
specified,  the  induced  sources  can  be  calculated  if  T  is  known. 

Since  bounded  yields  bounded  induced  source  distributions,  this 
implies  that  T  is  a  bounded  operator.  As  a  result,  T  can  be  repre¬ 
sented  in  integral  form  as 


?eq(^t)  =  jd*tt  (2St'5t)^i(2St) 


(2-2-20) 


The  kernel  t ( x t , x^)  is  analogous  to  the  time-varying  impulse  re¬ 
sponse  h(t,t).  One  can  show  that  t(xt,x£)  *  0  for  xfc  or  ^  outside 
of  Sp.  This  just  means  that  the  equivalent  sources  are  located 
within  the  boundaries  cf  the  scatterer. 

In  the  low-frequency  case,  where  the  wavelength  is  large  com¬ 
pared  to  the  characteristic  size  of  the  scatterer’ s  cross-section, 
the  equivalent  source  distribution  looks  like  a  line  current.  In 
this  case,  the  transition  operator's  kernel  can  be  written  as 


t(xt,x£)  *  ko<x6  (x*.)  6  (x^) 


— t  ' 


(2-2-21) 


The  introduction  of  the  coefficient  a  is  motivated  by  three- 
dimensional  scatterers  for  which  a  is  the  polarizability  [57]. 

The  scatterer  can  be  characterized  in  an  alternative  manner 
by  specifying  its  far-field  response  to  a  unit-amplitude  incident 
plane  wave.  If  the  incident  plane  wave  is  given  by 


*i<xt)  =  exp{-jk0i*xt) 


(2-2-22) 


then  the  scattered  field  in  the  radiation  zone  has  the  form 
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♦8<*t>  “ 


9 


(2-2-23) 


exp{-jk0xt) 


where  i  and  o  (=  xt/xt)  are  unit  vectors  in  the  direction  of  the 
incident  wave  and  the  observation  point  respectively  as  is  shown 
in  Figure  2-2.  Equation  (2-2-23)  serves  as  a  defining  equation 
for  the  scattering  amplitude , 

The  scattering  amplitude  and  the  transition  operator  are 
related  to  each  other.  The  relationship  is 


f (o,i)  =  yt(k0o,k0i)  ,  y  *  e“^  (2-2-24) 

—  —  «...  —  A  e 

where  t(kfc,k£)  is  the  Fourier  transform  of  t(xt,x£)  with  respect 
to  xfc  and  x£ ,  i . e . , 

f  “  ki-x') 

t(kt,k^)  =  ( 2 7T )  2  I  dxtdx£  t  (xt,x*)e  r  r  r  (2-2-25) 

The  relationship  between  f  and  t  is  derived  in  Appendix  B. 

In  developing  the  multiple-scattering  equations  the  transi¬ 
tion  operator  for  a  scatterer  not  located  at  the  origin  will  be 
needed.  Denote  the  transition  operator  kernel  for  the  jtl1  scat¬ 
terer  by  tj(xt,x£).  It  can  be  related  to  the  transition  operator 
of  the  scatterer  located  at  the  origin.  If  the  jth  scatterer’ s 
center  is  located  at  Xt^  where  Xfcj  is  measured  from  the  origin,  then 
by  a  simple  shift  of  the  incident  field  and  the  induced  charge,  one 
finds 

t  j  ( x ,  x^ )  *  t(xt-Xtj,x£-Xtj)  (2-2-26) 

Here  t(xt,x£)  is  the  transition  kernel  for  a  scatterer  located  at 
the  origin.  It  should  be  noted  here  t  is  a  function  of  0..  -  the 
orientation  angle  of  the  j*-*1  scatterer. 


2.2.3  Mean-Wave  Equation 

An  approximate  equation  for  the  mean  field  can  be  derived 
using  the  Foldy-Lax  method  [57] .  This  method  is  valid  in  the 
case  of  small  fractional  volume. 
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Consider  the  total  field  ¥  at  some  point  xt  to  consist  of  the 
incident  field,  and  the  various  scattered  fields  from  the  in¬ 
dividual  particles,  i.a. 

r  f  i ) 

V  «  +  I  M3'  ,  LTi  «  g  (2-2-27) 

j-i 

where  io  the  scattered  field  from  the  particle,  The  tvan- 

si  Mon  operator  for  the  particle  can  then  he  used  to  relate  the 
jth  scattered  field  to  the  incident  field  on  the  j**1  particle. 

This  is  expressed  mathematically  as 

-  Tjt(j)  (2-2-28) 

where  Tj  is  the  transition  operator  for  the  particle  and 
is  the  field  at  the  location  of  jth  particle  with  the 

particle  removed.  Substituting  Equation  (2-2-28)  into  Equation 
(2-2-27)  and  multiplying  from  the  left  by  L,  gives  the  following 
equation 

N  /-u 

LY  *  g  +  l  T.'?'3'  (2-2-29) 

j*l  3 

Averaging  the  above  equation  over  the  ensemble  of  ecatterer 
configurations  yields 


L<f>  -  g  +  l  <T.y*3)>  (2-2-30) 

j-1  3 

where  the  fact  that  L  is  deterministic  has  been  used.  To  obtain 
an  equation  for  the  mean  field  <y>,  the  Foldy-Lax  approximation 


y(:i)  «  <f>j  (2-2-31) 

is  introduced.  This  s^ys  that  the  field  incident  upon  the  j**1 
scattercr  is  approximately  equal  to  the  mean  field  at  the  jth 
acatterer.  It  can  be  shown  by  scaling  techniques  that  this  ap¬ 
proximation  is  valid  when  the  fractional  volume  occupied  by  the 
acatterer s  is  small.  Using  this  approximation 


<Tjj,(j)>  v 

3 


<Tj<y>> 


<Tj><y>j 


(2-2-32) 
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Putting  this  approximation  into  Equation  (2-2*30)  yields  the 
following  equation  for  the  mean  field 

L<f>  -  I<Tj><f3  "  9  (2-2-33) 

The  equation  for  the  mean  field  can  be  simplified  by  ex¬ 
plicitly  writing  cut  <Tj>  and  assuming  that  all  particles  are 
identically  distributed  in  location  and  orientation.  Doing  this, 
the  expression  for  <T^>  becomes 

<Tj>  -  |  dst  J^depts^ejT^s^e)  (2-2-34) 

s  0 

where  p(st,6)dstd0  is  the  joint  probability  that  the  position 
vector  Xt  will  take  on  a  value  close  to  st  and  simultaneously  the 
orientation  angle  0  will  take  on  a  value  close  to  6.  Note  that 
the  average  is  identical  for  all  particles  since  the  probability 
density  is  independent  of  j. 

Before  proceeding.  Equation  (2-2-34)  can  be  written  in  a 
somewhat  more  convenient  form  by  suppressing  the  orientation 
average  and  writing 

<T.,>  -  |  dstp(st)Tj(st)  (2-2-35) 

s 

where  Tia^) ,  the  orientation  average  over  0  at  a  particular  st, 
is  given  by 

r  2ir 

f(st)  -  I  d0p(0|8t)T(st,e)  (2-2-36) 

0 

A  simplified  equation  for  the  mean  can  now  be  written.  Using 
Equation  (2-2-35)  in  Equation  (2-2-33)  and  noting  that  the  terms 
of  the  start  are  identical  gives 

L<f>  -  |  dstp(st)f (st)<4'>  «  g 
e 

where  the  number  density  of  the  scatterers  p(s^)  is  just 


(2-2-37) 


The  notation  in  this  equation  can  be  simplified  by  defining  the 
mean  wave  operator  i 

X  «  L  -  j  (2-2-39) 

s 

Equation  (2-2-37)  then  becomes 

f <¥>  -  g  (2-2-40) 

The  mean  equation  as  given  above  is  in  its  most  compact 
form.  For  later  use,  it  will  prove  expedient  to  write  it  in  the 
more  explicit  form 

(Vt2+k| )  <V  (xt)  >+|dstp(s!|)  | dx£t  (xt-st ,x^-st)  <f  (x£)  >— g  (xt )  (2-2-41) 

s 

This  equation  is  an  integro-differential  equation;  however,  it 
is  deterministic  and  not  random.  Further,  the  equation  substan¬ 
tially  simplifies  in  the  case  of  dipole  scatterers  as  may  be  seen 
by  substituting  Equation  (2-2-21)  into  Equation  (2-2-41)  to  obtain 

[^2+k*(l  +  p(xt)a]<f(xt)>  -  -g(xt)  (2-2-42) 

where  it  is  understood  that  p(xt)  *0  for  outside  of  S.  This 
simplification  is  possible,  however,  only  at  low  frequencies 
(large  wave  lengths)  where  the  incident  field  is  essentially  con¬ 
stant  over  the  scatterer.  Since  this  low-frequency  approximation 
involves  the  assumption  that  the  incident  wavelength  must  be  large 
compared  to  the  sice  of  the  scatterers,  scattering  from  resonant- 
size  scatterers  must  be  treated  by  the  more  general  mean  Equation 
(2-2-41) . 


2.2.4  Space-Frequency  Correlation  Function 

A  basic  quantity  of  interest  for  characterizing  such  com¬ 
munications  channels  as  the  forest  medium  is  the  space- frequency 
correlation  function  of  the  scatterer  field,  <^(xt,w)4>*(xt,w) . 

This  is  the  ensemble  average  of  a  field  component  at  xfc  and  angu¬ 
lar  frequency  u>  times  the  conjugate  of  the  field  component  evaluated 
at  x^.  and  w.  At  a  fixed  point  (xt»xt) ,  the  resulting  frequency 
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correlation  function  indicates  the  effective  bandwidth,  and  its 
Fourier  transform  gives  the  delay  spread.  At  a  fixed  frequency 
(<»>■£>) ,  the  width  of  the  space  correlation  function  indicates  the 
minimum  separation  required  for  efficient  space  diversity. 

In  Section  2.1,  it  was  pointed  out  that  if  the  medium  was 
mostly  absorptive  in  nature  rather  than  scattering,  the  space- 
frequency  correlation  function  can  be  obtained  as  in  Equation 
(2-1-3)  by  way  of  the  mean  field  as  follows 

<f (x1:,w)**(xt,w)>  *  <»(xt,w)><f*(xt,w)>  (2-2-43) 

In  general,  hovrever,  as  the  frequency  increases,  scattering  be¬ 
comes  more  important  and  Equation  (2-2-43)  no  longer  holds.  The 
correlation  of  the  field  must  then  be  obtained  from  first  princi¬ 
ples  as  in  the  case  of  the  mean. 

An  equation  for  the  correlation  can  be  obtained  by  parallel¬ 
ing  the  development  of  the  mean  wave  equation.  The  equation  must 
again  be  derived  under  the  assumption  of  small  fractional  volume. 
However,  there  is  no  restriction  on  the  absorptive  or  scattering 
properties  of  the  particles.  Because  of  the  amount  of  technical 
detail  involved,  the  derivation  of  this  equation  has  been  rele¬ 
gated  to  Appendix  C. 

The  correlation  equation  obtained  in  Appendix  C  is 

£ -  |  dstp(st)T(st)T*(st)<fy*>  «  gg*  (2-2-44) 

s 

where 

r2it 

T(st)T*(stf  «  I  d0p(elst)Tj(8t,0)T*(st,e)  (2-2-45) 

0 

Note  that,  as  in  Equation  (2-2-43),  the  unstarred  and  starred 
quantities  are  functions  of  unhatted  and  hatted  quantities,  re¬ 
spectively.  The  mean  operator  £  is  given  in  Equation  (2-2-39) . 

The  correlation  equation  can  be  put  in  a  slightly  different  form 
by  multiplying  from  the  left  by  i”1****  where  the  inverse  mean 
operator  X”1  is  given  by 
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f’Vx*)  -  Jdx£G(xt,x£)iM2£> 


(2-2-46) 


and  the  kernel  of  the  integral  operator  is  the  mean  Green's  func¬ 
tion  which  satisfies 


£G(xt,x£)  - 


(2-2-47) 


The  result  is 


<ff*>  .  <*><¥*>  -  £"1£*“1|dstp(st)T(et)T*(st)<f1'*>  (2-2-4 


This  is  equivalent  to  the  correlation  equation  given  by  Ishimaru 
[55] ,  although  his  equation  is  given  for  the  three-dimensional 
scalar  case. 

As  in  the  case  of  the  mean.  Equation  (2-2-46)  will  now  be 
written  out  more  explicitly 

<* (xt)f*(xfc)>  -  <f(xt)><!*(xt)> 

¥ 

-  Jdx^dx^dx^dxJdstp(st)G(xt,x^)G*(xt,x^)  (2-2  49) 

* t  ^V-t » it’it }  (it“5t  *  *t~±t ]  * <f  ]  f  (*t ]  > 

In  the  above  equation  it  has  been  assumed  that  a >  ■  u>  and  then  this 
variable  has  been  suppressed. 

As  in  the  case  of  the  mean -wave  equation,  when  the  frequency 
is  low,  the  correlation  Equation  (2-2-49)  can  be  substantially 
simplified.  By  using  Equation  (2-2-21)  in  Equation  (2-2-49),  it 
is  found  that 


<¥(xt)¥*(xt)>  «  <f(xt)><1'*(xj)> 


-k$ | a | 2  j  dst p (st ) G (xt , st ) G* (xt , st ) • <f (st ) f * (st) 


(2-2-50) 


The  solution  to  this  equation  will  be  found  in  Chapter  7  for  an 
unbounded  medium  (S  ■  ™)  with  constant  density  p  and  x t  ■  ^ 
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2,3  Three-Dimensional  Scattering  Media 


2.3.1  Model  Formulation 

In  the  preceding  sect ion ,  equations  were  derived  for  the 
mean  and  correlation  of  the  field  in  a  medium  of  two-dimensional 
scatterers.  Attention  will  now  be  focused  on  ensembles  of  three- 
dimensional  scatterers.  These  scatterers  will  be  used  to  model 
leaves  and  branches.  Because  of  the  three-dimensional  character 
of  the  scatterers,  the  problem  can  no  longer  be  scalarized,  thus 
the  vector  wave  equation  must  be  used. 


The  configuration  of  scatterers  to  be  discussed  is  shown  in 

Figure  2-3.  There  are  N  identical,  nonaligned  scatterers  contained 

within  volume  V.  Each  scatterer  has  volume  and  relative  di- 

P 

electric  constant  e^.  A  deterministic  background  medium  is  also 
assumed  to  be  present.  The  relative  permittivity  of  the  composite 
medium  *8  given  by 


£  (x) 


N 

+  I  X,  (x) 

j-1  3 


(2-3-1) 


where  (x)  is  the  relative  permittivity  of  the  background  medium 
which  is  assumed  to  be  unity  inside  V  and  arbitrary  outside  V, 
i.e. 


1  ,  x  inside  V 

eb(x)  -  '  (2-3-2) 

cb(x)  ,  x  outside  V 


By  appropriately  choosing  V  and  E^fx),  a  half-space  or  slab  con¬ 
figuration  of  scatterers  can  be  constructed.  In  addition,  the 
ground  can  be  accounted  for  through  the  background  medium.  The 
Xj  (x)  used  in  Equation  (2-3-1)  is  the  susceptibility  of  the  j**1 
scatterer.  It  can  be  written  as 


Xj  (x) 


IX£  ,  x  inside  jth  scatterer 
0  ,  x  outside  jth  scatterer 


(2-3-3) 


Here,  as  in  the  case  of  the  two-dimensional  scatterers,  x£  is  the 
bulk  susceptibility  of  the  scatterer  material. 
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A  vector  wave  equation  for  the  electric  field  can  be  obtained 
by  using  Equation  (2-3-1)  in  Maxwell's  equations  as  given  by  Equa¬ 
tion  (2-2-1).  The  resulting  equation  is 


N 

VxVxE  -  kje.  (x)£  -  l  k Jx .  (x) E  -  -jwy.J  (2-3-4) 

°  ~  ~  j-1  3  "  ” 

An  operator  notation  is  again  introduced.  By  using  the  fol¬ 
lowing  notation 


L  -  Vx(VxI)  -  k|eb(x)I 

(2-3-5) 

Vj  »  ^Xj(x)  I 

(2-3-6) 

¥  «  E  £  *  “jwy.J 

(2-3-7) 

Equation  (2-3-4)  can  be  written  as 

<£  -  (je«>  •*  -  a 

■4  J 

(2-3-8) 

3 


The  double  underlined  notation  is  used  to  signify  a  dyadic;  thus 
L  and  v_.  are  dyadic  operators.  The  unit  dyadic  I  has  been  intro¬ 
duced  in  the  definition  of  L  and  V  ^ .  As  its  name  indicates,,  it 
has  the  property  that  I«T  *¥•!«¥  [see  Appendix  DJ .  This  pro¬ 
perty  is  used  to  reduce  the  abstract  form  of  the  vector  equation 
as  given  in  Equation  (2-3-8)  to  the  standard  form  as  given  in 
Equation  (2-3-4). 

One  of  the  principal  advantages  of  the  abstract  operator  no¬ 
tation  is  now  evident.  A  comparison  of  the  vector  wave  equation 
as  given  by  Equation  (2-3-8)  with  the  two-dimensional  scalar  wave 
equation  as  given  by  Equation  (2-2-10)  shows  that  the  two  equations 
have  basically  the  same  form.  Thus  derivations  carried  out  in  the 
vector  case  will  parallel  those  already  performed  in  the  scalar 
case. 


2.3.2  Single-Scattering  Problem 

Before  a  mean  equation  can  be  derived  for  the  field,  the 
scatterers  must  be  elect romagnetically  characterized.  The  dyadic 
transition  operator  and  the  dyadic  scattering  amplitude  can  be 


used  for  this  purpose.  The  development  in  this  subsection  paral¬ 
lels  the  development  for  the  scalar  two-dimensional  problem  as 
given  in  2.2.2. 

Consider  a  three-dimensional  acatterer  located  at  the  origin. 
The  electric  field  £  satisfies  the  equation 

{&  -  V)  •  ¥  -  a  (2-3-9) 

where 

V  -  k*X(x)£  (2-3-10) 

and  L,  £  are  defined  in  Equations  (2-3-5)  and  (2-3-7)  respectively. 
Here,  as  before,  X(x)  is  X&  inside  the  scatterer  and  zero  outside 
the  scatterer.  The  orientation  of  the  scatterer  as  specified  by 
the  polar  angle  6  and  the  azimuth  angle  ♦  is  contained  implicitly 
in  the  definition  of  X(x)  ,  i.e.,  X(x)  *X(x,0,$). 

To  define  the  transition  operator,  the  field  is  broken  up  into 
a  free-space  component  and  a  scattered  component  ^  so  that 

t  *  ±i  ♦  ^  (2-3-11) 

where  *  g  and  Here  as  in  the  scalar  case,  V*^  can 

be  viewed  as  equivalent  sources  gCg  that  create  the  scattered 
field.  The  dyadic  transition  operator  can  now  be  introduced.  This 
operator  acts  on  the  incident  field  to  produce  the  equivalent 
Bources ;  thus 

Seg  »  T  *  ij  (2-3-12) 

Since  the  relationship  is  linear  and  bounded.  Equation  (2-3-12)  can 
be  written  explicitly  as 

^qU)  =  jdx't (x,x’) (x')  (2-3-13) 

where  t(x,x‘)  is  called  the  transition  operator  kernel.  The  know¬ 
ledge  of  t  completely  specifies  the  scattered  field  once  the  inci¬ 
dent  field  is  given. 

The  relationship  between  ^  and  c^g  is  substantially  simpli¬ 
fied  in  the  low-frequency  limit.  In  this  limit  the  induced  source 
distribution  can  be  approximated  by  an  electric  dipole.  Thus 


J 


W 
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(2-3-14) 


Seq  *  “  *  - jwu. I jw£fi (x) 1 

where  £  is  the  electric  dipole  moment  of  the  induced  source  dis¬ 
tribution  and  6(x)  is  the  three-dimensional  Dirac  delta  function. 

The  dipole  moment  £  can  be  directly  related  to  the  incident  elec¬ 
tric  field  ^  by  the  dyadic  polarizability  gu  This  relationship 

is 

£  *  e«a  *  £0  (2-3-15) 

Comparing  this  equation  with  Equations  (2-3-12)  and  (2-3-13)  the 
transition  operator  and  its  kernel  can  be  written  as 

T  -  kjjafi  (x)  (2-3-16) 

and 

t(x,x* )  -  k*a6(x)5(x')  (2-3-17) 

The  far-field  behavior  of  the  scatterer’ s  response  to  a  plane 
wave  can  als~)  be  used  to  characterize  the  scatterer.  Consider  a 
plane  wave 

*  aie"jko~*-  (2-3-18) 

incident  upon  the  scatterer.  Here  q?  is  a  unit  polarization  vector 
and  i  is  a  unit  vector  in  the  direction  of  propagation.  The  scat¬ 
tered  field  in  the  far-field  zone  of  the  scatterer  is  written  as 

*  jkoX 

jjj—  (x)  %  f (o,i)  •  2?  -  -  —  ,  o  *  (2-3-19) 

—b  —  x  x  —  |  x  | 

where  f(o,i)  is  the  dyadic  scattering  amplitude.  Here  o  is  a  unit 
vector  in  the  direction  of  observation  and  x*  |x| . 

The  scattering  amplitude  f  can  be  understood  more  easily  in 
terms  of  its  components.  To  do  this,  incident  and  scattered  polar¬ 
ization  vectors  need  to  be  defined.  An  arbitrary  incident  polar¬ 
ization  vector  2?  can  be  decomposed  into  two  orthogonal  components; 
call  them  h?  and  vV  indicating  horizontal  and  vertical  polariza¬ 
tion,  respectively.  Similarly,  the  scattered  field  can  be  decom¬ 
posed  into  two  mutual  orthogonal  components  hf  and  v*.  The  no- 
tation  has  been  employed  that  a  superscript  zero  indicates  a  unit 
vector. 
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By  using  the  unit  polarization  vectors  introduced  above  the 
dyadic  scattering  amplitude  can  be  written  as 

*  '  fhh^-i  +  fhv^i  +  fvh4«  +  <2-3-J0> 


where  the  components  f  are  given  by 


f pq  “  Es  *  -  *  '  p,qe{h,v)  (2-3-21) 

Thus  if  the  incident  polarization  is  horizontal,  <j|  -  h  ? ,  the 
scattered  far  field  will  be 


^ '  <fhh« +  fvhi;>  <2-3*22> 

The  dyadic  scattering  amplitude  f  can  be  directly  related  to 
the  Fourier  transform  of  the  transition  operator.  This  relation¬ 
ship  is  derived  in  Appendix  A  of  Lang  [57].  It  is 

f  (o,i)  -  2tt2  (I-o  o)  •  t  (keo,k0i)  •  (I-i  i)  (2-3-23) 

where 

t  (k#k  )  «  — fdxdx't(xrx')ej  (2-3-24) 

(2tt)  J 

The  relationship  between  f  and  t  as  given  in  Equation  (2-3-23) 
i*s  in  dyadic  form.  By  using  the  polarization  unit  vectors  already 
defined,  the  unit  dyadic  can  be  decomposed  as  follows 

I  -  h|hj  +  vjv«  +  oo 

or  (2-3-25) 

I  *  h|h?  +  +  i  i 


This  is  a  result  of  the  fact  that  I  can  be  decomposed  in  the  above 
form  by  any  system  of  mutually  orthogonal  unit  vectors.  The  scalar 
form  of  the  relationship  can  be  found  by  dotting  Equation  (2-3-22) 
from  the  right  by  h?  or  v?  and  by  dotting  from  the  left  by  h|  or 
v*  and  then  using  Equation  (2-3-25) .  The  result  is 

f__ (o,i)  *  2ir2t  (k0o,k0i)  ,  p,qe{h,v)  (2-3-26) 

pq  —  —  pq  ~  — 

Where  *pq  "  * 

Both  the  transition  operator  and  scattering  amplitude  forms  of 
characterizing  the  scatterer  are  used.  Usually  the  transition 


operator  is  employed  in  the  derivation  of  approximate  equations, 
while  the  scattering  amplitude  appears  when  these  equations  are 
solved. 

2.3.3  Mean -Wave  Equation 

An  approximate  equation  for  the  mean  field  will  now  be 
given.  The  method  for  obtaining  this  mean  field  equation  exactly 
parallels  the  two-dimensional  procedure  as  given  in  Section  2.2.3. 
As  a  result,  the  derivation  will  not  be  given  here. 

Starting  with  the  vector  wave  equation  as  given  by  Equation 
(2-3-8),  introducing  the  dyadic  transition  operator,  averaging 
and  using  the  Foldy-Lax  approximation,  the  three-dimensional 
vector  mean -field  equation  is  found  to  be 

£<¥>  -  |  dsp (s)T(s)  •  <t>  «  £  (2-3-27) 

v 

where  £(s)  is  the  orientational  average  for  a  particle  located  at 
s,  i.e. 

7r  2ir 

T(s)  -  |  dej  d$p(0,$|s)T(s,6,$)  (2-3-28) 

0  0 

where  p(0,$)  is  the  angular  probability  density  function  given  the 
particle  is  located  at  s. 

As  before,  it  will  be  of  use  to  write  the  mean-field  equation 
explicitly.  Equation  (2-3-27)  becomes 

VxVx<E(x)>  -  k|eb(x)<E(x)  > 

-  j  ds(s)pdx,t(x~s,x*-s*) *<E(x* )>  *  -jwy0J(x)  (2-3-25) 

v 

Equation  (2-3-29)  differs  from  the  vector  Helmholtz  equation  only 
because  of  the  added  integral  term  arising  as  a  result  of  the 
scatterers.  A  similar**  scatterer-related  integral  term  is  also 
observed  in  the  case  of  two-dimensional  scatterers  [Equation 
(2-2-41) ] . 
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The  mean-field  equation  can  be  substantially  simplified  in 
the  low-frequency  (Rayleigh)  regime.  This  is  easily  shown  by 
substituting  Equation  (2-3-17)  into  Equation  (2-3-29)  and  obtain¬ 
ing 

VxVx<E>  -  kple^fx/l  +  paJ*<E>  =  -ju)y0J(x)  (2-3-30) 

In  this  low-frequency  (Rayleigh)  regime  where  the  wavelength  is 
large  relative  to  the  size  of  the  scatterers,  it  is  apparent  from 
Equation  (2-3-30)  that  solutions  to  the  mean-field  equation  can 
be  obtained  using  classical  techniques. 

2.4  Three-Dimensional  Scattering  in  Two-Dimensional  Media 

In  this  section,  mean-wave  propagation  in  a  medium  of  two- 
dimensional  scatterers  with  three-dimensional  sources  £(x)  is 
considered.  The  equation  for  the  mean  field  is  derived  from  the 
three-dimensional  mean-field  equation  given  in  Equation  (2-3-29). 
The  scattering  being  considered  differs  from  the  two-dimensional 
waves  considered  in  Section  2.2  because  it  doesn't  reduce  simply 
to  a  scalar  treatment;  in  general  the  complete  vector  problem 
must  be  considered.  This  results  from  the  fact  that  three- 
dimensional  sources  produce  waves  that  travel  at  oblique  angles 
with  respect  to  the  scatterers. 

The  scatterers  are  characterized  by  the  dyadic  transition 
operator  t(x,x' )  as  defined  in  Equation  (2-3-13).  Since  the 
scatterers  do  not  vary  with  z,  translations  of  the  incident  field 
and  of  the  equivalent  or  induced  sources  by  the  same  amount  in 
any  direction  will  have  no  effect  on  the  transition  operator; 
thus  t  (xt,x^,z,z ' )  is  just  a  function  of  z-z'  or 

t(x,x')  *  t (xt,x£,z  -  z' )  (2-4-1) 

This  is  analogous  to  the  impulse  response  of  a  time-invariant 
filter;  one  obtains  the  same  response  independently  of  when  the 
impulse  was  applied. 

Now  using  this  two-dimensional  form  of  t  in  Equation 
(2-3-39) , 
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(2-4-2) 


VxVx<E(xt/z)>  -  k2<S(xt,z)> 


-  |  d»t|dszp (s) J^x^dz '^(xt“8t ,x^-st,z-z ' ) 


*  <E(xt,z)>  =  2(xt#2) 

Here  it  is  assumed  Ejj(x)  “  1,  and  for  all  x  that  the  %otal  volume  V 
under  consideration  is  an  infinite  z-directed  cylinder  having 
cross-section  S.  An  examination  of  this  mean  equation  shown  that 
it  is  Fourier-transformable  in  the  z  direction.  Introducing  the 
Fourier  representation. 


>*  is|dkz<^(-t#kz)  >e 


-jV 


into  Equation  (2-4-2)  yields 


$x$x<E(xt,kz)>  -  k2<|(xt,kz)> 

-  |  dstp(st)|dx^£(xt-st,x^.-stfkz) 


(2-4-3) 


*  <®(2£t,k*)>  *  a(xt,k2) 


(2-4-4) 


where  V =  V.  -  jk„z° 
t  z— 


Jt  *  xB3/3x  +  £°3/3y 


£^xt#x^,kz)  *  jdzt (xt,x£,z)i 


jV 


(2-4-5) 


The  hat  notation  has  been  used  to  indicate  a  Fourier  transform 
with  respect  to  z.  The  transverse  density  p(st)  is  related  to 
p(s)  as  follows 


Jdszp(s)  *  Njdszp(s)  -  Np(st)  =  P(£t) 


(2-4-6) 


Thus  p(s^)  is  measured  in  units  of  particles  per  area  rather  them 
particles  per  volume. 

As  in  the  case  of  scalar  two-dimensional  fields  and  vector 
three-dimensional  fields,  the  fields  scattered  by  an  isolated 
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two-dimensional  scatterer  can  be  characterized  by  a  scatt  ring 
amplitude.  The  difference  between  this  vector  two-dimens  onal 
scattering  amplitude  and  the  scalar  case  is  that  the  direction 
vector  i  of  the  incident  plane  wave  of  unit  amplitude  and  polar¬ 
ization  2°  is  not  necessarily  perpendicular  to  the  generating 
element  of  the  cylindrical  scatterer  (z-axis).  The  dyadic 
scattering  amplitude  is  defined  in  terms  of  the  transformed 
scattered  field  as  follows 

‘jktxt“jkz2 

es(x.,z)  'v  f  (o,i)*qa  - -  (2-4-7) 

/5t 

where 

i  =  i t  +  iz£°  #  i2  «  i  •  z°  (2-4-S) 

and 

-  "  -t  +  iz-°  (2-4-9) 

The  wavenumber  in  the  z  direction  is  k  *k0i„.  This  is  deter- 

z  z 

mined  wholly  by  the  incident  wave.  Because  of  the  two-dimen¬ 
sional  property  of  the  scatterer,  the  wave  number  k_  is  preserved 
by  the  scattering  process;  thus  the  scattered  wave  in  direction 
o  must  also  have  a  z  component  i  .  This  means  that  scattering  is 
restricted  to  lie  on  a  cone  of  angle  as  shown  in  Figure  2-4. 

Since  both  f  and  t  describe  the  single-scattering  process, 
they  should  be  related  in  a  similar  manner  to  the  f  and  t  of  the 
three-dimensional  scatterers.  In  Appendix  B,  it  is  shown  that 
this  is  indeed  the  case.  The  result  is 

f  (o,i;k2)  *  y  (I-oo)  •  t  (kt£t  ,ktit  ;kz)  •  (I-i  i)  (2-4-10) 

where  y  is  given  by  Equation  (2-2-24)  and 

x  if-  j (k.»x.  -  k' *x') 

t(k  ,k' ;k  )  =  — — — -  dx.  dx* t (x. ,x* ;k  )e  *  (2-4-11) 

-  -t  -t  Z  (27r)  2  j  -t  ~t-  -t  Z 


i/a 


2-24 


3.0  Dyadic  Scattering  Amplitudes 

The  theory  of  electromagnetic  scattering  presented  in  Section 
2  assumed  the  knowledge  or  availability  of  the  dyadic  transition 
operators  (£)  or  dyadic  scattering  amplitudes  (f ) .  It  is  the  pur¬ 
pose  of  this  section  to  provide  the  dyadic  scattering  amplitudes 
for  trunks ,  branches  and  leaves. 

For  both  two-  and  three-dimensional  objects ,  the  scattering 
amplitude  is  determined  by  illuminating  the  object  with  a  unit- 
amplitude  plane  wave 

E^x)  ■  a°e"jk*(-‘-)  (3-0-1) 


where  i  is  a  unit  vector  in  the  direction  of  propagation  and  <£° 
i»  the  unit  polarization  vector.  The  -i  is  often  specified  by 
the  spherical  angles  0^  and  ^  measured  with  respect  to  a  polar 
axis  z.  One  can  then  write: 

-i  ■  x#sin0i  cos^  +  £#sin0i  sin^  +  zccos0^  (3-0-2) 

Incident  waves  with  both  horizontal,  SL#*h®,  and  vertical  g*»v|, 
polarizations  are  considered.  Here  the  horizontal  polarization 
vector  h|  is  taken  parallel  to  the  x-y  plane  of  the  forest  floor. 
More  specifically, 

i  x  z# 

h?  -  -  (3-0-3) 

-i  I  i  x  z  •  | 

The  vertical  polarization  vector  is  taken  perpendicular  to  both 
i  and  h?;  thus 

vj  -  h[  x  i  (3-0-4) 


For  three-dimensional  scatterers  the  dyadic  scattering  am¬ 
plitude  is  defined  by 


-jk, (o  *  x) 

E^x)  %  f  (o,i)  -2°  - - ~ - 


(3-0-5) 


where  o  is  a  unit  vector  in  the  direction  of  the  observation 
point  x.  The  o  vector  can  be  described  by  angles  0  and  $ 
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measured  with  respect  to  the  z  and  x  axes,  respectively.  In  terms 
of  these  angles,  o  is  written  as 

o  *  x°sin0g  cos$#  +  y°sin0g  sin$g  +  z#cos0g  (3-0-6) 

The  scattered  wave  in  general  can  have  components  in  the  horizontal 
h*s  and  the  vertical  v|  directions.  These  scattered  polarization 
vectors  are  defined  by 


o  x  z* 
h*  .  _Z - 

Is  x  *°l 

(3-0-7) 

v?  »  h?  x  o 
— b  — e  — 


As  pointed  out  in  Section  2,  f (o,i)  has  only  four  components 
since  o*f(o,i)  *0  and  f(o#i)»i.*0.  Thus, 

£<2-i>  -  ♦  fvhyxi  ♦  *vhv;y  +  fwv;vj  0-0-8) 

The  dyadic  scattering  amplitude  for  two-dimensional  scatterers 
is  identified  as  the  coefficient  of  a  cylindrically-expanding  wave: 

-jk0 (o«x) 

E  (x)  *  f  (o, i)  *q#  ? - =-=-  (3-0-9) 

-®  -  -  /* 

where  0g  ■  tt-0^.  This  restriction  on  the  angle  0g  results  from 
the  requirement  that  scattered  fields  have  the  same  variation  in 
a  as  the  incident  wave.  This  means  that  o  must  lie  on  the  sur¬ 
face  of  a  cone  of  angle  0^. 


Infinitely-Long,  Circular  Cylinder  (Tree  Trunk) 


In  this  section,  the  dyadic  scattering  amplitude  for  an 
infinite  circular  cylinder  of  radius  a  and  dielectric  constant 
e£  is  given.  Assume  that  the  axis  of  the  cylinder  is  the  polar 
z  axis  [Refer  to  Figure  3-1].  There  is  no  need  to  find  the 
scattering  amplitude  from  first  principles  since  the  required 
results  are  given  by  Wait  [82]  and  summarized  by  Ruck  et  al. 
[78].  These  results  are: 


3-2 


IT  ej  n/4  r  n  dd  ( 

, (o»i)  «J-  _e-- - l  (-l)nCppe  1  8  ,  pe{h,v)  (3-1-1) 

1 7T  /k0sine^  n*-«° 


1 7f  /kosin0^  n**-® 


where 


n_ 

V 


Vn-qnJn(x0’Hnal(x0)Joixl> 

PnVl<JnH,5!)(x0)Jn(x1)]1 


(3-1-2) 


(3-1-3) 


MnNn-gnJn^0|Hn)(x0)Jn(xl) 

pnNn-',5nHnI,(x0>Jn(xl>)2 


C  -zil 

n  irx„ 


F _ 80^n(xl) _ 1 

Lp„Vl%Hn'’(xO>Jn<xl>l2J 


(3-1-4) 


(3-1-5) 


xn  -  k0a 


(3-1-6) 


*  k0a/ei-cos2ei 


(3-1-7) 


=  _ i_l 

n  k0a  Le£"CO820^  sin20^J 


(3-1-8) 


»lJn(~  A(xl>  '  Vn'VW 


(3-1-9) 


rlHn2)(x0,JA,xl»  ‘  Vn^'W*!’  0-1-10) 
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,¥r 


(3-1-11) 


Nn  =  ■1H,;,)(x0)ji(x1)  -  .^'(vvv 


Mn  "  rlJn<xO,JA(xl)  *  VA'VW 


(3-1-12) 


50  ~  sin 6^ 


/e^-cos26^ 


(3-1-13) 


rl  * 


(3-1-14) 


'eJt-cos20i 


with  Jn(*)  being  the  Bessel  function  of  order  n  and  H^2^  (•)  being 
the  Hankel  function  of  the  second  kind  of  order  n.  The  primes 
over  the  cylinder  functions  mean  derivatives  with  respect  to  their 
arguments. 

An  examination  of  the  results  shows  that ,  in  general ,  scatter¬ 
ing  from  the  cylinder  gives  rise  to  depolarization  since  fhv*  °* 
It  should  be  noted  that  this  depolarization  does  not  exist  when 
the  incident  wave  is  normal  to  the  cylinder  (6^  *  ir/2)  because 
then  qn«0.  The  cross-polarized  terms  are  also  zero  in  the  for¬ 
ward-scattered  direction.  In  this  case,  9  »6.+ir  which  leads  to 

i _ \  # _ x  HI 


<£*>  «-C(vh)  (n  *  1,2,...)  and  Cn*0.  Thus, 
n  — n  u 

Wi'i>  -  -fvh(i'i>  '  0 


(3-1-15) 


3.2  Finitely-Long,  Circular  Cylinder  (Branch) 


Consider  a  planar  electromagnetic  wave  to  be  incident  upon  a 
finitely-long,  circular  cylinder  of  complex  relative  permittivity 
.  without  loss  of  generality,  the  longitudinal  axis  of  the 
cylinder  can  be  taken  as  inclined  to  the  z-axis  by  an  angle  9, 
and  the  plane  defined  by  the  cylinder  axis  and  the  z-axis  as  ro¬ 
tated  an  angle  $  from  the  x-z  plane  (refer  to  Figure  3-2).  An 
incident  electric  field  of  unit  amplitude  can  be  represented  as 
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(3-2-1) 


where 


^(x)  =  a#e“^ko*-*-* 


kc  »  u>/n0e, 


(3-2-2) 


-i  =  x’cos^  sinS^  +  y®sin$^  sin0^  +  2°cos0^ 


(3-2-3) 


x  »  x°x  +  £#y  +  z°z 


(3-2-4 ) 


and  0^  is  the  angle  between  the  direction  of  the  incident  wave, 
and  the  z-axis,  and  is  the  angle  between  the  plane  defined 
by  i  and  the  z-axis  and  the  x-z  plane. 

If  the  cylinder  is  sufficiently  long,  and  thin  in  comparison 
to  the  wavelength  within,  the  induced  electric  field  within 
the  cylinder  may  be  ascertained  using  quasi-static  techniques. 
Under  this  approximation,  the  electromagnetic  boundary  condition 
requiring  the  continuity  of  the  tangential  field  components  (E 
and  H)  can  be  employed  to  show  that  the  induced  electric  field 
component  directed  along  the  axis  of  the  cylinder  is  given  by 


(3-2-5) 


whereas  the  induced  electric  field  component  normal  to  the  axis 
of  the  cylinder  is  given  by 


r°  x  E^  =  r°  x 


eAh^) 


(3-2-6) 


where  r°,  a  unit  vector  directed  along  the  axis  of  the  cylinder, 
is  given  by 


r°  =  x®cos$  sin0  £°sin$  sin0  +  z°cos0 


(3-2-7) 


As  a  consequence  of  Equations  (3-2-1),  (3-2-5),  and  (3-2-6),  the 
induced  electric  field  within  the  cylinder  may  be  shown  to  be 


ea(x)  -  (a* •£•)£•-♦  ^  +  r^a8-r»)r»i 


(3-2-8) 


As  a  consequence  of  the  radiation  condition,  the  electric 
field  E  in  free-space  attributable  to  the  current  distribution  J 
within  the  closed  volume  is  given  by 


^  _■  m _ ,  m  •  -  j  - 


(3-2-9) 


E (x)  *  -jwy0J  dx'G(x,x' ) • J(x' ) 


where  G(x,x')  is  the  free-space  dyadic  Green's  function  [126]. 
Although  the  total  current  induced  within  the  finite  cylinder 

is 


J  *  j(oe0e.E, 


t- i 

that  part  of  the  total  current  responsible  for  the  scattered 
field  Eg  is  only 


(3-2-10) 


Jg  *  ju>e0  (ej-DEj,  »  jweoX0E 


l-i 


(3-2-11) 


where  X£  is  the  susceptibility  of  the  cylinder.  Substitution  of 
Equation  (3-2-11)  into  Equation  (3-2-9)  provides  the  following 
relation  between  the  field  E_,  scattered  by  the  cylinder  and  the 

*■5 

field  E£  induced  within  the  cylinder 


Eg (x)  =  kSx£|  dx'G(x,x') *E£(x') 


(3-2-12) 


The  far-field  free-space  dyadic  Green's  function  is  given 


by 


G(x,x' )  *  (I  - o  o) 


k-jk0 [o* (x-x') ] 
~  . .  4ttx  ~  ~ 


(3-2-13) 


with 


o  *  x°cos4>s  sin9g  +  y°sin$a  sin0g  +  z°cos6 


s  s 


s 


(3-2-14) 


Substitution  of  Equations  (3-2-8)  and  (3-2-13)  into  Equation 
(3-2-12)  yields 


§s(x)  =  f  •  a* 


e 


-jk0x 


-S  X 

where  f,  the  dyadic  scattering  amplitude,  is  given  by 
f(o,i)  -S(I-OO).  [r-r-  ♦  (jjVt) <* -  £*£* >J  ~ -jgf  1 


(3-2-15) 


(3-2-16) 


c  “ 


-4 


and 

kL  -  kft(o  -  i)  •  r*  (3-2-17) 

p  «  (irX£)  (*a2)  (k./27r)2  (3-2-18) 

3.3  Circular  Disc  (Leaf) 

Consider  a  planar  electromagnetic  wave  to  be  incident  upon  a 
circular  disc  of  radius  a,  thickness  t,  and  complex  relative  per¬ 
mittivity  e^.  Without  loss  of  generality*  the  orientation  of  the 
disc  can  be  defined  in  terms  of  the  two  eulerian  angles  6  and  4 
shown  in  Figure  3-3.  An  incident  electric  field  of  unit  amplitude 
can  be  represented  as 

E^x)  »  aoe-3k0(i.x)  (3-3-1) 

where 

k0  «  w/Upe0  (3-3-2) 

-i  *  x'cos^k  sin0i  +  £°sin$^  sin0i  +  z°cos9i  (3  3-3) 

X  «  x°x  +  £°y  +  z°z  (3-3-4) 

and  6^  is  the  angle  between  the  direction  of  the  incident  wave  and 
the  z-axis,  and  4^  is  the  angle  between  the  x-z  plane  and  the 
plane  defined  by  i  and  the  z-axis. 

If  the  disc  is  relatively  thin  (a>>t)  and  the  radius  large  in 
comparison  with  the  wavelength  (a>>A) *  the  induced  electric  field 
E£  within  the  disc  may  be  approximated  by  th<»  electric  field  in 
an  unbounded  slab  having  the  same  orientation  as  the  disc.  Under 
this  approximation,  the  electromagnetic  boundary  condition  requir¬ 
ing  the  continuity  of  the  tangential  field  components  (E  and  H) 
across  an  arbitrary  interface  can  be  employed  to  show  that  the 
induced  electric  field  within  the  disc  is  given,  approximately,  by 

-jk0  (i*.**) 

E4(x)  *  E (n)e 


(3-3-5) 


(3-3-6) 


IT-" 


mr^s\ 


where , 


and 


E(n)  =  9+eqe^Kn  +  ale“e“jlcn 


n°  =  x°sin$  sin6  -  y®cos$  sin6  +  z®cos0 


-n  =  n°  •  x 


i t  =  i  -  (n°  •  i°)n° 


<  =  k0  (e.  -  i^)' 


litl 


(3-3-7) 

(3-3-8) 

(3-3-9) 

(3-3-10) 


The  unit  polarization  vectors  appearing  in  Equation  (3-3-6)  depend 
upon  the  relative  orientation  of  the  disc.  As  a  consequence 


<i±  *  «#) 

1 

h: 

1 

I1!  X  — °  I 

II 

oV 

* 

V® 

|  1 

h?  x  i° 

—X 

1 

_3  — _ 

and 


i±  =  lit  -  -°(e£  “  it)**3/y^I 


(3-3-li;t 

(3-3-12) 

(3-3-13) 


The  amplitudes  of  the  induced  field  components,  e^,  are  given  by 


t  G-jM»c-k0)t/2 

+  _  q _ 

q"  1  -  rie^51^ 


e’  =  r  eVjKt 

q  q  q 


(3-3-14) 


(3-3-15) 


where 


J  rv 


(Ej-q)15  -  d-i£) ** 


(3-3-16) 


(3-3-17) 
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1 

•A 


■  A 


2(l-i£) 


(e^i*)’  +  (1-i*)* 

2/£7  »-*{>* 

+  e.d-i*)1 


(3-3-18) 


(2-3-19) 


As  a  consequence  of  the  radiation  condition,  the  electric  field 
scattered  by  the  disc  and  the  electric  field  induced  within 
the  disc  E£  are  related  by  the  Kirchhoff  integral 

E^x)  «  k^XjJ  dx'G(x,x')*E£(x')  (3-3-20) 


where  G(x,x')  is  the  free-space  dyadic  Green's  function  and  X£ 
is  the  susceptibility  of  the  disc  (X£  *  e-  -  1) .  The  far-field 
free-space  dyadic  Green's  function,  needed  to  find  the  dyadic 
scattering  amplitude  f  which  has  been  defined  so  that 


etx)  «  t  •  a° 


-jk0x 


is  given  by 


G(x,x')  -  (I-oo) 


,-jk0 [o* (x-x* ) ] 


(3-3-21) 


(3-3-22) 


If  Equation  (3-3-20)  is  generalized  to  accommodate  a  dyadic  field 
representation  so  that 


£g(x)  -  k.XjJ  dx'G(x,x')  *E£(x') 


(3-3-23) 


then  it  is  apparent  from  Equations  (3-3-21),  (3-3-22),  and  (3-3-23) 
that 

f  (o,i)  -  (itX£)  (k0/2ir)  2  j  dx*  (I  -  o  o)  •  E£  (x’  )ejk®  ]  (3-3-24) 


and  where,  as  a  consequence  of  Equation  (3-3-5) 

jj,  (x)  -  +  eVJKnJe-jk«(if5> 


(3-3-25) 


3-12 


with 


e*  *  ©wh?h#  +  ef.v?v° 
™  n — i—  v — x — 


(3-3-26) 


Substituting  Equation  (3-3-25)  into  Equation  (3-3-24)  and  intro¬ 
ducing  the  identity 


o  *  ot  +  (o*n°)n° 


yields 


f (o,i)  «  (irXA)  (ke/2ir)*J  dx'(l-oo)*  [e+ej 


*n  +  e-e'lKnl 


(3-3-27) 


where 


iit  ■  9.t  -  it 


(3-3-28) 


(3-3-29) 


For  a  circular  disc  of  radius  a  and  thickness  t ,  Equation 
(3-3-28)  can  be  evaluated  directly  to  obtain 


f(o,i)  -  (tXA)  (k0/2Tr)2((a/vt.)J1(2Travt)] 

*  (£“22)  *  Ie+cinc0+  +  e“sinc0”] 
where  ( • )  is  the  Bessel  function  of  order  unity  and 

e*  -  (k,t/2)  [(o-n'J  +  Uj-i*)*] 

|vt|  *  vt  =  (k0/2ir)(£2  +  n2)1* 

£  “cos0[sin0.  sin ($-$.)  +  sin0_  sin ($-#_)  J 

1  X  S  8 

-  sir.0[cos@.  +  cos0_] 
x  s 

ri  »  sin0^  cos(4>-$^)  +  sin0g  cos($-$g) 


(3-3-30) 


(3-3-31) 


(3-3-33) 


(3-3-34) 


The  dyadic  scattering  amplitude  of  the  circular  disc  as 
given  by  Equation  (3-3-30)  simplifies  substantially  in  the  case 
of  an  electrically  thin  disc  where 
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k0/e^t  «  1 


(3-3-35) 


In  this  case,  the  phase  angles  xt  and  6  are  small  and  Equation 
(3-3-30)  can  be  approximated  by 

f(o,i)*a®  *  (ffX£)  (h./2iT)  2tl  (a/vt)  Jx  (2Travt)  ] 

•  (I -oo)  •  (2*  -  !—-(!»•. 2*)n*l  (3-3-36) 

SL 

For  this  thin  disc  case  in  the  low-frequency  (Rayleigh)  regime 
where,  in  addition  to  inequality  (3-3-35), 

a  «  X  (3-3-37) 

then  avfc  =  0  and 

X# 

f(o,i).2°  *  X£t(M/2)J(I-oo)‘  [2°  -  (n#-3#)n*l  (3-3-38) 

I 

Note  that  in  the  direction  of  forward  scatter  where  o  «  i, 

6s  “  ff  “  0i  '  *  *1  +  *  (3-3-39) 

*  0,  and  Equation  (3-3-36)  for  the  electrically  thin  disc 
reduces  to  the  Rayleigh  result  of  Equation  (3-3-38). 
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Coherent  Forest  Scatterin' 


A  physically  appealing  representation  for  the  mean  field  can 
be  obtained  by  postulating  that  the  mean  (or  coherent)  field  sat¬ 
isfies  Maxwell's  equation  "in  the  mean”  and  that  the  ensemble  of 
discrete  scatterers  can  be  replaced  by  an  equivalent  continuous 
medium  described  by  an  effective  dyadic  permittivity  t,  or  alter¬ 
natively  ,  by  an  effective  dyadic  susceptibility  Because  £  is 
found  to  be  directly  proportional  to  the  fractional  volume  occupied 
by  the  scatterers/  it  is  easily  scaled  with  respect  to  the  forest 
density  and  therefore  preferable  to  e  for  characterization  of  the 
equivalent  continuous  medium.  In  Section  4.1,  general  expressions 
are  derived  which  relate  £  directly  to  the  dyadic  scattering  am¬ 
plitudes  f.  In  subsequent  sub-sections,  these  expressions  are 
employed  to  determine  for  tree  trunks,  for  branches  and  leaves, 
specific  expressions  for  their  respective  effective  dyadic  suscep¬ 
tibilities.  These  specific  expressions  are  exemplified  by  calcu¬ 
lations  which  are  then  compared  with  experiment.  It  may  be  noted 
that,  as  all  forest  constituents  are  assumed  uniformly  distributed 
about  the  vertical,  the  equivalent  continuous  medium  is  uniaxially 
anisotropic . 

Plane-wave  propagation  within  an  unbounded,  equivalent  con¬ 
tinuous  medium  is  considered  in  Section  4.2.  Here,  a  general 
dispersion  relation  is  derived  relating  the  plane-wave  propagation 
constants  ic  of  the  mean  (or  coherent)  field  and  the  dyadic  scat¬ 
tering  amplitudes  f.  In  subsequent  sub-sections,  this  dispersion 
relation  is  employed  to  determine,  for  plane-wave  propagation 
through  an  unbounded  forest  of  tree  trunks,  branches,  or  leaves, 
specific  expressions  for  the  wave-propagation  constant.  The  imag¬ 
inary  part  of  the  wave -propagation  constant,  the  specific  attenua¬ 
tion,  is  numerically  evaluated  and  compared  with  experiment. 

Experimental  verification  of  the  electromagnetic  model  and, 
ultimately,  prediction  of  radio  system  performance  requires  the 
identification  of  measurable  quantitative  parameters  to  charac¬ 
terize  the  forest.  In  Section  4.3,  some  of  these  parameters  are 
identified  and  related  directly  to  the  specific  attenuation. 
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4.1  Effective  Dyadic  Permittivity  of  the  Forest 

The  mean  electric  field  propagating  within  an  unbounded  en¬ 
semble  of  three-dimensional  scatterers  can  be  determined ,  as  shown 
in  Section  2.3,  from  the  mean  wave -equation 

L*<E(x)>  -  Jdsp (s) Jdx' t (X”8,x'-s) •<E(x* ) >  *  -jwy0J(x)  (4-1-1) 


where 


L  *  VxVxI  -  k|l 


(4-1-2) 


From  the  physical  viewpoint,  a  more  appealing  representation  of  the 
mean  field  can  be  obtained  by  postulating  that  the  mean  field  sat¬ 
isfies  Maxwell's  equations  in-the-mean  and  that  the  ensemble  of 
discrete  scatterers  can  be  replaced  by  an  equivalent  continuous 
medium  described  by  an  effective  dyadic  permittivity  e.  If 
Faraday's  law 


Vx<E(x)>  *  -jmy0<H(x)>  i 

jls  substituted  into  the  first  term  of  Equation  (4-1-1) ,  then 

Vx<H (x) >  *  jwe0 (I*<E(x) >  +  Jdsp (s) Jdx' t(x-s,x'-s) 


(4-1-3) 


•  <E(x')>  +  J(x)  ] 
and,  by  Ampere's  law 


Vx<H(x)>  •  jw<D(x)>  +  J  (x) 


(4-1-4) 


(4-1-5) 


it  follows  that 


<D(x)>  »  e0[I»<E(x)>  +  jdsp (b)  jdx' t (x-s  ,x'-b) 


•  <E(x')>] 


(4-1-6) 


Taking  the  Fourier  transform  with  respect  to  x  transforms  Equation 
(4-1-6)  into 

<D(k)>  -  e*(I-<|(k)>  jd£P(£)  j  dx'£(k,x'-s)  •  <E(x' )  >ej-‘-  (4-1-7) 


4-2 


The  second  integral  appearing  in  Equation  (4-1-7)  may  be  recog¬ 
nized  as  a  convolution  with  respect  to  the  parameter  s;  as  a 
consequence 


<D(k)>  *  e0[l*<E(k)>  +  [  dsp(s)  |dk't(k,k') 

k0  J  ~  J 


•  <E(k')>ej 


(4-1-8) 


If  the  number  density  p(s)  is  independent  of  location  so  that 

p(s)  -  p  (4-1-9) 

then  the  order  of  integration  in  Equation  (4-1-8)  may  be  inverted 
with  the  result 


<D(k)> 


e.[J  +  -  -,-p  £(k,k)]-<i(k)> 
k“ 


The  constitutive  relation 


(4-1-10) 


D  *  ec£  *  E  (4-1-11) 

indicates  that  Equation  (4-1-10)  may  be  re-written  in  the  form 

<D(k)>  -  ect  •  <E(k)>  (4-1-12) 

and  that 


£  *  I  +  t(k,k) 


(4-1-13) 


is  the  effective  dyadic  permittivity  of  the  equivalent  continuous 
medium.  Note  that  £  depends  upon  the  direction  of  propagation 
through  the  wave  vector  k;  such  a  medium  is  termed  spatially  dis- 
persive*.  Because  t  is  a  dyadic,  the  medium  is  termed  anisotropic. 
It  is  also  interesting  to  note  that,  far  from  the  source  where  the 
mean  field  <E(x)>  is  essentially  planar,  only  those  components  of 
£  which  are  orthogonal  to  the  direction  of  propagation  will  be  of 
significance.  This  is  consistent  with  the  relation  established  in 
Section  2.3  between  t  and  the  dyadic  scattering  amplitude  of  f 
[refer  to  Equation  (2-3-23)]. 

In  the  low-frequency  (Rayleigh  scattering)  limit,  where  each 
scatterer  can  be  represented  as  an  electric  dipole  with  dyadic 
polarization  a,  the  Fourier  transform  of  Equation  (2-3-17) 


*Landau  and  Li fs chit*. 
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(4-1-14) 


£(k,k')  -  k*a/(2n)* 


can  be  substituted  into  Equation  (4-1-15)  to  yield 

e  -  l  +  pa  (4-1-15) 

in  the  Rayleigh  regime.  Note  that  in  the  Rayleigh  regime, 

£  is  no  longer  dependent  upon  the  direction  of  the  wave  vector 
k;  however,  the  medium  can  be  anisotropic. 

The  effective  dyadic  susceptibility  can  be  defined  as 

X  •  e  -  I  (4-1-16) 

so  that,  in  general,  as  a  consequence  of  Equation  (4-1-13) 


(4-1-17) 


However,  in  the  Rayleigh  regime,  as  a  consequence  of  Equation  (4-1* 

X  “  pa  (4-1-18) 

It  is  apparent  from  Equations  (4-1-17)  and  (4-1-18)  that  the  ef¬ 
fective  dyadic  susceptibility  is  directly  proportional  to  the  num¬ 
ber  density  of  the  discrete  scatterers.  For  this  reason,  x  is 
often  preferred  over  t  for  the  characterization  of  the  equivalent 
medium.  Because  of  the  assumed  exp{jut>  time-dependence  of  all 
field  quantities,  it  proves  convenient  to  define  the  real  and 
imaginary  parts  of  x  »o  that  x  *  X~}X* 

In  general,  the  effective  dyadic  sysceptibility  x  can  be 
written  explicitly  in  terms  of  the  unit  vectors  h® ,  v°,  and  i®. 
Here,  i°  is  a  unit  vector  in  the  direction  of  propagation  and  h° 
and  v®  are  the  polarization  vectors.  In  component  form,  x 


•  l  XoBa*@' 

o®go 


o*,6e{h®,v®,i®} 


(4-1-19) 


There  are  nine  components.  However,  only  x^h,  *hv,  xvh,  and  *w 
are  important.  The  other  components  are  either  zero  or  do  not 
contribute  at  this  level  of  approximation.  These  four  principal 
components  of  the  effective  dyadic  susceptibility  can  be  deter¬ 
mined  directly  from  the  components  of  the  dyadic  scattering  ampli¬ 
tude  £.  Equations  (2-3-2G) ,  (4-1-17) ,  and  (4-1-19)  yield 


't 


p,qe{h,v) 


(4-1-20) 
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The  development  presented  above  leading  to  the  concept  of  an 
effective  dyadic  permittivity  was  predicated  upon  the  mean-wave 
equation  for  three-dimensional  scattering  from  three-dimensional 
scatterers  [Equation  (4-1-1) ]„  An  analogous  result,  the  mean- 
wave  equation  for  three-dimensional  scattering  from  two-dimensional 
scatterers,  was  developed  in  Section  2.4  [Equation  (2-3-29)].  As 
a  consequence,  an  effective  dyadic  susceptibility  for  a  continuous 
medium  equivalent  to  an  unbounded  ensemble  of  two-dimensional  scat¬ 
terers  can  be  analogously  defined  with  the  effective  dyadic  sus- 
ceptiblility  given  by 

X  *  12lAL£|(kt,kt;k2)  (4-1-21) 

k0 

Here  again  only  the  Xhh,  Xhv,  Xvtl  and  components  are  important. 
They  can  be  directly  related  to  the  four  components  of  the  two- 
dimensional  dyadic  scattering  amplitude  by  employing  Equation  (2- 
4-10}  in  component  form.  The  result  is 

xna  *  flinT IT  •  p,qe(h, v)  (4-1-22) 

***  Yko  ** 

where  y  is  given  in  Equation  (2-2-24). 


4.1.1  Effective  Dyadic  Susceptibility  of  Trunks 


The  radiowave  propagation  model  developed  in  this  report 
views  the  forest  as  a  random  ensemble  of  tree  trunks,  branches  and 
leaves  having  prescribed  location  and  orientation  statistics.  Tree 
trunks  are  modelled  as  vertical  lossy  dielectric  circular  cylinders 
of  infinite  length.  The  salient  scattering  properties  of  individual 
cylinders  can  be  characterized  in  terms  of  the  dyadic  scattering 
amplitude  £  which  was  determined  earlier  in  Section  3.1.  The 
effective  dyadic  susceptibility  for  an  unbounded  forest  of  tree 
trunks  can  be  found  using  Equation  (4-1-22). 


jp4 

(k0) 2sin0£ 


C(PP) 


) 


9 


pe{h,v> 


(4-1-1-1) 


The  C(PP)  are  given  by  Equations  (3-1-3)  and  (3-1-4).  As  was 
n 

pointed  out  in  Section  3,  there  is  no  depolarization  of  the  elec¬ 
tromagnetic  wave  in  the  forward  scattering  direction  by  vertical 
circular  cylinders,  i.e.  fhv(i.#i)  “  fhv(i,i)  *  thus 


*hv  *  Xvh 


0 


(4-1-1-2) 


and  the  equivalent  continuous  medium  for  vertical  trunks  is  uni- 
axially  anisotropic. 

The  effective  dyadic  susceptibility  of  tree  trunks  as  ex¬ 
pressed  by  Equation  (4-1-1-1)  is  exemplified  in  Figures  4-1  and 
4-2  for  the  case  of  a  radiowave  propagating  parallel  to  the  forest 
floor  (0^=90°).  These  calculations  are  based  upon  a  tree  trunk 
radius  of  10  centimeters  and  a  trunk  number  density  of  1  trunk  per 
square  meter.  Three  models  are  employed  to  describe  the  permit¬ 
tivity  of  the  wood  [refer  to  Appendix  A] ,  but  only  Cyber Com  model 
III  can  be  considered  realistic.  The  horizontal  components  of  the 
effective  dyadic  susceptibility  (denoted  by  X^,  in  general,  but 
here  for  0^=90°  by  Xt)  are  shown  in  Figure  4-1;  the  vertical  com¬ 
ponent  of  the  effective  dyadic  susceptibility  (denoted  by  Xw,  in 
general,  but  here  for  0. *90°  by  X_)  shown  in  Figure  4-2. 

It  is  seen  from  these  figures  that  the  effective  dyadic  sus¬ 
ceptibility  £  is  relatively  insensitive  to  the  Cyber Com  model 
employed  for  wood  permittivity.  Further,  although  the  resonant 
response  apparent  in  these  figures  is  not  shown  in  detail  (calcu¬ 
lations  having  been  made  at  100  MHz  intervals) ,  it  is  clear  that 
for  tree  trunks  of  this  (10  centimeter)  radius,  resonance  plays  a 
major  role  in  the  frequency  band  200  -  2000  MHz.  In  this  frequency 
band  the  real  parts  of  the  dyadic  susceptibility  (X£,X^)  decrease 
roughly  as  the  square  of  the  frequency;  the  imaginary  part  of  the 
horizontal  component  (X£)  decreases  roughly  as  the  two-thirds 
power  of  the  frequency;  and  the  imaginary  part  of  the  vertical 
component  (X")  decreases  roughly  as  the  five-thirds  power  of  the 
frequency. 

Tamir  [107],  concerned  only  with  frequencies  below  200  MHz, 
has  suggested  that  the  forest  could  be  represented  by  an  effective 
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scalar  permittivity  (s')  of  about  1.1  and  an  effective  conductivity 
(c)  of  about  10"Hmhos/meters  [refer  to  Figure  1-1] .  The  relations 

X'  =  e*  -  1  (4-1-1-3) 

X"  *  o/we0  (4— 1—1— 4 ) 

can  be  used  to  deduce  the  effective  scalar  susceptibilities  sug¬ 
gested  by  Tamir;  these  susceptibilities  are  also  plotted  in  Figure 
4-1  and  4-2,  It  is  apparent  from  these  figures  that  the  frequency- 
independent  value  of  X'  suggested  by  Tamir  differs  significantly 
from  the  inverse  frequency-squared  dependent  of  X£,  X^  found  by 
Cyber Com;  the  inverse  frequency -dependent  value  of  X"  suggested 
by  Tamir  is  about  cne  order  of  magnitude  smaller  than  the  corres¬ 
ponding  values  of  X£,  X£  found  by  Cyber Com.  Near  agreement  between 
X"  and  X£,  X£  could  be  achieved  by  decreasing  the  assumed  tree 
trunk  number  density  from  1  tree  trunk  per  square  meter  to  0.2  tree 
trunk  per  square  meter.  Although  not  shown  in  this  figure,  the 
low-frequency  (Rayleigh)  behavior  of  the  vertical  component  (X£, 

X^)  agrees  with  Tamir' s  model  for  Cyber Com  wood  permittivity  models 
II  and  III.  This  Rayleigh  behavior  is  also  anticipated  in  the 
frequency  band  200  -  2000  MHz  for  tree  trunks  (and  branches)  of 
extremely  small  radius  (less  than  1  centimeter) . 


4.1.2  Effective  Dyadic  Susceptibility  of  Branches 

Branches  are  modeled  as  lossy-dielectric  circular  cylinders 
of  radius  a  and  length  £.  The  salient  scattering  properties  of 
individual  finite-length  cylinders  can  be  characterized  in  terms  of 
their  dyadic  scattering  amplitude  f  which  is  given  in  Section  3.2. 
The  component  susceptibilities  for  an  unbounded  forest  of  branches 
can  be  found  by  employing  Equation  (4-1-20) . 

The  four  components  of  the  dyadic  scattering  amplitude  in  the 
forward  direction  o=i  can  be  obtained  from  Equation  (3-2-16)  by 
dotting  it  from  the  left  and  right  with  the  horizontal  and  vertical 
polarization  vectors.  Thus, 


,,q(i'i>  "(“r)  (e°*£°)  (£°*sl0) 

2+X«  ) 


The  mean  scattering  amplitudes  are  given  by 
-  fi t  1 2n 


-  tv  t£v 

fpq(i,i)  =  j  dej  d*p(0^)fpq(i,i) 


0  0 


(4-1-2-1) 


(4-1-2-2) 


where  p(0,<j>)  is  the  probability  density  function  of  branch  inclin¬ 
ations.  If  it  is  assumed  that  the  branches  are  distributed  uni¬ 
formly  in  the  azimuthal  *ngle  $  and  that  all  have  the  same  polar 
angle  0fa  then 

6(6-6.) 

«  /  n  r  \  _  /  i  i  ^ 


P(6,<j>)  = 


(4-1-2-3) 


Now,  using  i  [as  given  in  Equation  (3-0-2))  in  Equations  (3-0-3) 
and  (3-0-4),  the  following  expressions  can  be  obtained  for  h? 


and  v? 

A 


h?  *=  -x°sin4>-  +  jr°cos^. 

™*~  J-  "*  A  A 


(4-1-2-4) 


v?  *  -x°cos4»i  cosSj^  -  £°sin^  cos6i  +  z°6in6i 


(4-1-2-5) 


In  addition  r°,  which  is  a  unit-vector  directed  along  the  branch, 
can  be  written  as 

r"  =  x°sin4>  sin6  +  £°sin$  sine  +  z°cos0  (4-1-2-6) 

When  these  expressions  for  the  unit  vectors  h? ,  v?  and  r°  are 
used  in  Equation  (4-1-2-1)  and  the  result  averaged  over  the  prob¬ 
ability  density  function  given  in  Equation  (4-1-2-3) 


fhh (- 


i,i)  -  - — -  sin20.  + 

\  2  /  £  L2(2+XA)  b 


(4-1-2-7) 
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I 


<# 


(c?M. 


A  £ 

-  (cos20.  sin20. 

2 (2+X^)  1  D 


+  2s.in26^  cos20j3)  + 


(4-1-2-8) 


are  obtained.  The  depolarized  forward-scattering  amplitudes  both 
average  to  zero  so  that 


fhv<i'i>  *  fvh(i'i>  *  0 


(4-1-2-9) 


The  four  principal  components  of  the  effective  dyadic  suscepti¬ 
bility  are  then  found  by  substituting  the  above  results  into 
Equation  (4-1-20)  and  obtaining 


xhh  *  ’V3’1 


2 (2+Xj) 


sin20.  + 


(4-1-2-10) 


irX^pa2  Z 


.2(2+X£) 


(cos20i  sin20b  +  2sin20i  cos26b) 


(4-1-2-11) 


XHW  ■  Xvh  "  0  (4-1-2-12) 

Thus  the  equivalent  continuous  medium  for  an  azimuthally  symmetric 
distribution  of  branches  is  uniaxially  anisotropic. 

The  effective  dyadic  susceptibility  of  the  branches  as  ex¬ 
pressed  by  Equations  (4-1-2-10)  and  (4-1-2-11)  is  exemplified  in 
Figure  4-3  for  the  case  of  a  radiowave  propagating  parallel  to 
the  forest  floor  (0^»9O°).  These  calculations  are  based  upon  a 
branch  radius  of  1  centimeter,  a  branch  length  of  1  meter,  and  a 
branch  number  density  of  1  branch  per  cubic  meter.  All  branches 
are  assumed  to  be  inclined  45  degrees  with  respect  to  the  vertical, 
but  uniformly  oriented  in  azimuth.  Three  models  are  employed  to 
describe  the  permittivity  of  the  wood  [refer  to  Appendix  A] ,  but 
only  model  III  can  be  considered  realistic.  The  horizontal 
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Effective  Dyadic  Susceptibility  (Xt»x£/*x£»X 


and  vertical  components  of  the  effective  dyadic  susceptibility 
(denoted  by  Xhh  and  X^i  in  general,  but  here  for  0O  =  90®  by  Xfc 
and  X  ,  respectively)  are  shown  in  Figure  4-3.  Assumptions  im- 
plicit  in  the  development  of  the  thin-branch  scattering  model 
preclude  the  consideration  of  branches  of  significantly  greater 
radius  over  the  frequency  band  200  -  2000  MHz. 

In  the  thin-branch  frequency  regime,  the  frequency  behavior 
of  the  effective  dyadic  susceptibility  is  determined  by  the  fre¬ 
quency  behavior  of  the  susceptibility  of  a  single  scatterer.  For 
branches,  this  is  clearly  evident  from  Equations  (4-1-2-10)  and 
(4-1-2-11) ,  where  the  only  frequency-dependent  parameter  is  the 
susceptibility  of  the  wood  X^.  Further,  since  all  three  models 
for  wood  permittivity  satisfy  the  inequality  (refer  to  Appendix 
A] 

XI  <<c  X£  *  39  (4-1-2-13) 

it  is  apparent  from  Equations  (4-1-2-10)  and  (4-1-2-11)  that  for 
branches 

X£h  -  si”!eb  +  w]  (4-1-2-14) 

X^y  5  ffX"pa2££j  (cos20i  sin20fa  + 2sin20i  cos20b)+j^J  (4-1-2-15) 

These  equations  explain  why  the  frequency  behavior  of  X"  and  X" 

t  z 

shown  in  Figure  4-3  so  closely  reflects  the  frequency  behavior  of 
X£  shown  in  Figure  A-2  of  Appendix  A.  Note  that,  in  contradis¬ 
tinction  to  resonant  scattering,  scatterer  dimensions  do  not  affect 
the  frequency  behavior  of  X. 

The  effective  scalar  susceptibilities  suggested  by  Tamir  [refer 
to  Section  4.1.1]  are  also  shown  in  Figure  4-3.  These  values  appear 
to  be  about  one  order  of  magnitude  greater  than  the  computed  values 
of  X.  This  disparity  may  be  attributed  to  the  fairly  low  branch 
number  density  assumed  for  the  calculations  (1  branch  per  cubic 
meter)  and  also  to  the  fact  that  Tamir* s  value  includes,  not  only 
the  effects  of  branches,  but  all  other  vegetative  components  as 
well. 
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Leaves  are  modelled  as  lossy-dielectric  circular  discs. 

The  salient  scattering  properties  of  individual  discs  can  be 
characterized  in  terms  of  their  dyadic  scattering  amplitude  f 
which  is  given  in  Section  3.3.  The  component  susceptibilities 
for  an  unbounded  forest  of  leaves  can  be  found  by  employing 
Equation  (4-1-20) . 

According  to  Equation  (3-3-30) ,  the  four  components  of  the 
dyadic  scattering  amplitude  in  the  direction  of  forward  scatter 
are 

f^(i,i)  *  (k0/2r)  2  (e*  sinc0+  +  e~  since”)  (4-1-3-1) 

Pm  —  ~  *•  pq  pq 

where , 

e*  =  E°  ‘  *  2°  '  p,qe{h,v}  (4-1-3-2) 

F'd 

and  all  other  parameters  are  defined  in  Section  3.3.  The  mean 
scattering  amplitude,  averaged  over  leaf  orientation  and  required 
in  Equation  (4-1-20) ,  is  given  by 


fpq(i,i)  -  |7rde|2Wd^p(e,^)fpq(i,i)  (4-1-3-3) 

0  0 

where  p(6,<|>)  is  the  probability  density  function  of  the  leaf 
inclinations.  Unfortunately,  even  under  the  assumption  that 
the  leaves  are  distributed  uniformly  in  azimuth  so  that 

P<e,*)  *  p(0)  (4-1-3-4) 


the  expression 

TT  f  2TT 

dep(6)j  eWfpqU,!.)  (4-1-3-5) 

)  0 


is  difficult  to  evaluate  analytically  because  of  the  complex  de¬ 
pendence  of  fpq(i,i)  on  the  azimuthal  angle  $.  Although  it  can 
be  shown  on  the  basis  of  symmetry  that  the  cross-polarized  com¬ 
ponents  must  average  to  zero  so  that 
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fhv(±'i)  ■  *  0  (4-1-3-6) 

the  co-polarized  components  have  been  obtained  only  by  using  nu¬ 
merical  integration. 

Explicit  expressions  for  the  averaged  co-polarized  components 
can  be  found,  however,  when  the  discs  are  electrically  thin  and 
satisfy  the  inequality 

k 0/tl  t  «  1  (4-1-3-7) 

Then,  according  to  Equation  (3-3-38) 

*pp<i«i>  *  -  raq  <E*-n*> *]  (4-1-3-8) 

Using  Equations  (3-0-2),  (3-0-3)  and  (3-0-4)  to  recast  Equation 
(4-1-3-8)  in  terms  of  the  angular  variables  ,  and  substituting 
the  result  into  Equation  (4-1-3-5)  yields  the  following  explicit 
expressions  for  the  co-polarized  components 


-  (¥) 

1  - 

Xl  I 

T(T*xJT  Jlj 

(4-1-3-9) 

fw<i-±>  ■  (¥)  *** 

1  - 

571+X^T  (IlCos  eA  +  2I2sin  ei)J 

(4-1-3-10) 

where , 

f  Tf 

I  d0p(0)sin20 

1 0 

(4-1-3-11) 

*2  ' 

< 

r* 

d0p(0)cos20 

>0 

(4-1-3-12) 

The  effective  dyadic  susceptibility  for  an  unbounded  forest 
of  electrically  thin  leaves  can  be  obtained  by  substituting 
Equations  (4-1-3-6),  (4-1-3-9)  and  (4-1-3-10)  into  Equation  (4- 
1-20)  and  finding 

-  itX. pa3t |1 - — —  1.1  (4-1-3-13) 

1  l  2(1+Xt)  lJ 

Again  it  is  seen  that  the  equivalent  continuous  medium  for  an  azimu- 
thally-symmetric  distribution  of  leaves  is  uniaxially  anisotropic. 
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V  : 


irXApa2t 


[- 


2_a+Xir  f1!008*0!  +  2I2sin26j 


(4-1-3-14) 


r'1 

(4-1-3-15) 


subject  to  the  condition  k0i/F^  t  <<  1.  Typically,  leaves  are  no 
more  than  about  1  millimeter  thick;  according  to  Appendix  A,  £^=40 . 
As  a  consequence.  Equations  (4-1-3-13),  (4-1-3-14)  and  (4-1-3-15) 
should  prove  valid  over  the  entire  frequency  band  200  -  2000  HMz . 

The  effective  dyadic  susceptibility  of  leaves  in  the  thin- 
disc  approximation  [Equations  (4-1-3-13)  and  (4-1-3-14)]  is  exem¬ 
plified  in  Figure  4-4  for  the  case  of  a  radiowave  propagating 
parallel  to  the  forest  floor  (6^«90*).  These  calculations  are 
based  upon  a  leaf  radius  of  5  centimeters,  a  leaf  thickness  of  1 
millimeter,  and  a  leaf  number  density  of  200  leaves  per  cubic 
meter,  fairly  typical  values.  The  random  orientation  of  the  leaves 
is  described  by  a  probability  density  function  assumed  to  be  uni¬ 
form  in  azimuth,  and  uniform  in  elevation  angle  over  the  range 
0-30  degrees  (the  leaves  tend  to  be  horizontal) .  Three  models  are 
employed  to  describe  the  permittivity  of  the  leaves  [refer  to 
Appendix  A],  but  only  model  III  can  be  considered  realistic.  The 
horizontal  and  vertical  components  of  the  effective  dyadic  sus¬ 
ceptibility  (denoted  by  and  X^,  in  general,  but  here  for 
6±  «  90°  by  Xt  and  Xz,  respectively)  are  shown  in  Figure  4-4. 

The  effective  scalar  susceptibilities  suggested  by  Tamir 
[refer  to  Section  4.1.1]  are  also  shown  in  Figure  4-4.  The  agree¬ 
ment  between  Tamir’ s  values  and  those  for  X”  based  upon  leaf 
permittivity  models  II  and  III  for  frequencies  below  500  MHz  is 
remarkable.  At  higher  frequencies,  relaxation  losses  attributable 
to  polarization  of  the  water  molecule  [refer  to  Appendix  I]  are 
likely  to  be  appreciable,  thereby  invalidating  both  the  Tamir 
values  and  those  computed  on  the  basis  of  leaf  permittivity  model 
II. 

A  comparison  of  Figures  4-3  and  4-4  reveals  that  for  the 

leaves  >  X.  whereas  for  the  branches  X,.  <  X..  This  difference 
t  z  t  z 


i 


w  ‘.-i 


4-16 


Effective  Dyadic  Susceptibility  (x' 


4-17 


is  explained  by  noting  that,  for  the  leaf  and  branch  orientation 
distributions  assumed  in  the  calculations,  the  leaves  present 
primarily  a  horizontal  profile  to  a  propagating  radiowave  (discs 
with  axes  at  0°  to  30°  from  vertical),  whereas  the  branches 
present  primarily  a  vertical  profile  (rods  at  45°  that  appear 
vertical  when  in  plane  with  viewer).  Under  these  conditions,  a 
horizontally-polarized  radiowave  tends  to  be  more  susceptible  to 
the  effects  of  the  leaves,  whereas  a  vertically -polarized  radio¬ 
wave  tends  to  be  more  susceptible  to  the  effects  of  the  branches. 


4.2  Plane-Wave  Propagation 


4.2.1  Dispersion  Relation 

The  mean  electric  field  propagating  within  an  unbounded 
ensemble  of  three-dimensional  scatterers  can  be  determined,  as 
shown  in  Section  2.3,  from  the  mean  wave  equation 

L-cE(x)  >  -  Jdsp  (s)  Jdx't (x-s,x'-s)  *<E(x* )  >  *  -ju>y0J(x)  (4-2-1) 


where 


L  =  VxVxI  -  k*I 


(4-2-2) 


Consider  the  plane-wave 


<E(x) >  -  E  e"^-’- 


(4-2-3) 


to  be  propagating  in  a  source-free  region  where 

J(x)  =  0 


(4-2-4) 


Substitution  of  Equation  (4-2-3)  and  (4-2-4)  into  Equation  (4-2-1) 
yields,  upon  taking  the  Fourier  transformation  with  respect  to  the 
variable  x 

£X  (k  X  E)  +  jdsp  (s )  jdx't  (£.»  X ' -s )  E  -  0  (4-2-5) 


Recognizing  the  second  integral  appearing  in  Equation  (4-2-5)  as 
a  convolution  with  respect  to  the  parameter  s  and  assuming 
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the  number  density  of  the  scatterers  p(s)  to  be  independent  of 
location  so  that 


P(s)  -  P 

Equation  (4-2-5)  can  be  written  as 


icx  (jcxE)  +  k*  £l  +  (2tt)  s  t  (ic ,  <)  j  •  E  =  0 


(4-2-6) 


(4-2-7) 


Because  the  fractional  volume  occupied  by  the  scatterers  is 
assumed  small  [see  Section  4.3],  the  second  term  in  the  square 
brackets  of  Equation  (4-2-7)  is  also  small  and  perturbation  tech¬ 
niques  can  be  used  to  find  an  approximate  solution  for  E.  For 
ordering  purposes,  it  is  then  convenient  to  replace  p  in  Equation 
(4-2-7)  by  6p  and  expand  both  E  and  tc  as  power  series  in  6  so  that 


§  -  Eq  +  Ex  + 


(4-2-8) 


v 


K  *  (K0  +  +  . . . . )  i^ 


(4-2-9) 


Substituting  Equations  (4-2-8)  and  (4-2-9)  into  Equation  (4-2-7) 
and  setting  to  zero  the  coefficient  of  each  power  of  6  yields 


and 


•cjlix  (ixlg)]  +  *51^  -  0 


2*0*1  [ix  +  K0[ix  fixIi)] 


(4-2-10) 


+  koE^  +  (2tt)  3  p  t  (<0i,K0i)  •  E1  = 


(4-2-11) 


Equation  (4-2-10)  can  be  recognized  as  the  free-space  wave  equation; 
as  a  consequence 

Eq  •  i  *  0  (4-2-12) 


<0  * 


(4-2-13) 


Equation  (4-2-11)  can  be  simplified  by  expanding  the  vector  triple¬ 
cross  products  and  using  Equations  (4-2-12)  and  (4-2-13)  to  obtain 

-2k0x1E^  +  k|  (Ex  •  i) i  +  (2ir)  *  p  t  (k.i  ,k0i)  -  =  0  (4-2-14) 
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This  vector  equation  is  equivalent  to  three  linear  algebraic 
equations  and,  in  order  for  a  non-trivial  solution  to  exist,  the 
determinant  of  the  coefficients  must  be  zero.  The  vector  wave 
numbers  (propagation  constants)  satisfying  this  condition  can 
be  determined  by  expressing  Eg  as 

Eg  -  (Eg-h®)h°  +  (Eg*v0)v°  (4-2-15) 

and  substituting  Equation  (4-2-15)  into  Equation  (4-2-14)  to  obtain 
the  following  dispersion  relation  for 


-2k0x1  +  (2ir)  *  p  tj 


(2n)  3  p  thv 


(2ir)  3  p  i 


-2k*x1  +  (2ir)  3  p 


(4-2-16) 


under  the  assumption  that  (E^  •  i)  *0.  The  dispersion  relation 
can  also  be  written  in  the  form 


k|x^  -  Bko^  +  C  *  0 


where 


B  =  2irp[: 


hh  +  ^ 


(4-2-17) 


(4-2*18) 


C  *  (2irp)  2 


(4-2-19) 


and  Equation  (2-3-23)  has  been  employed  to  express  t  in  terms  of 
the  dyadic  scattering  amplitude  f .  As  a  consequence  of  Equation 
(4-2-17),  the  allowable  propagation  constants  are 


±  B  ± 


(4-2-20) 


'M?[ 


(?hh  +  ?vv)  ±  V(fhh  "  *vv)2  + 


hv  vh 


(4-2-21) 


For  classes  of  scatterers  for  which 


vh  hv 


(4-2-22) 
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(such  as  branches  and  leaves  having  azimuthally  uniform  orienta¬ 
tion  distributions) , 


(4-2-23) 


« l  5  *lh  -  <4-2*24> 

Thus,  if  6p  is  now  replaced  by  p  (or,  equivalently,  6  is  set  equal 
to  unity).  Equations  (4-2-9),  (4-2-13),  (4-2-23)  and  (4-2-24)  can 
be  used  to  show  that  to  first  order  in  the  perturbation  parameter 
6, 


Kp  “  k°  +  *pp  '  pe{h,v}  (4-2-25) 

whenever  Equation  (4-2-22)  is  satisfied. 


For  three-dimensional  wave  propagation  in  two-dimensional 
media  [refer  to  Section  2.4],  a  similar  development  can  be  pur¬ 
sued  to  show  that,  in  general. 


k  *  £t+kz— °  ,  k„  *  kocos0<  ,  *  tc^i 


— t  t-t 


(4-2-26) 


where 


-  kt  +  ^  [<£hh  +  fw>  ±  V<fhh-£w>i+4fhvfvh]  <4-2-27’ 


and 


If 


kt  =  koSine^^ 


(4-2-2  3) 


(4-2-29) 


(as  in  the  case  of  tree  trunks)  then  Equation  (4-2-27)  can  be 
simplified  to 

Ktp  “  +  Ippf5^r  •  pe<h,v)  (4-2-30) 
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4.2.2  Specific  Attenuation 


In  general ,  the  wave  propagation  constant  t<  has  both  real 
and  imaginary  components  so  that 


<_  ■  ic'  -  jx" 
P  p  p 


pe(h,v> 


(4-2-31) 


and,  as  a  consequence. 


<E(x)>  =  E  e 


(4-2-32) 


The  real  part  of  the  propagation  constant  (k^)  is  expressed  in 
radians  per  meter;  the  imaginary  part  (x*) ,  also  called  the 
specific  attenuation,  is  expressed  in  nepers  per  meter  or,  alter¬ 
natively,  in  decibels  per  meter  through  the  relation 


»  -I 


a  *  (201og1  e)x£  2  8.686  x£  (4-2-33) 

Because  the  fractional  volume  occupied  by  the  vegetation  is  small, 
the  real  part  of  the  propagation  constant  (x^)  is  dominated  by 
the  free-space  component  (k„) ;  this  is  shown  in  Figure  4-5.  The 
effect  of  the  imaginary  part  (the  specific  attenuation  x£) ,  how¬ 
ever,  can  be  appreciable  and  is  considered  further  in  the  following 
sub-sections. 

Some  investigators  (34,  35,  44,  45]  have  reported  that  the 
measured  specific  attenuation  decreases  with  increasing  path 
length.  The  theory  developed  in  this  report  predicts  no  such 
behavior  for  the  mean  field  of  a  radiowave  propagating  directly 
through  an  unbounded  forest.  Such  behavior  might  conceivably 
arise,  however,  as  a  consequence  of  lateral-wave  propagation  modes 
along  the  air- forest  interface  [refer  to  Section  5] ,  non-homogen- 
eous  transmission  paths  through  non-uni formly  forested  regions, 
and/or  random  (non -coherent)  scattered  field  behavior  [refer  to 
Section  7 ] . 

4. 2. 2.1  Specific  Attenuation  of  Tree  Trunks 

The  radiowave  propagation  model  developed  in  this  report 
views  the  forest  as  a  random  ensemble  of  tree  trunks,  branches 
and  leaves  having  prescribed  location  and  orientation  statistics. 
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Figure  4-5 i  Propagation  Constant  (Real  Part) 

[Frat-^pact ;  Tree  Trunks ,  Branches,  leaves] 


Tree  trunks  are  modelled  as  vertical  lossy  dielectric  circular 
cylinders  of  infinite  length.  The  salient  scattering  properties 
of  individual  cylinders  can  be  characterized  in  terms  of  the  dyadic 
scattering  amplitude  f  which  was  determined  earlier  in  Section  3.1. 
The  specific  attenuation  for  the  mean  field  propagating  through  an 
unbounded  forest  of  tree  trunks  can  be  found  from  the  imaginary 
part  of  Equation  (4-2-30). 

The  pro*.  _  cion  constants  for  the  horizontally-  and  vertically- 
polarized  components  of  the  mean  field  propagating  through  an  un¬ 
bounded  forest  of  parallel  tree  trunks  can  be  obtained  by  substi¬ 
tuting  Equation  (3-1-1)  into  Equation  (4-2-30).  Recognizing  that 
for  a  forward-scattered  radiowave  propagating  parallel  to  the 
forest  floor. 


8.  =  it/2 


♦  i  "  ♦b  * 


(4-2-1-1) 


then 


C(PP)  =  C(PP) 
n  — n 


(4-2-1-2) 


and  the  equation  for  the  propagation  constants  simplifies  to 


where , 


Kp  =  k0 


(PP)  + 

A  ” 


l  sT)] 


-  Mn/Pn  *  (P-h) 


(4-2-1-3) 


(4-2-1-4) 


(  "  Vn/Nn  '  (P*V) 

and  all  other  parameters  are  defined  in  Section  3.1. 


(4-2-1-5) 


The  specific  a-:  tenuation  attributable  to  tree  trunks  as  ex¬ 
pressed  by  the  imaginary  part  of  Equation  (4-2-1-3)  is  exemplified 
in  Figure  4-6  for  the  case  of  a  radiowave  propagating  parallel  to 
the  forest  floor  (8^*90°).  These  calculations  are  based  uuon  a 
tree  trunk  radius  of  10  centimeters  and  a  trunk  number  density  of 
1  trunk  per  square  meter.  Three  models  are  employed  to  describe 
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Figure  4-6 s  Specific  Attenuation  -  Trunks 


the  permittivity  of  the  wood  [refer  to  Appendix  A] ,  but  only  model 
III  can  be  considered  realistic. 

It  is  apparent  from  Figure  4-6  that,  in  agreement  with  exper¬ 
iment  [33,  36,  39,  40],  vertically  polarized  radiowaves  are  atten¬ 
uated  more  severely  than  are  horizontally  polarized  radiowaves. 

It  is  further  apparent  from  Figure  4-6  that  for  horizontal  polar¬ 
ization  the  specific  attenuation  is  relatively  insensitive  to  the 
particular  choice  of  the  wood  permittivity  model.  This  is  also 
true  for  vertical  polarization  above  70  Megahertz.  This  behavior, 
as  well  as  the  oscillatory  behavior  evident  at  still  higher  fre¬ 
quencies  and  shown  in  greater  detail  in  Figure  4-7  for  several 
values  of  0O,  can  be  attributed  to  resonance  effects.  It  is  also 
apparent  that  above  200  Megahertz  resonance  plays  a  major  role  in 
tree  trunk  scatter  models  and  essentially  precludes  the  utilization 
of  simple  dipole  models  at  UHF  for  all  but  the  smallest  trees. 

In  the  VHF  band,  when  the  radius  of  the  tree  trunk  is  small 
relative  to  the  radio  wavelength  within  the  trunk.  Equation  (4-2- 
1-3)  reduces  to  the  dipole  approximation 

Kh  =  k«^1  +  P^a2X£/(2+XJl)]  (4-2-1-6) 

for  horizontal  polarization,  and  to 

kv  *  k„[l  +  pTta2X£/2]  (4-2-1-7) 

for  vertical  polarization.  These  equations  have  been  compared 
with  those  derived  by  Brown  and  Curry  [25].  Unfortunately,  a 
direct  comparison  was  not  possible  because  Brown  and  Curry  averaged 
their  results  over  a  uniform  distribution  of  tree-trunk  inclination 
angles  and  did  not  provide  intermediate  results.  However,  when 
their  technique  was  used  to  find  the  propagation  constant  appro¬ 
priate  to  an  array  of  parallel  cylinders,  the  derived  expressions 
agreed  exactly  with  Equations  (4-2-1-6)  and  (4-2-1-7) .  This  com¬ 
parison  provided  an  independent  check  of  the  CyberCom  results. 
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4.2. 2.2  Specific  Attenuation  of  Branches 


Branches  are  modelled  as  lossy  dielectric  circular  cylinders 
of  finite  length.  The  salient  scattering  properties  of  individual 
cylinders  can  be  characterized  in  terms  of  the  dyadic  scattering 
amplitude  f  which  was  determined  earlier  in  Section  3.2.  The  spe¬ 
cific  attenuation  of  the  mean  field  propagating  through  an  unbounded 
forest  of  branches  can  be  found  from  the  imaginary  part  of  Equation 
(4-2-25) . 

The  propagation  constants  of  the  mean  scattered  field  propa¬ 
gating  through  an  unbounded  forest  of  branches  can  be  obtained  by 
substituting  Equations  (4-1-2-7)  and  (4-1-2-8)  into  Equation  (*-2- 
25)  and  finding,  for  the  horizontally-polarized  component. 


<h  =  k0|l  +  (p/2)  (-rra24)X4 
and,  for  the  vertically-polarized  component, 

<v  -  k0|l  +  ( p/2 )  ( Tra  2  & )  X & (cos26i  sin20fa 

+  2sin20.  cos 2 6.  )  +  - 

x  b  2+X^ 


572+ X, 


sin 


6b  + 


2 

[+X 


;]| 


(4-2-2-1) 


(4-2-2-2) 


The  specific  attenuation  attributable  to  blanches  as  ex¬ 
pressed  by  the  imaginary  part  of  Equations  (4-2-2-1)  and  (4-2-2-2) 
is  exemplified  in  Figure  4-8  for  the  case  of  a  radiowave  propaga¬ 
ting  parallel  to  the  forest  floor  (0^  *90°).  These  calculations 
are  based  upon  a  branch  radius  of  1  centimeter,  a  branch  length 
of  1  meter,  and  a  branch  number  density  of  1  branch  per  cubic 
meter.  All  branches  are  assumed  to  be  inclined  45  degrees  with 
respect  to  the  vertical,  but  uniformly  distributed  in  azimuth. 
Three  models  are  employed  to  describe  the  permittivity  of  wood 
[refer  to  Appendix  A] ,  but  only  model  III  can  be  considered  real¬ 
istic.  Assumptions  implicit  in  the  development  of  the  thin  branch 
scattering  model  preclude  the  consideration  of  branches  of  signi¬ 
ficantly  greater  radius  over  the  frequency  band  200  -  2000  MHz. 
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In  the  thin-branch  frequency  regime,  the  frequency  behavior 
of  the  specific  attenuation  is  determined  by  the  frequency  be¬ 
havior  of  the  susceptibility  of  a  single  scatterer.  For  branches, 
this  is  clearly  evident  from  Equation  (4-2-2-1)  and  (4-2-2-2) , 
where  the  only  frequency-dependent  parameter  is  the  susceptibility 
of  the  wood  X^.  Further,  since  all  three  models  for  wood  permit¬ 
tivity  satisfy  the  inequality  [refer  to  Appendix  A] . 

XJ  <<  X£  *  39  (4-2-2-3) 

it  is  apparent  from  Equations  (4-2-33) ,  (4-2-2-1)  and  ( 4—2— 2— 2 ) 
that  for  radiowaves  propagating  parallel  to  the  forest  floor 
(0±«  90°) 

Ob  z  (8.686)Jc0  (p/4)  (7Ta2£)X£sin26b  (dB/m)  (4-2-2-4J 

and 

otv  z  (8.686)k0(p/4)  (tra2£)X£cos2eb  (dB/m)  (4-2-2-5) 

so  long  as  @b  is  not  too  near  90  degrees. 

The  empirically-derived  behavior  predicted  by  Saxton  and  Lane 
[refer  to  Figure  1-1]  is  also  shown  in  Figure  4-8.  Although  the 
model  predictions  agree  reasonably  well  with  experiment  so  far  as 
order  of  magnitude  is  concerned,  only  computations  based  upon  wood 
permittivity  model  I  reflect  a  similar  frequency  dependence.  How¬ 
ever,  because  wood  permittivity  model  I  does  not  properly  account 
for  ohmic  losses  within  the  wood,  it  must  be  considered  nonrealis- 
tic  and  the  near  agreement  between  the  branch  model  employing  it 
and  experiment  only  coincidental.  This  is  not  to  say,  however, 
that  the  thin-branch  model  is  invalid.  It  must  be  recognized  that 
the  empirically-derived  behavior  predicted  by  Saxton  and  Lane  re¬ 
fers  to  attenuation  through  a  forest  (tree  trunks,  branches  and 
leaves)  and  not  solely  to  attenuation  through  branches.  Thus,  for 
example,  if  tree  trunks  were  the  dominant  scatterers,  the  sub¬ 
resonant  behavior  contributed  by  the  trunks  [refer  to  Figure  4-6] 
could  easily  mask  the  essentially  frequency-independent  behavior 
predicted  here  for  the  branches  using  wood  permittivity  models  II 
and  III.  Further  study  is  required. 
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Figure  4-8  s  Specific  Attenuation  -  Branches 


4.2.2. 3  Specific  Attenuation  of  Leaves 

Leaves  are  modelled  as  lossy-dielectric  circular  discs. 

The  salient  scattering  properties  of  individual  discs  can  be  char¬ 
acterized  in  terms  of  the  dyadic  scattering  amplitude  f  which  was 
determined  earlier  in  Section  3.3.  The  specific  attenuation  for 
the  mean  field  propagating  through  an  unbounded  forest  of  leaves 
can  be  found  from  the  imaginary  part  of  Equation  (4-2-25) . 


The  propagation  const  mts  of  the  mean  scattered  field  propa¬ 
gating  through  an  unbounded  forest  of  electrically  thin  leaves  can 
be  obtained  by  substituting  Equations  (4-1-3-9)  and  (4-1-3-10)  into 
Equation  (4-2-25)  and  finding,  for  the  horizontally-polarized  com¬ 
ponent 

Kh  -  k.jl  +  (P/2)  (*a*t)Xt[l  -  iJJ  (4-2-3-1) 

and,  for  the  vertically-polarized  component 


kv  «  k0  1+  (p/2)  (Tra2t)XA|l 
where , 


I 


(6)  sin20 


2TI Ajt  (Iicos!ei  +  2I2sin'6 


A)  (4-2-3-2J 


f" 

I2  *  I  d0pA (0)cos26  (4-2-3-3) 

0 


and  p^(0)  is  the  probability  density  function  of  the  leaf  inclina¬ 
tion  angle. 

In  the  frequency  band  200  -  2000  Megahertz,  leaves  may  be  con¬ 
sidered,  nearly  always,  to  be  electrically  thin  [refer  to  Section 
4.1.3] . 

The  specific  attenuation  attributable  to  leaves  as  expressed 
by  the  imaginary  part  of  Equations  (4-2-3-1)  and  (4-2-3-2)  is 
exemplified  in  Figure  4-9  for  the  case  of  a  radiowave  propagating 
parallel  to  the  forest  floor  (0^«9O#).  These  calculations  are 
based  upon  a  leaf  radius  of  5  centimeters,  a  leaf  thickness  of  1 
millimeter,  and  a  leaf  number  density  of  200  leaves  per  cubic 
meter.  The  random  orientation  of  the  leaves  is  described  by  a  pro¬ 
bability  density  function  assumed  to  be  uniform  in  azimuth,  and 
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uniform  in  elevation  over  the  range  of  0  -  30  degrees  (the  leaves 

tr 

tend  to  be  horizontal) .  Three  models  are  employed  to  describe  the  j 

permittivity  of  the  leaves  (refer  to  Appendix  A] ,  but  only  model 
III  can  be  considered  realistic. 


For  electrically  thin  leaves,  the  frequency  behavior  of  the 
specific  attenuation  is  determined  solely  by  the  frequency  be¬ 
havior  of  a  single  scatterer.  This  is  clearly  evident  from  Equa¬ 
tions  (4-2-3-1)  and  ( 4— 2— 3— 2 )  where  the  only  frequency  dependent 
parameter  is  the  susceptibility  of  the  leaf  X£.  Further,  since  all 
three  permittivity  models  satisfy  the  inequality  [refer  to  Appendix 
AJ 

Xi  <<  Xi  "  39  ( 4— 2 —3— 4 ) 

it  is  apparent  from  Equations  (4-2-33) ,  (4-2-3-1)  and  (4-2-3-2) 
that  for  leaves 


* 

•8 


iS 

•18 


ah  =  (8.686)k0(p/2)  (ira2t)Xj[(l  -  H^) 


(4-2-3-5) 


av  s  (8 . 686)k0  (p/2)  (wa2t) X£ (1  -  I2)  (4-2-3-6J 

Especially  noteworthy  is  the  strong  similarity  between  the  fre¬ 
quency  behavior  of  the  leaves  (Figure  4-9]  and  that  of  the  bran¬ 
ches  [Figure  4-10].  This  similarity  is  a  consequence  of  using 
quasi-static  boundary  conditions  for  both  the  leaf  and  the  branch 
in  developing  expressions  for  their  dyadic  scattering  amplitudes 
(refer  to  Section  3] . 


•! 
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The  empirically-derived  behavior  predicted  by  Saxton  and  Lane 
[refer  to  Figure  1-1]  is  also  shown  in  Figure  4-9.  Although  the 
model  predictions  agree  reasonably  well  with  experiment  so  far  as 
order  of  magnitude  is  concerned,  only  computations  based  upon  leaf 
permittivity  model  I  reflect  a  similar  frequency  dependence.  How¬ 
ever,  because  leaf  permittivity  model  I  does  not  properly  account 
for  ohmic  losses  within  the  leaves,  it  must  be  considered  nonreal- 
istic  and  the  near  agreement  between  the  electrically  thin  leaf 
model  employing  it  and  experiment  only  coincidental.  This  is  not 
to  say,  however,  that  the  chin  leaf  model  is  invalid.  It  must  be 
recognize  that  the  empirically-derived  behavior  refers  to  atten- 
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uation  through  a  forest  (tree  trunks,  branches  and  leaves)  and  not 
solely  to  attenuation  through  leaves.  Thus,  for  example,  if  tree 
trunks  were  the  dominant  scatterers,  the  sub-resonant  behavior 
contributed  by  the  trunks  (refer  to  Figure  4-6)  could  easily  mask 
the  essentially  frequency-independent  behavior  predicted  here  for 
the  leaves  by  using  leaf  permittivity  models  II  and  III.  A  similar 
caveat  was  discussed  in  the  case  of  branches. 

♦•3  Salient  Forest  Descriptive  Parameters 

Experimental  verification  of  the  electromagnetic  forest  model 
requires  the  identification  of  measurable  quantitative  parameters 
to  describe  the  forest.  Some  of  these  parameters  are  microscopic 
(e.g.,  the  size  and  relative  permittivity  of  individual  tree  trunks, 
branches  and  leaves);  other  parameters  are  macroscopic  (e.g.,  the 
number  of  trees  per  acre  and  the  number  of  leaves  per  unit  volume 
of  forest).  Several  of  these  parameters  (microscopic  and  macro¬ 
scopic)  have  been  identified  earlier  in  this  report  and  are  dis¬ 
cussed  below. 

For  tree  trunks,  the  specific  attenuation  (a)  experienced  in 
the  low-frequency  (Rayleigh)  regime  by  a  horizontally  or  verti¬ 
cally  polarized  radiowave  propagating  through  an  unbounded  forest 
can  be  determined  from  Equations  (4-2-1-4)  and  (4-2-1-5) .  Using 
these  equations,  CyberCom  has  been  able  to  derive  the  following 
explicit  relations  between  the  specific  attenuation  and  select 
salient  parameters  of  the  trunks 


Ojj  -  (1.13)aapfXJ/|XJl| 

(dB/m) 

(4-3-1) 

av  «  (0 . 286)a2pf X£ 

(dB/m) 

(4-3-2) 

where  a  is  the  trunk  radius  (meters),  p  is  the  trunk  number  den¬ 
sity  (trees  per  square  meter),  f  is  the  frequency  (Megahertz,  and 
X^  is  the  susceptibility  of  the  wood.  Note  that  the  specific 
attenuation  is  directly  proportional  to  the  area  occupied  by  the 
trunks  and  to  the  frequency.  Analogous  expressions  for  specific 
attenuation  have  not  yet  been  determined  (outside  the  Rayleigh 
regime, at  higher  frequencies). 
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For  branches.  Cyber  Coin  has  been  able  to  derive  from  Equations 
(4-2-2-4)  and  (4-2-2-S)  the  following  explicit  relations  between 
the  specific  attenuation  due  to  electrically-thin  branches  and 
select  salient  parameters: 


“h  - 

(0. 143)a2 Apfxj  sinaeb 

(dB/m) 

(4-3-3) 

av  - 

(0. 286) a2£pfX£  cosa6b 

(dB/m) 

(4-3-4) 

where  a  is  the  branch  radius  (meters),  l  is  the  branch  length 
(meters) ,  p  is  the  branch  number  density  (branches  per  cubic  meter) , 
f  is  the  frequency  (Megahertz),  X^  is  the  susceptibility  of  the 
wood,  and  0b  is  the  angle  that  the  branch  makes  with  the  vertical. 

It  is  apparent  from  these  equations  that  the  specific  attenuation 
depends  upon  the  square  of  the  branch  radius,  but  only  linearly 
upon  branch  length,  number  density,  and  X£. 

For  leaves,  CyberCom  has  been  able  to  derive  from  Equations 
(4 -2-3-5)  and  (4-2-3-6)  the  following  explicit  relations  between 
the  specific  attenuation  due  to  electrically  thin  leaves  and 
select  salient  foliage  parameters: 

«  (0.286)a2tpfx£(l  -  Ix)  (dB/m)  (4-3-5) 

av  «  (0. 286) a2tpf X£ (1  -  I2)  (dB/m)  (4-3-6) 

where  a  is  the  leaf  radius  (meters),  t  is  the  leaf  thickness 
(meters),  p  is  the  leaf  number  density  (leaves  per  cubic  meter), 
f  is  the  frequency  (Megahertz) ,  and  X^  is  the  susceptibility  of 
the  leaves.  The  parameters  1^  and  I2  are  related  to  the  proba¬ 
bility  density  function  of  the  leaf  inclination  angles  and  are 
defined  in  Equation  (4-2-3-3) .  It  is  apparent  from  these  equations 
that  the  specific  attenuation  depends  upon  the  square  of  the  leaf 
radius,  but  only  linearly  upon  the  leaf  thickness,  number  density, 
and  X£. 


It  may  be  noted  that  each  of  the  above  expressions  is  of 
the  form 

a  =  kVpfG  (4-3-7) 

where 

k  a  constant 
V  *  volume  of  element 
p  *  density  (no./ms) 
f  =  frequency  in  MHz 
G  =  geometrical  factor 
Vp  *  fractional  volume 

As  it  is  required  in  several  developments  in  this  report  that 
fractional  volume  be  small,  a  preliminary  evaluation  of  this  quan¬ 
tity  has  been  made  using  Reference  15.  This  lists  for  forests  in 
many  parts  of  the  world  the  following: 

A  =  basal  area  in  m2 /hectare  (10*m2) 

*w  *  dry  mass  of  stem  wood 
"b  *  dry  mass  of  branches 
Mf  *  dry  mass  of  foliage 

A  cursory  average  of  values  for  the  United  States  gives  the 
following  fractional  volumes: 

Trunks  A  x  10”1*  *  .0064 

Branches  .0064  x  A*b/Mw  *  .0013 

Leaves  .0064  x  M^/Mw  =  .00016 

These  values  are  clearly  very  small. 
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5.0  Anisotropic  Forest  Slab  Model 

Results  were  obtained  above  for  an  unbounded  continuous  medium 
equivalent  to  an  infinitely  high  forest  of  model  trunks ,  branches  or 
leaves.-  In  this  section  the  medium  is  bounded  at  a  forest-air  up¬ 
per  interface  to  provide  a  more  realistic  overall  model.  This 
necessitates  starting  again  from  Maxwell's  equations,  inserting 
the  dyadic  permittivities  from  Section  3.  The  resulting  inhomo¬ 
geneous  wave  equation  for  an  anisotropic  medium  is  converted  to  a 
homogeneous  equation  for  an  isotropic  medium  by  the  substitution 

z*  «  /e“/e“  z  (5-0-1) 

The  vector  potentials  are  found  by  an  asymptotic  evaluation 
using  integrations  in  the  complex  plane.  The  reflected  component 
is  found  to  include  a  lateral  (tree-top)  component.  These,  plus 
the  direct  wave  are  the  potentials  used  to  find  the  corresponding 
vertical  components  of  the  E  field  received  by  a  vertical  dipole. 

5.1  Model  Formulation 

The  basic  slab  model  for  the  forest  is  shown  in  Figure  5-1. 

The  trees  and  vegetation  are  assumed  to  be  distributed  uniformly 
between  a  smooth  forest  floor  and  the  air  interface  at  height  h. 
Earlier  efforts  [93,  104,  107]  have  shown  that  the  presence  of  the 
ground  complicates  the  model  significantly.  The  complications  can 
be  avoided,  however,  by  allowing  the  ground  plane  to  recede  to 
z  -*•  — ®  so  that  the  model  reduces  to  the  half-space  representation 
shown  in  Figure  5-2.  This  simplification  provides  a  very  good 
approximation  to  the  basic  slab  model  if  neither  the  transmitting 
or  receiving  antenna  is  located  too  close  to  the  ground  [107]. 

Consider  a  transmitting  antenna  which  is  representable  as  a 
vertical  electric  (Hertzian)  dipole  having  a  time-harmonic  current 
moment  Idl  •  exp{jwt)  and  immersed  a  distance  d  below  the  inter¬ 
face  of  two  semi-infinite  media  [refer  to  Figure  5-2].  The  iso¬ 
tropic  upper  half-space  (z  >  0)  represents  the  air  and  is  charac¬ 
terized  by  the  permittivity  e0,  and  permeability  y0  of  free- 
space.  The  electrically  anisotropic  lower  half-space  (z  <  0) 


represents  the  forest  and  is  characterized  (at  least  so  far  as  the 
mean  fields  are  concerned)  by  the  effective  dyadic  permittivity 
£e0  And  the  free-space  permeability  p0.  The  relation  between  the 
relative  effective  dyadic  permittivity  £  and  the  biophysical  para¬ 
meters  of  the  forest  has  been  described  in  Section  4  of  this  report. 


5.2  Mean  Fields 


In  any  charge-free  medium  where  the  electromagnetic  fields 
vary  harmonically  as  exp{jwt}.  Maxwell's  equations  may  be  written 
in  the  form 


V  x  E  =  -jwB 


V  •  D  =  0 


V  x  H  *  jdiD  +  J 


7  •  B  =  0 


(5-2-1) 


If  the  medium  is  electrically  anisotropic  and  characterized  by  the 
relative  permittivity  dyadic  £,  the  constitutive  relations 

D  =  e0e  •  E  B  =  y0H  (5-2-2) 


may  be  introduced  and  Maxwell's  equations  re-written  as 

V  x  E  =  -j(i)U0H  V  •  £  •  E  *  0 

V  x  H  =  ja)e0,e  »E+JC  7  *  H  =  0 


(5-2-3) 


Because  the  forest  may  be  considered  uniaxially-anisotropic  with 
respect  to  the  mean  fields  (refer  to  Section  4] , 

£  =  etx°x°  4  et^°  +  £2£°£0  (5-2-4) 

can  be  used  to  represent  the  effective  relative  complex  permit¬ 
tivity  of  the  forest,  and  E,  H  the  mean  fields. 

Because  for  any  vector  A, 

7  •  (V  x  A)  =  0  ,  (5-2-5) 

the  relation 

7  •  H  -  0  (5-2-6) 
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Figure  5-2 :  Half-Space  Model  for  the  Forest 


(5-2-7) 


# 

suggests  that  there  is  a  vector  A  satisfying  the  relation 

H  ■  j0)£oEt  (V  x  A) 

If  this  expression  is  introduced  into  Maxwell's  equations,  it 
may  be  shown  that  A  satisfies  the  equation 

72A  +  k2E  •  A  +  e“l(e-etI)  *7(7*  A)  «  (-jweeet) " 1 3^  (5-2-8) 

and  that  the  mean  electric  and  magnetic  fields  may  be  derived  from 
A  using  the  .relations 

E  *  k|etA  +  7(7  •  A)  (5-2-9) 

H  -  jwe0et(7  x A)  (5-2-10) 

For  an  electrically  isotropic  medium  where 

£t  ■  e2  *  £  (5-2-11) 

Equation  (5-2-8)  reduces  to  the  inhomogeneous  vector  wave 
equation 

72A  +  k|cA  =  (-ju>E.e)“lA  (5-2-12) 

and  A  is  known  as  the  electric  Hertz  potential  [111] . 

In  Region  I  (above  the  forest)  where  there  are  no  electro¬ 
magnetic  sources 

-  0  (5-2-13) 

the  vector  potential  (A1)  satisfies  the  homogeneous  vector  wave 
equation 

72AJ  +  k2AJ  -  0  (5-2-14) 

In  Region  II  (within  the  forest)  where  a  vertical  electric 
(Hertzian)  dipole  of  current  moment4  Idl  is  immersed  a  distance 
d  below  the  top  of  the  forest  canopy,  the  vector  potential  (A11) 
satisfies  the  inhomogeneous  vector  equation 


*  The  scalar  coefficient  jwe#£t  has  been  included  for  sathematical  convenience. 
+  The  harmonic  factor  exp{ jut}  has  been  assumed  and  suppressed. 
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7 2 A11  +  k^e-A11  +e“l  (£  -  etI)  ‘  7(7  •  A11) 


=  (-ja)e0et)“lIdld(x)6(y)6(2 +d)2? 


(5-2-15) 


Because  the  current  moment  Idl  is  directed  along  the  z-axis, 
a  solution  A11  will  be  sought  having  only  a  z-directed  component, 
i.e. 


(0,0,a“) 


(5-2-16) 


Equation  (5-2-16)  permits  Equation  (5-2-15)  to  be  reduced  to  the 
inhomogeneous  scalar  wave  equation 


32A^^  32A^ 

- L_  +  - L_ 

3x2  3y2 


—  +  (“) 
3y2  Vt/ 


32AJI 

- —  +  k2e  A*1 

3z2  22 


(-ju>c0et)”  Idl6  (x)  6  (y)  6  (z+d) 


By  introducing  the  anisotropy  factor 


(5-2-17) 


a  -  Ut/ez) 


and  the  change  of  variable 


z  •  *  az 


(5-2-18) 


(5-2-19) 


Equation  (5-2-17)  can  be  re-written  as 


V^A*1  +  kje^j1  -  (-jwe.et)’1Idl‘6(x)6(y)6(z'  +d')  (5-2-20) 


where  7 '  represents  the  del  operator  transformed  according  t.o 
Equation  (5-2-19) „  Note  that  as  a  consequence  of  the  change  of 
variable,  the  inhomogeneous  scalar  wave  equation  for  the  aniso¬ 
tropic  medium  (Equation  (5-2-17) J  has  been  transformed  into  the 
homogeneous  scalar  wave  equation  for  an  equivalent  isotropic  me¬ 
dium  (111] .  The  possibility  of  such  a  transformation  for  unbounded 
anisotropic  media  apparently  was  first  suggested  by  Clemmow  (121] . 

The  solution  to  Equation  (5-2-20)  consists  of  two  parts:  the 
complementary  solution  (A*^  )  to  t  homogeneous  scalar  wave 
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equation,  and  the  particular  solution  (A**)  to  the  inhomogeneous 
scalar  wave  equation.  The  particular  solution  can  be  obtained  by 
first  introducing  the  three-dimensional  Fourier  transform  pair 

Vi>  =  jAz(r)exp{jg  •  r}dr  (5-2-21) 

Az(r)  *  (27t)”SA2(B)exp{-je  •  r}d6  (5-2-22) 

and  subsequently  using  Jordan's  Lemma  to  reduce  the  particular 
solution  to  the  two-dimensional  integral 


.II  /  * 
AZ.P<£> 


A“P<£'Z) 


-Idl  f  exp{-j(Bt*6  +  i  Jz+d|)} 

—5 - J - - - i -  dBt  (5-2-23) 


8tt  uie0e 


where 


T2  =  a(k^z  -  .  T.  _  j  I  T 2  |  (5-2-24) 

The  condition  that  t*  <  0  is  necessary  to  ensure  the  convergence 
of  the  integral  as  |z+d| 

Recognizing  that  the  particular  solution  [Equation  (5-2-23)] 
is  a  two-dimensional  integral,  complementary  solutions  to  Equations 
(5-2-14)  and  (5-2-20)  will  be  sought  utilizing  the  two-dimensional 
Fourier  transform  pair 

A2  (Bt,z)  *  jAz  (£#z)exp{  jJ3t  •  £>d£  (5-2-25) 

Az(£,z)  «  (2tt)”2|az  (£t,v3:exp{-jBt  .  p}dBt  (5-2-26) 


with  appropriate  consideration  given  to  Mie  transformation  of 
Equation  (5-2-19).  Substitution  of  Equation  (5-2-26)  into  Equations 
(5-2-14)  and  (5-2-20)  yields  the  complementary  solutions 

i5;«Ut..)  *  C+'II(It)exp{  jTi  ,2z}  +cf,II(3t)exp{-jt1  2z}  (5-2  *27) 

where , 

T,  “  ^  “  3  I T  J  i  (5-2-28) 
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T2  «  a(k*e2  **  B*)**  *  t2  -  j  |  |  (5-2-29) 

To  ensure  that  the  electromagnetic  fields  remain  finite  as  |z| 
it  is  necessary  that 

cJ(Bt>  =  0  (5-2-30) 

C^I(Bt)  =  0  (5-2-31) 

II 

Following  a  normalization  of  the  coefficients  C*  (B*.)  and 

I  + 

C_(Bt)  and  the  addition  of  the  complementary  and  particular  so¬ 
lutions  associated  with  the  wave  equation  for  Region  11,  the  po¬ 
tential  functions  for  Region  1  and  11  may  be  written,  respectively, 
in  the  following  forms: 

A*(£,z)  *  -  [fKB^expC-jTjZjexpt-jB  •  £}dBt  (5-2-32) 

Bit  <*>e0e2  * 

A2T(£,z)  *=  — —  f  l  (l/t2)exp{-jT  |z+dj ) +N(Bt)exp{  jt2z}l 

we0e„  J 

z 

•  exp{-jBt  •  £>d£t  (5-2-33) 

Because  the  tangential  components  of  the  electromagnetic 
fields  must  be  continuous  across  the  air-forest  interface. 

Equations  (5-2-9) ,  (5-2-10)  and  (5-2-16)  require 

aJ(£,0)  -  etA”(£,0)  (5-2-34) 

3lA^(£,0))/3z  -  3(A2I(£,0)]/3z  (5-2-35) 

Substitution  of  Equations  (5-2-32)  and  (5-2-33)  into  Equations 
(5-2-34)  and  (5-2-35)  reveals  that 

M(Bt)  *  I2et/(T2  + etTj) )exp{-jT2d}  (5-2-36) 

N(Bt)  *  t (t2  -  etx1)/(T2  +  etT1)T2]exp{-jx2d)  (5-2-37) 
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so  that  the  potential  functions  for  Region  I  and  Region  II  are, 
respectively, 

A*(£,z)  «  — —  f - - -  exp{-jx  d}exp{-jx  z) 

*  87r2a)e0e2  J(T2+etTl) 


*  exp{- jBt  •  £}dgt 


(5-2-38) 


-Idl  f 

r  -jx2 |z+d| 

•  -  +1 

Vi) 

jx2  (z-d)-. 

2  1 
8n  u)e0e„/ 
z 

l  ^  1 

Vt2  +  etTj> 

t,  J 

•  exp{-j$t  •  £}dBt 


(5-2-39) 


Equations  (5-2-38)  and  (5-2-39),  in  conjunction  with  Equations 
(5-2-9)  and  (5-2-10),  constitute  the  formal  solution  for  the 
mean  electromagnetic  fields  of  a  vertical  electric  (Hertzian) 
dipole  immersed  in  an  anisotropic  forest  half-space. 

5.3  Asymptotic  Evaluation 

The  integral  representations  for  the  potential  functions 
afforded  by  Equations  (5-2-38)  and  (5-2-39)  are  not  amenable  to 
exact  analytic  evaluation.  Fortunately,  analytic  asymptotic 
approximations  can  be  derived  which  will  prove  adequate  for  most 
engineering  applications.  However,  before  proceeding  with  the 
development  of  these  asymptotic  approximations,  it  will  prove 
expedient  to  introduce  the  transformations 


Bx  *  X  cos  \p 


8y  “  X  sin 


Px  ■  P  cos  <t> 


Py  “  P  sin  $ 


(5-3-1) 


so  that  with  the  help  of  the  identity 

1 

Jo(Xp)  *  j--  I  exp{-jXp  cos  (i|*-$)}dij/ 


(5-3-2) 
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the  potential  function  within  the  forest  can  be  written  as 


*”<£.*>  -  “Idl 


4itwe 


°ez  0  L 


-jT4I*-d| 


+  r(t 


jtjtz-dn 

J 


J^(Xp) 


XdX 


(5-3-3) 


where 


r{Ti»V 


T*  “  etT» 


(5-3-4) 


Ta  *  etTi 


may  be  recognised  as  the  Fresnel  reflection  coefficient  associated 
with  the  air-forest  interface  and,  as  before. 


-  i I t~ |  -  (k*  -  xM1* 


(5-3-5) 


*r2  *  Tj  -  j  | t; |  -  a(k2ez  -  X1)** 


(5-3-6) 


If  the  identity  [111,  Equation  (5.35a)] 

e-jic|s+d| 


exp{-jko/eR} 


J.  (Xp)XdX 


(5-3-7) 


where 


ic  °  <k2e  -  X1)' 


R  -  IP*  +  (s-d)2]** 


(5-3-8) 

(5-3-9) 


is  analytically  continued  into  the  complex  z -plane  using  the 
transformation  of  Equation  (5-2-19) ,  the  potential  function  in 
the  forest  can  be  more  conveniently  written  in  the  form 


♦  A'*’ 


(5-3-10) 


where 


(5-3-11) 


.«s>_  -lidl  . 

4trwe#c8  aR^ 

Rd  -  [p1  +  a1  (z+d) 1 J^*  (5-3-12) 

is  the  component  of  the  potential  function  associated  with  the 
direct  wave,  and 

jTa(z-d) 

Air>  *  -2-1-—  I  r(T1,TJ)  5 - •  J0(Xp)XdX  (5-3-13) 

4irweae.J5  Tj 


is  the  component  of  the  potential  function  associated  with  the 
reflected  wave. 

The  asymptotic  evaluation  of  is  most  readily  achieved  by 

first  analytically  continuing  the  real  integration  variable  X  into 
the  complex  plane.  Unfortunately,  the  integral  presented  by  Equa¬ 
tion  (5-3-13)  presents  a  mathematical  problem:  the  path  of  inte¬ 
gration  begins  at  the  origin.  This  can  be  eliminated  by  expressing 
the  Bessel  function  J0(Xp)  in  terras  of  the  Hankel  functions  of  the 
first  and  second  kinds  (h11}(Xp)  and  H*2)(Xp)I,  and  then  using 
analytic  continuation  to  express  the  Hankel  function  of  the  first 
kind  in  terms  of  one  of  the  second  kind  1111)  so  that  the  electro¬ 
magnetic  waves  exhibit  the  proper  asymptotic  behavior  at  infinity 
as  might  be  expected  with  the  harmonic  time  dependence  exp{jwt}. 

As  a  consequence  of  these  mathematical  manipulations,  the  poten¬ 
tial  associated  with  the  reflected  wave  can  be  recast  into  the 
form 


It  is  especially  important  to  note  that  without  the  con¬ 
straints  imposed  by  Equations  (5-3-5)  and  (5-3-6)  on  the  imaginary 
parts  of  tj  and  t2,  the  integrand  of  Equation  (5-3-14)  would  not 
be  uniquely  determined.  Generally  speaking,  the  integrand  is  four¬ 
valued  (corresponding  to  the  four  possible  combinations  of  signs 
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in  th«  square  roots  defining  Tt  and  t#),  and  its  Riemann  surface 
consists  of  four  sheets.  To  insure  the  convergence  of  the  integral 
at  infinity  and  the  uniqueness  of  the  integral  everywhere,  the  in¬ 
tegration  path  in  the  complex  X-plane  must  be  constrained  to  the 
permissible  (upper)  sheet  defined  by  Equations  (S-3-S)  and  (5-3-6). 
This  can  be  achieved  by  joining  the  branch  points  (i.e.,  those 
points  where  x,  and  t2  are  sero)  by  two  (essentially  arbitrary) 
branch  cuts  and  insuring  that  the  branch  cuts  are  not  intersected 
by  the  path  of  integration. 

An  asymptotic  approximation  for  the  potential  function  in  the 
forest  suitable  at  relatively  large  distances  from  the  transmitter 
(k0p  >>  1)  can  be  obtained  by  evaluating  Equation  (5-3-14)  using  the 
method  of  steepest  descent  [111].  This  approach  involves  the  con¬ 
formal  re-mapping  of  the  integrand  from  the  X-plane  to  the  w-plane 
using  the  transformation 

X  «  k./e^sinw  (5-3-15) 

to  eliminate  the  branch  cut  associated  with  t2;  the  subsequent 

deformation  of  the  Integration  contour  to  coincide  with  the  path 

of  steepest  descent  passing  through  the  saddle  point  associated 

with  the  exponential  factor;  and,  finally,  the  approximate  evalua¬ 
te 

tion  of  the  integral  defining  A'  . 

z 

The  complex  v-plane  associated  with  the  transformation  of 
Equation  (5-3-15)  is  shown  in  Figure  5-3.  Under  this  transfor¬ 
mation,  the  coordinate  axes  of  the  X-plane  are  mapped  into  lines 
in  the  X-plane  defined  by* 

tanh  w"  *  (X"/2) tan  w*  ,  (X*-axis)  (5-3-16) 

Z 

and 

tanh  w*  =  ~ (2/X") tan  w»  ,  (X"-axis)  (5-3-17) 

z 

The  branch  points  associated  with  xx  are  located  at 

wB  =  ±(2n +  1)ti/2  +  ,  (n-1,2,...)  (5-3-18) 


*  Predicated  upon  the  condition  that  |x  I  ■  | X*  -  jX"|  «  1. 

Z  Z  Z 1 


If  the  associated  branch  cuts  are  defined  so  that 

Im{ t  j }  «  Im{k„  (1  -  ezsin2w)*}  =  0  (5-3-19) 

then  from  the  branch  point  to  the  origin, the  branch  cut  is 
defined  by* 

tanh  w"  a  (X^/2)tan  w'  ,  (w'w”  > 0)  (5-3-20) 

z 

and  from  the  origin  to  infinity  by* 

tanh  w"  a  -(2/X")tan  w'  ,  (w'w"  <  0)  (5-3-21) 

z 


Note  that  the  branch  cut  proceeds  from  the  branch  point  along  the 
X'-axis  to  the  origin,  experiences  a  clockwise  angular  rotation  of 
it/2  radians,  and  then  proceeds  from  the  origin  along  the  X-axis  to 
infinity. 

Poles,  associated  with  the  denominator  of  the  reflection 
coefficient  r(Tj,T2),  occur  wherever 

cos  w  +  a(l  -  ezsin2w)^*  =  0  (5-3-22) 

These  poles  can  be  found  on  the  top  sheet  where* 

wp  a  +U/4  +  jX£/4)  (5-3-23) 

and  on  the  bottom  sheet  where* 

wp  *  ±(3tr/4  +  jX£/4)  (5-3-24) 

Anticipating  an  asymptotic  evaluation,  the  Hankel  function 
in  Equation  (5-3-14)  can  be  replaced  by  its  approximation  for 
large  arguments 


h!z) (Xp> 


exp{-j  (Xp  -  tt/4 )  ) 


(5-3-25) 


Introducing  the  transformation  defined  by  the  mapping  of 
Equation  (5-3-15)  and  the  geometric  transformation  suggested  by 

*  Predicated  upon  the  condition  that  |  X  I  «*  |  X ’  -  jX"|  <<  l. 

z  z  z 
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Equation  (5-2-19),  vi2. 

p  *  R  sin  6 
r  r 

-a(z-d)  *  Rr  cos  6r 

recasts  the  potential  function  into  the  form 


-Idl  ej  tt/4 

4TTwe0e„ 

z 


(k  /FT)  **  r 

- ? — - - c  r (w)  /slnw 

a  (2irRr  sin  0r)*  * 


(5-3-26) 


•  exp{-jk0/e~  R_  cos  (w  -  0_)  }dw  (5-3-27) 

z  r  r 


When 

k0|/e“Rir|»l  (5-3-28) 

z  r 

Equation  (5-3-27)  can  be  evaluated  asymptotically  by  defining  the 
integration  path  P  (which,  in  the  X-plane,  coincides  with  the 
X'-axis)  into  the  steepest-descent  path  SDP  (which,  in  the  w-plane, 
passes  through  the  saddle  point  w»0r)  defined  by 

Re(/e^  Rr  cos  (w-0r)}  *  Re(/e^  Rr)  (5-3-29) 


Integration  of  A 


(r) 

z 


along  the  steepest-descent  path  yields 


Idl 


4TTwe0e„ 
*  z 


r  (©r) 


exp{-jko/e^  Rr) 


aR. 


(5-3-30) 


where 

cos0_  -  /e7(l  -  e_sin*0_) ** 

r < ©r)  *  - - - - - ? - —  (5-3-31) 

cos0_  +  /e7(l  -  e-Sin2©-)1* 
r  t  z  r 

is  the  Fresnel  reflection  coefficient  associated  with  the  wave 
specularly  reflected  from  the  air-forest  interface. 


The  validity  of  Equation  (5-3-30)  is  predicated  not  only 
upon  the  condition  that 


5- 


9 


(5-3-32) 


k0|/e^  Rrl  >>  1 


but  also  upon  the  condition  that  the  value  of  the  reflection  coef¬ 
ficient  r(w)  appearing  in  Equation  (5-3-27)  does  not  vary  appre¬ 
ciably  along  the  steepest  descent  path  in  the  vicinity  of  the  sad¬ 
dle  point,  i.e. 

r*(w0)  s  0  (5-3-33) 

This  assumption  will  prove  justified  if  [87] 

k*|Rr(er  -  0C)2|  »  1  (5-3-34) 

where 

6c  -  Arc  sin  (1//F^)  (5-3-35) 

For  the  forest,  the  effective  permittivity  e  is  close  to  unity. 
Further,  because  the  distance  p  between  the  transmitter  and  re¬ 
ceiver  will,  nearly  always,  be  much  greater  than  their  height 
differential  |z-d],  the  angle 

t?r  *  Arc  tan  [p/a|z-d|]  (5-3-36) 


will,  nearly  always,  be  close  to  ir/2.  As  a  consequence,  the  re¬ 
quirement  expressed  by  Inequality  (5-3-34)  is  usually  much  more 
stringent  than  that  expressed  by  Inequality  (5-3-32) . 

Because  6r  will,  nearly  always,  be  close  to  tt/2,  the  defor¬ 
mation  of  the  contour  P  into  the  steepest-descent  path  SDP  nearly 
always  results  in  the  capture  of  the  branch  point  w0.  In  order 
to  avoid  crossing  the  branch  cut  when  the  branch  point  is  cap¬ 
tured,  an  additional  line  integral  about  the  branch  cut  must  be 
introduced  [refer  to  Figure  5-3] .  Integration  about  the  branch 
cut  yields  the  asymptotic  result  that  [87] 


U) 


Idl 

2ttwe0 


1 _ 


k0 (e2-l)*/p 


exp(-jkp  [p+a(ez-l)l*|z-d|  ] } 
[  (e.-D^p  -  a|z-d|  ]  #/2 


(5-3-37) 
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subject  to  the  constraint  that 


k.  1  (*x~1)Hq  “  al  *“dl  I 

(k,eta|**d|  1  ** 


»  1 


(5-3-38) 


p  »  | a |  Is-dj/U^-l)*1 


(5-3-39) 


then 


<*>  *  JM3L  . 


x  exp{-jk#(p +a(ex-l)^|s-d|) } 


2Trwe0  k. (e^-1) 


(5-3-40) 


Actually,  the  validity  of  Equation  (5-3-37)  is  predicated,  not 
only  upon  Inequality  (5-3-38) ,  but  also  upon  the  condition  that 
the  reflection  coefficient  f(w)  appearing  in  Equation  (5-3-27) 
does  not  vary  appreciably  along  the  steepest-descent  path  in  the 
vicinity  of  the  branch  point,  i.e. 


r'(vB)  *  0 


This  assumption  is  justified  if 


0co|Rr|  le^-ll)*8  »  1 


(5-3-41) 


(5-3-42) 


Observing  that  when  a*l,  the  phase  factor  appearing  in  the 
exponents  of  Equations  (5-3-37)  and  (5-3-40)  can  be  written  in 
the  form 


k0lp+  (e  -1)  **)*-£  |  ]  -  k0 /£T  ( d/cos  0) 


+  k, Ip  -  d(sin6_/cos6_)  -  I z | (sin0„/cos8,J  ) 


+  k,/e^(  |z |/cos6c) 


(5-3-43) 


where 


0C  «  Arc  sin  (l//e~) 


(5-3-44) 
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The  three  terms  of  Equation  (5-3-43 )  suggest  that  the  vector  po¬ 
tential  can  be  associated  with  a  so-called  "lateral"  wave 

which  propagates  from  the  transmitter  up  through  the  forest  at  the 
critical  angle  0_  to  the  air-forest  interface,  through  the  air 

w 

along  the  air-forest  interface,  and  down  through  the  forest  at  the 
critical  angle  from  the  air-forest  interface  to  the  receiver. 

As  a  consequence  of  this  asymptotic  evaluation,  a  radiowave 
propagating  within  the  forest  may  be  considered  to  consist  of 
three  components  -  the  direct  *fave ,  the  reflected  wave,  and  the 
lateral  wave  -  and  the  vector  potential  within  the  forest  can  be 
written  in  the  form 


(d) 


+  A. 


(r) 


+  A, 


(i) 


(5-3-45) 


where  A^d^  is  the  vector  potential  associated  with  the  direct 
wave  (Equation  (5-3-11)],  A^r^  is  the  vector  potential  associated 
with  the  reflected  wave  [Equation  (5-3-30) ] ,  and  A^*^  is  the 
vector  potential  associated  with  the  lateral  wave  [Equation  (5-3- 
37)]. 


5.4  Electromagnetic  Fields 

The  electromagnetic  fields  (E  and  H)  within  the  forest  can 

be  derived  from  the  vector  potential  A*1  [Equation  (5-3-45)]  by 

z 

using  Equations  (5-2-9)  and  (5-2-10).  Because  the  electromagne¬ 
tic  fields  fire  related  linearly  to  the  vector  potential,  just  as 
the  vector  potential  exhibits  three  asymptotic  components  (the 
direct,  reflected,  and  lateral  waves),  so  too  will  the  electro¬ 
magnetic  fields.  For  example,  the  vertically-polarized  components 
of  the  electric  field  vector  E,  derived  under  the  asymptotic 
approximation  from  Equation  (5-2-9)  upon  substitution  of  Equations 
(5-3-11),  (5-3-30)  and  (5-3-37),  are 
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Idl  sin2  6^ 


exp(-jk#/e^  Rd) 


(5-4-1)  j 
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(5-4-2) 


,(r) 


exp{-jk#/e^  Rr) 


(wy  \  exp{-jk 

— - )  Idl  *in*6  r (8  )  - - 

«»  /  r  r  Rr 


E*l)  »  60  Idl 
z 


a  vexp{-jk#Ip  +  a (cx-l) ** I *-d I } 


(*) 


where, 

Rd  -  [p2  +  a2(s+d)2]* 

Rr  -  IP2  ♦  a2  (*-d)  21I* 

8d  *  Arc  tan 
6_  *  Arc  tan  - 2 - 

r  l  • 

co*8_  -  /Erd-c.iin'e.)^ 

r  (er)  -  - £ - i - 5 - £— 

coaer  +  /ej(l-exsin26r)^ 
a  *  /et/,£z 


(5-4-3) 

(5-4-4) 

(5-4-5) 

(5-4-5) 

(5-4-7) 


(5-4-8 


(5-4-9) 


Note  that  in  thf»  limiting  case  of  an  isotropic  forest,  the  aniso¬ 
tropy  factor  a  +  1,  and  the  above  results  reduce  to  those  obtained 
by  Staiman  and  Tamir  [114,  Equation  (28)  with  a«-ir/2].  In  this 
limiting  isotropic  case,  the  parameters  Rd,  Rf,  6d  and  6r  acquire 
the  geometric  interpretation  shown  in  Figure  5-2. 

The  remaining  component  of  the  electromagnetic  field,  the 
magnetic  field  intensity,  can  be  obtained  similarly.  In  this  re¬ 
port,  however,  attention  will  be  focused  upon  the  vertical  electric 
field  E  ,  the  only  component  received  by  a  vertical  whip  antenna. 


6.0  Forest  Pulse  Response 

The  anisotropic  half -space  model  of  the  forest  developed  in 
Section  5  of  this  report  is  a  time-harmonic  model.  The  trans¬ 
mitting  antenna  was  assumed  to  be  a  vertical  electric  (Hertzian) 
dipole  having  a  time-harmonic  current  moment  Id£*exp{ jut) .  How¬ 
ever,  because  the  effective  dyadic  permittivity  characterizing  the 
equivalent  continuous  medium  of  the  mean  wave  is  linear,  the  model 
can  be  extended  using  linear  system  theory  to  accommodate  the 
transmission  of  arbitrarily-modulated  waveforms.  This  extension 
is  especially  important  if  the  model  is  to  be  employed  in  analysis 
and  evaluation  of  wideband  spread-spectrum  radio  communication 
systems . 

Following  a  brief  review  of  linear  system  theory  in  Section 
6.1,  the  forest  transfer  function  F(f,t)  is  identified  in  Section 
6.2  and  employed  in  Section  6.3  to  define  the  forest  pulse  re¬ 
sponse.  In  Section  6.4,  the  forest  pulse  response  is  evaluated 
for  the  important  practical  case  of  a  broadband,  rectangular, 
r-f  pulse. 

6.1  Linear  System  Theory 

A  system  is  a  collection  of  interrelated  components  or  objects 
for  which  there  is  specified  a  set  of  dynamic  variables  called 
excitations,  or  inputs,  and  another  set  called  responses,  or  out¬ 
puts  [145] .  The  radio  transmission  channel  through  the  forest 
represents  such  a  system.  The  tree  leaves,  branches,  and  trunks 
constitute  tl.e  collection  of  components  or  objects  which  are  inter¬ 
related  by  mutually-induced  currents  excited  within  them  by  electro¬ 
magnetic  fields.  The  dynamic  variable  representing  the  input  of 
this  system  can  be  the  time-harmonic  electric  current  momemt  Id£ 
of  an  electrically  small  (Hertzian)  dipole  transmitting  antenna. 

The  dynamic  variable  representing  the  output  of  this  system  can 
be  the  electric  field  E  at  the  location  of  the  receiving  antenna. 

The  objective  of  systems  analysis  is  to  determine  how  such  a 
collection  of  components  of  objects  behaves  when  subjected  to  an 
arbitrary,  but  specified,  excitation.  For  the  analysis  of  linear 


6-1 


systems  (such  as  the  forest  transmission  channel) ,  which  (by  de¬ 
finition)  satisfy  the  principle  of  superposition  so  that 


F*alxl  +  a2x2*  "  ajFfx^  +  a2P[x2] 


(6-1-1) 


where  F[x]  represents  the  output  response  of  the  system  to  the 
input  excitation  x,  there  are  two  basic  approaches.  The  first 
of  these,  which  might  be  termed  the  "direct"  approach,  consists 
in  directly  solving  the  input-output  relation 

y  -  F [x]  (6-1-2) 

subject  to  a  known  set  of  initial  conditions  or  boundary  condi¬ 
tions.  In  the  "indirect"  approach,  the  input  x  is  first  resolved 
into  a  set  of  elementary  functions  all  of  which  are  similar  in 
form.  The  response  of  the  system  to  each  elementary  component 
(presumably  determined  more  easily  than  the  response  of  the  system 
to  an  arbitrary  input)  is  then  obtained  and  the  responses  to  all 
the  elementary  components  of  the  input  added  to  obtain  (by  virtue 
of  system  additivity)  the  output  corresponding  to  the  input. 


Any  one  of  a  number  of  elementary  components  may  be  used  for 
the  decomposition  of  the  input.  Two  often-employed  choices  for 
the  elementary  components  are  the  Dirac  del;  _a  function  and  the 
complex  exponential.  Resolution  of  the  input  x  into  a  continuum 
of  Dirac  delta  functions  may  be  achieved  with  the  help  of  the 
sifting  integral 


x(t)  “  J  x(t) 6 (t— T)dT 


(6-1-3) 


If  the  response  to  the  linear  system  at  time  t  to  a  Dirac  6- 
function  applied  at  time  x[6(t-x)J  is  f(t,t),  the  superposition 
principle  [Equation  (6-1-1)]  and  the  sifting  integral  [Equation 
(6-1-3) ]  dictate  that  the  response  of  the  linear  system  to  the 
arbitrary  input  x(t)  is 

y(t)  »  |  x (t) f (t , t) dx  (6-1-4) 


The  function  f(t,x)  is  termed  the  time-variant  impulse  response. 
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Alternatively,  the  input  x  may  be  resolved  into  a  continuum  of 
complex  exponentials  using  the  Fourier  transform  pair 


x ( t )  -  |  X(f)exp{j2aft)df 
•00 

(6-1-5) 

X(f)  »  |  x(t)exp{-j2nft}dt 

(6-1-6) 

The  function  X(f)  is  called  the  amplitude  pectrum  of  the  input 
x(t) .  If  the  response  of  the  linear  system  at  time  t  to  the  com¬ 
plex  exponential  of  unit  amplitude  axp{ j2nft)  is  F(f ,t) ,  the  super¬ 
position  principle  [Equation  (6-1-1))  and  the  inverse  Fourier  trans¬ 
form  [Equation  (6-1-5) ]  dictate  that  the  response  of  the  linear 
system  to  the  arbitrary  input  x(t)  is 

y (t)  -  J  X(f)F(f,t)exp{j2xft}df  (6-1-7) 

The  function  F(f,t)  is  termed  the  time-variant  transfer  function. 

6.2  Forest  Transfer  Function 

The  time-variant  transfer  function  of  a  linear  system  repre¬ 
sents  the  response  of  that  linear  system  at  time  t  to  the  complex 
exponential  of  unit  amplitude  exp{  j2irft) .  Inasmuch  as  the  trans¬ 
mitting  antenna  of  the  anisotropic  half-space  model  of  the  forest 
is  taken  to  be  a  vertical  electric  (Hertzian)  dipole  having  a  time- 
harmonic  current  moment  Idl*exp{ jwt)  [refer  to  Section  5.1] ,  the 
mean  component  of  the  vertically-polarized  electric  field  found  in 
response  to  that  current-moment  [refer  to  Section  5.4],  may  be  con¬ 
sidered  to  be  the  forest  transfer  function  F(f,t)  after  normalizing 
the  expressions  for  the  mean  field  by  Idl.  Further,  as  a  consequence 
of  the  relationship  between  the  expressions  for  the  mean  field  and 
the  forest  transfer  function,  the  forest  transfer  function  F(f,t) 
may  be  considered  to  consist  of  three  components  corresponding  to 
the  direct  wave  [Fg(f,t)l,  the  reflected  wave  (Fr(f,t)),  and  the 
lateral  wave  [F£(f,t)]. 

Inasmuch  as  a  real  input  (excitation)  to  any  linear  system 
is  always  observed  to  produce  a  real  output  (response),  it  may 


be  shown  as  a  consequence  of  Equations  (6-1-6)  and  (6-1-7)  that 
the  real  part  of  F(f,t)  must  be  an  even  function  of  frequency, 
i.e. 

F* (f,t)  -  F‘(-f,t)  (6-2-1) 

and  the  imaginary  part  of  F(f,t)  must  be  an  odd  function  of  fre¬ 
quency,  i.e. 

F"(f,t)  -  -F"(-f,t)  (6-2-2) 

Further,  inasmuch  as  d-c  transmission  through  the  forest  is  not 
possible, 

F(0,t)  -  0  (6-2-3) 

These  relations  have  proved  helpful  in  evaluating  the  field  ex¬ 
pressions  of  Section  5.4  at  negative  frequencies  as  required  by 
Equation  (6-1-7) . 

It  is  apparent  from  Equation  (6-1-7)  that  the  forest  transfer 
function  F(f,t)  should  be  defined  over  the  entire  frequency  range 
from  minus  infinity  to  plus  infinity.  In  practice,  however,  the 
forest  transfer  function  need  only  be  defined  over  the  spectral 
range  occupied  by  the  input  signal  X(f).  This  is  fortunate  be¬ 
cause  the  forest  transfar  function  (the  normalized  vertically- 
polarized  electric  field)  is  expressed  oa  the  basis  of  a  high- 
frequency  asymptotic  evaluation  [refer  to  Section  53]  which  may 
give  incorrect  results  below  one  Megahertz. 

Because  each  of  the  three  forest  transfer-function  components 
may  be  considered  to  be  a  complex  phasor,  the  received  signal  will 
be  either  enhanced  or  degraded  depending  upon  their  relative  strength 
and  phase.  Their  relative  strengths  and  phases  are  clearly  fre¬ 
quency-dependent  [refer  to  Section  5.4]  and,  therefore,  the  forest 
transfer  function  is  frequency-selective. 

As  a  consequence  of  the  frequency-dependent  factor  (wy0/4Tr) 
appearing  in  Equations  (5-4-1)  and  (5-4-2),  the  effect  of  the 
direct-  and  reflected-wave  components  of  the  forest  transfer  func¬ 
tion  on  the  transmitted  signal  will  be  primarily  one  of  differen¬ 
tiation.  This  is  also  true  in  frae-space  transmission  as  can  be 
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clearly  seen  by  setting  et *  ez  *  1  in  Equation  (5-4-1).  Except  for 
its  exponential  factor  (which  is  primarily  one  of  phase  delay) , 
the  lateral-wave  component  of  the  forest  transfer  function  is 
nearly  frequency  independent.  This  means  that  if  an  equalizer 
having  a  response  inversely  proportional  to  frequency  (f“l)  is 
employed  at  the  receiver  to  ensure  distortionless  free-space  trans¬ 
mission,  direct-  and  reilected-wave  pulses  are  received  nearly 
distortion  free,  whereas  wide-band  (40C  Ifiz)  lateral-wave  pulses 
could  be  severely  distorted. 


6.3  Forest  Pulse  Response 

Consider  the  current  moment  of  the  electrically  small  (Hertz¬ 
ian)  dipole  transmitting  antenna  to  be  the  rectangular  r-f  pulse 

x (t )  -  p(t>  •  sin (2irfet)  (6-3-1) 


where  f0  is  the  r-f  carrier  frequency,  and  p(t)  is  the  pulse 
envelope  defined  here  to  be 


P(t) 


It |  <  L/2 
|t |  >  L/2 


(6-3-2) 


where  L  is  the  pulse  length  (duration) . 


The  amplitude  spectrum  of  the  pulse  envelope  and  of  the  rec¬ 
tangular  r-f  pulse  can  be  determined  from  Equation  (6-1-6)  to  be, 
respectively, 

P(f)  *  L  sinc(irLf)  (6-3-3) 

and 

X(f)  -  >*jI^sinc{TrL(f  +  f 0)  }  -  sinc{wL(f  -  f 0)  } |  (6-3-4) 

The  spectral  bandwidth  occupied  by  the  central  lobe  about  the 
carrier  frequency  f0  is  2/L. 

The  forest  pulse  response ,  if  defined  to  be  the  vertically- 
polarized  component  of  the  mean  electric  field  arising  in  response 
to  a  rectangular  r-f  input  pulse  of  current  moment,  can  be  obtain¬ 
ed  from  Equation  (6-1-7)  upon  the  substitution  of  Equation  (6-3-4) 
for  the  amplitude  spectrum  of  the  input  X(f)  and  the  use  of 
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Equations  (5-4-1),  (5-4-2)  and  (5-4-3)  for  the  forest  channel 
transfer  functions. 

€.4  Numerical  Evaluation 

The  evaluation  of  Equation  (6-1-7)  is  not  easily  effected 
analytically  due  to  the  complex  frequency  dependence  exhibited 
by  the  effective  dyadic  permittivity  of  the  forest  £.  As  a  con¬ 
sequence,  numerical  techniques  based  upon  the  fast  Fourier  trans¬ 
form  (FFT)  have  been  employed  to  expedite  its  evaluation.  A  re¬ 
view  of  the  important  properties  of  the  FFT  and  special  consider¬ 
ations  bearing  upon  its  application  to  the  numerical  evaluation  of 
Equation  (6-1-7)  is  presented  in  Appendix  E. 

To  exemplify  the  numerical  evaluation  of  the  forest  pulse 
response,  CyberCom  has  elected  to  consider  (1)  a  transmitting- 
antenna  current  moment  representable  by  a  rectangular  r-f  pulse; 

(2)  a  forest  transfer  function  developed  asymptotically  on  the 
basis  of  an  anisotropic  half-space  consisting  entirely  of  elec- 
trically-thin  leaves;  and  a  receiver  with  or  without  equalisa¬ 
tion  for  distortionless  pulse  transmission  in  free  space.  The 
channel  model  is  illustrated  in  Figure  6-1. 

The  density  of  the  leaves  representing  the  forest  has  been 
taken  as  200  leaves  per  cubic  meter  with  each  leaf  represented 
as  a  dielectric  disc  5  centimeters  in  radius  and  1  millimeter 
in  thickness;  the  leaves  have  been  assumed  to  be  randomly  posi¬ 
tioned  in  asimuth  and  location,  and  inclined  from  the  horisontal 
with  a  probability  density  function  uniform  over  the  range  0  to 
30  degrees.  The  transmitter  has  been  positioned  4  meters  0^) 
and  the  receiver  6  meters  (hr)  below  the  air-forest  interface; 
they  have  been  separated  by  1  kilometer  (p)  [refer  to  Figure  5-2]. 

The  forest  transfer  functions  corresponding  to  the  direct- 
wave  and  lateral-wave  components  have  been  determined  for  the 
postulated  forest  parameters  from  Equations  (5-4-1)  and  (5-4-3) , 
respectively,  and  are  shown  (magnitudes  only)  in  Figure  6-2  for 
both  unequalized  and  free-space  equalized  receivers.  The  nearly 
linear  frequency  response  of  the  unequalized  direct-wave  component 
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Figure  6-1:  Forest  Channel  Model 


Frequency  in  MHz 


apparent  below  600  Megahertz  will,  as  noted  previously,  lead  to 
differentiation  and  distortion  of  a  pulse  propagated  through  the 
forest.  The  transfer  function  of  the  direct-wave  component 
equalized  for  distortionless  free-space  transmission  is  nearly 
independent  of  frequency  below  600  Megahertz.  However,  the 
transfer  function  of  the  equalized  lateral-wave  component 
appears  now  to  behave  as  a  low-pass  filter.  A  comparison  of 
the  unequalized  and  equalized  forest  transfer  functions  suggests 
that,  for  distortionless  transmission  through  arbitrary  forests, 
adaptive  equalization  parametric  upon  path  length,  antenna  heights, 
frequency,  and  forest  susceptibility  may  be  required. 

The  rectangular  r-f  input  pulse  of  current  moment  x(t)  is 
shown  in  Figure  6-3 (a).  This  short  pulse,  only  5  nanoseconds 
long,  consists  of  3-cycles  of  a  600-Megahertz  r-f  carrie*  and 
occupies  a  spectral  bandwidth  (central  lobe)  of  400  Megahertz. 

The  unequalized  r-f  output  pulse  responses  y(t)  corresponding 
to  the  principal  components  of  the  forest  transfer  function  are 
shown  in  Figure  6-3 (b)  to  (e) .  These  responses  have  been  com¬ 
puted  for  the  forest  channel  model  of  Figure  6-1  using  a  256- 
point  (28-point)  FFT  having  a  time  resolution  of  0.20833... 
nanoseconds  and  a  frequency  resolution  of  53.33...  Kilohertz. 

The  direct-  and  reflected-  wave  pulses  are  distorted  more  than 
the  lateral-wave  pulse  for  reasons  previously  noted.  For  the 
equalized  case  shown  in  Figure  6-4,  the  direct-  and  reflected- 
wave  pulses  are  distorted  less  than  the  lateral-wave  pulse. 
Distortion  of  the- lateral-wave  pulse  is  primarily  attributable 
to  the  poor  low-frequency  response  of  the  asymptotic  model , 
although  some  distortion  is  also  attributable  to  the  equalizer. 
Although  the  pulse  distortion  introduced  by  the  forest  appears 
to  be  relatively  small  (at  least  for  the  forest  parameters 
employed  here),  multipath  and  intersymbol  interference  (ISI) 
could  be  substantial  for  multipath  transmissions. 

At  the  r-f  carrier  frequency  of  600  Megahertz,  the  am¬ 
plitude  of  the  reflected-wave  relative  to  that  of  the  direct- 
wave  is  0.986  (0.12dB  down)  and  its  phase  almost  directly 
opposite.  Although  the  differential  phase  shift  between 
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direct-  and  reflected  "waves  attributable  to  path  geometry 


-  72p^ht^g  (degreea)  (6-4-1) 

is  only  about  35  degrees,  the  phase  shift  imparted  upon  reflection 
from  the  air- forest  interface  is  about  172  degrees.  Thus,  the  rel¬ 
ative  amplitude  of  the  direct-plus-reflected  wave  is  only  0.107  or 
about  19  dB  below  the  direct.  In  the  configuration  considered, 
which  may  not  be  realistic  because  the  effects  of  tree  trunks  and 
branches  are  not  yet  included,  although  the  amplitude  of  the  lateral- 
wave  relative  to  that  of  the  direct-wave  is  only  0.22  (13  dB  down), 
it  is  approximately  6  dB  above  the  direct-plus-reflected  component 
propagating  within  the  forest.  The  differential  time  of  arrival 
(DTOA)  between  pulses  arriving  via  the  direct  and  lateral  waves 
is  approximately 

(PXi/2  -  (ht+hr)  i£J"]300  *  5.9  nanoseconds  (6-4-2) 
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7.0  Non-Coherent  Forest  Scattering 

In  the  preceding  three  sections  propagation  through  the  for¬ 
est  has  been  described  neglecting  fluctuations  about  the  mean  field. 
This  section  will  show  that  these  become  important  as  the  frequency 
and  propagation  distance  increase.  In  this  section,  the  effect  of 
the  fluctuations  will  be  taken  into  account  for  one  specialized 
forest  configuration.  The  results  will  give  an  indication  of  the 
importance  of  the  fluctuating  field  in  the  forest  environment. 

In  order  to  keep  the  analysis  tractable,  the  scalar  case  of 
propagation  through  an  infinite  medium  of  vertical  trunks  will  be 
considered.  The  trunks  will  be  assumed  to  be  identical  vertical 
circular  cylinders  of  radius  a  and  relative  dielectric  constant 
e^.  In  addition,  their  diameter  will  be  chosen  small  enough  so 
that  the  Rayleigh  criterion  is  satisfied.  Although  the  Rayleigh 
assumption  is  somewhat  restrictive,  physically  meaningful  results 
can  still  be  attained  over  a  portion  if  the  frequency  range  of 
interest.  The  field  is  excited  by  an  electric  line  source  of  unit 
magnitude  which  is  parallel  to  the  trunks.  As  has  been  pointed 
out  in  Section  2.2.1,  the  assumed  source  and  scatterers  give  rise 
to  an  electric  field  vector  wholly  in  the  z  direction.  The  pro¬ 
blem  can  be  exactly  scalarized  in  terms  of  this  component. 

The  strength  of  the  random  component  will  be  examined  by 
determining  the  intensity  of  the  field 

I(xt)  =  <4'(xt)’i'*(xt)>  (7-1) 

and  relating  it  to'  the  intensity  of  the  mean  component 

I«(xt)  =  <^(xt)><4'*(xt)>  (7-2) 

For  two-dimensional  Rayleigh  scatterers,  the  intensity  satisfies 
Equation  (2-2-50)  with  xt*xt.  The  equation  is 

I(xt)  *  jG(xt)  |2  -  ejdst|G(xt  -st) |2I(st)  (7-3) 

where  &»k£|a|2p  with  a  being  the  polarizability  constant  appearing 
in  Equation  (2-2-21) .  For  a  circular  cylinder,  the  polarizability 
is  a*-XpTra2  with  X^*^-!.  In  obtaining  Equation  (7-3)  from 


Equation  (2-2-50)  it  has  been  assumed  that  the  trunk  medium  is 
unbounded  (s*»)  and  the  density  p  is  constant.  These  homogeneity 
assumptions  lead  to  a  translationally  invariant  Green's  function: 

G(*t'-t*  "  G*2it“-t*  (7-4) 

The  mean  Green's  function  satisfies  Equation  (2-2-47).  In 
the  Rayleigh  regime  this  equation  reduces  to 

[V|  +  k2(l  +  pa)]G(xt)  *  6(xt)  (7-5) 

The  solution  is  given  in  terms  of  the  zeroth  order  Hankel  function 
of  second  kind.  It  is 

G(xt)  -  J  H !2) (k)  ,  (7-6) 

and 

k  «  k0/l  ♦  pa  (7-7) 


Before  proceeding  it  should  be  pointed  out  that  the  intensity 
Equation  (7-3)  was  derived  under  the  assumption  of  small  fractional 
area  (two-dimensional  analog  of  small  fractional  volume  used  for 
three-dimensional  problems) .  Analytically,  this  condition  states 
that  ira2p<<l.  For  a  forest  having  a  »  10cm  and  p  *  .1/m,  the  con¬ 
dition  is  met,  since  7ra2p  «  0.003  <<  1.  Since  pa  «  X^(7ra2p)  and 
| | =40 ,  one  can  also  conclude  pa<<l. 

The  exact  solution  to  Equation  (7-3)  can  be  obtained  by  con¬ 
volution  techniques;  however,  the  inverse  transform  involves  a 
double  integral.  J£o  simplify  the  form  of  the  solution,  the  Green's 
function  in  Equation  (7-3)  is  replaced  by  its  far-field  approxima¬ 
tion.  One  obtains 


C(xt) 


j  -j  (KXt  -  tt/4) 
2/2TTKXt  e 


(7-8) 


where  |  kx^.  |  > >1 .  Further  simplification  is  obtained  by  using  the 
fact  that 


k  z  k0  (1  +  p  a/2) 


(7-9) 
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in  the  asymptotic  expression  for  G.  This  result  is: 


1  -j  (k.x«.  -k.p  a/2  -  ir/4) 

G(xt)  -  —  ZZZ '  '  e  (7-10) 

c  2*'2irk0xt 

Now  using  Equation  (7-10)  in  Equation  (7-3) ,  the  approximate 
intensity  equation  is 

I(xt)  -  K(xt)  -  ejdstK(xt  -it)I(*t)  (7-11) 


where 


k0pa"xt 

K(x.  )  -=  5 - 

8irk0xt 


(7-12) 


The  solution  to  this  equation  can  be  obtained  in  terms  of  the 
single  quadrature  since  now  the  Fourier  transform  of  K(xfc)  can 
be  explicitly  evaluated. 


Before  solving  Equation  (7-11),  the  physical  restrictions 
imposed  on  the  solution  when  the  asymptotic  form  of  the  Green's 
function  is  used  are  addressed.  First,  the  term  |G(xt)|2  appear¬ 
ing  on  the  right  hand  side  of  Equation  (7-3)  has  been  replaced 
by  its  asymptotic  approximation.  This  requires  that  x|xel>>l  or 
that  the  receiving  antenna  be  in  the  far  field  of  the  transmitter. 
Second,  the  Green's  function  in  the  integrand  of  Equation  (7-3) 
has  been  replaced  by  its  asymptotic  expansion.  This  requires  that 
x|x-s|>>1  or  that  the  scatterers  be  in  the  far  field  of  the  obser 
vation  point.  Thus  scatterers  that  are  in  the  near  field  of  the 
observation  point  are  not  treated  correctly.  This  inaccuracy  is 
small,  however,  if  the  far  field  distance  is  not  too  large.  In 
this  case,  there  are  many  more  scatterers  in  the  far  field  of  the 
observation  point  than  in  the  near  field. 


Integral  Equation  (7-11)  can  be  solved  exactly  using  Fourier 
transform  theory.  If  the  Fourier  transform  of  I(x)  is  denoted  by 
I(k)  where 


I(x*) 


dk.  I (k. )e 


•x. 


t  -t 


(7-13) 


then  the  transform  of  Equation  (7-11)  becomes: 

I(kt)  -  K(kt)  +  lK(kt)l(kt)  ,  g  -  (2ir)  2g 


where 


(7-14) 


K(kt)  -  K(kt)  -  (4ir')*jt#  *  +  (k.a**p) » 

with  kt *  |kt|.  Solving  Equation  (7-14)  for  i(kt)  gives 
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Kkt) 


K(kt) 

1  -  3K(kt) 


(7-16) 


Now  using  this  result  in  Equation  (7-13) ,  the  intensity  is 


K(kfc) 


f  RlKt 

(x. )  -  dk  - ~ 

r  J  z  1  -  6K( 


(7-17) 


8K(kt) 


Since  I(kt)  only  depends  on  the  absolute  value  of  k^,  the  angular 
depend'  'ce  in  Equation  (7-17)  can  be  integrated  explicitly.  Let 
6  be  the  angle  between  kt  and  xt;  then 


*t*-t  *  fctxtcos0 


(7-18) 


and  dkfc  «  ktdktd8.  Equation  (7-17)  becomes 


rm  r 2tt  k*K(k.)  -jk.  x  cos6 
(x.)  -  I(xt)»  dk.  de  -  --_■■■-  e  z  z 
*  L  L  1  -  BK(k. ) 


(7-19) 


BK(kt) 


1  «st>  -  sr  IX 


htJ0(ktxt) 


(7-20) 


The  expression  for  the  intensity  given  in  Equation  (7-20)  is 
a  function  of  x^,  k0 ,  a  and  3.  This  dependence  on  four  parameters 
can  L>e  reduced  to  two  parameters  by  appropriate  scaling  and  nor¬ 
malization.  For  this  purpose,  the  two-dimensional  scattering 

cross  section  o_  and  total  cross  section  a .  are  introduced.  These 

s  t 

quantities  are  related  to  the  transform  of  the  two-dimensional 
transition  operator  as  follows  (Appendix  F) 
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°s  “  TT  J  de|t(k0o,k,i) 


at  ■  -  -  Im  t (k0o,k.i) 


For  Rayleigh  scatterers,  t(k0o,k0i)  becomes 

t (k0o,k.i) 


a  v2 

r  *• 


<2w> 

Using  this  in  the  expressions  for  and  ofc  gives 

k0|a|* 

°s  * 


(7-22) 


(7-23) 
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ot  «  -  k#a"  (7-25) 

The  albedo,  W0 ,  and  the  optical  distance,  t,  can  be  introduced 
using  the  cross  sections  just  defined.  They  are: 


W0 


t  -  potxt 


(7-26) 


The  albedo  is  a  ratio  of  the  scattered  power  to  the  total  power, 
with  the  total  power  defined  as  the  scattered  power  plus  the  power 
absorbed  by  the  particle  (Appendix  F) .  The  albedo  lies  between  0 
and  1.  For  W„  ■  0  there  is  no  scattering  (all  absorption)  and  for 
W0  «  1  there  is  no  absorption  (all  scattering) . 


The  optical  distance  is  a  measure  of  the  attenuation  of  the 
mean  intensity,  I0(xt).  Thus,  in  the  far  field  of  the  source 


I*  <xt) 


|G(xt)|2  -  * 


-potxt 


8trk0xt 


(7-27) 


where  Equations  (7-10)  and  (7-25)  have  been  used.  From  Equation 
(7-27)  it  is  seen  that  when  T*potxt»l,  the  exponential  has  an 
argument  of  -1.  Now  Equation  (7-20)  con  be  normalized.  Let  q  = 
kt/pot  and  define  a  normalized  intensity 


*i 


m 


I(xt)  -  I(xt)/I#(xt) 


(7-28) 


The  expression  for  the  normalized  intensity  is  then  given  by 

-  f*  dqqJ.  (qx) 

I(x)  -  xeT  |  -  (7-29) 

'0  /qs  +  I  -  W„ 

Thus  the  normalized  intensity  T  is  only  a  function  of  x  and  W0 . 

The  expression  of  T(x)  given  in  Equation  (7-29)  is  exact  but 
slowly  converging.  As  q  becomes  large ,  the  integrand  is  propor¬ 
tional  to  cos(qx)//q.  This  convergence  makes  numerical  evaluation 
difficult.  Fortunately,  by  deforming  the  integral  into  the  com¬ 
plex  plane,  a  more  rapidly  converging  representation  can  be  obtained 
(see  Appendix  G) .  The  result  is 

T(t)  »  2TW.eTK0(T/l -Ml)  eTf  dpp2  *■? (T-A±£lL  (7-30) 

"  J0  p1  +  w. 

The  integrand  is  now  exponentially  decaying  for  large  p.  As  a 
check,  it  is  interesting  to  note  if  W„*0,  then  I(x)  «1.  This 
means  that  if  there  is  no  scattering,  the  mean  intensity  is  exactly 
equal  to  the  total  intensity. 

The  coherent  and  noncoherent  intensities  have  been  computed. 
Some  representative  results  of  the  computations  are  shown  in  Figure 
7-1.  Because  of  the  need  to  satisfy  the  Rayleigh  criterion, 
k0/e  a  <<  1,  the  small  trunk  radius  of  0.01  meters  was  chosen. 

Using  this  trunk  radius  and  Cyber Com  permittivity  modal  III,  the 
Rayleigh  number  was  computed  to  be  0.4  at  300  MHz.  It  is  also 
seen  from  Figure  7-1  that  a  density  of  0.1/m2  has  been  used.  This 
represents,  on  the  average,  one  tree  in  every  10  square  meter  area. 

It  can  be  seen  that  the  coherent  intensities  (related  to  the 
powers  by  a  constant)  decrease  almost  linearly  with  distance  as 
expected  in  the  assumed  two-dimensional  medium.  The  noncoherent 
intensities,  however,  fall  off  much  more  slowly  with  distance, 
thereby  becoming  an  increasing  fraction  of  the  total  intensity. 

At  30C  MHz,  the  noncoherent  intensity  becomes  greater  than  the 
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coherent  intensity  beyond  about  230  meters.  This  behavior  is  simi¬ 
lar  to  that  in  Figure  8  of  Schwering  et  al  [28a] ,  which  also  uses 
transport  theory.  If  the  total  intensity  is  replotted  on  a  linear 
distance  scale  as  in  Schwering  et  al,  it  is  seen  to  fall  off  at  a 
decreasing  rate  with  distance. 

In  Figure  7-2  are  plotted  the  normalised  intensities  (total/ 
coherent) .  An  examination  of  the  curves  shows  that  the  normalized 
intensity  increases  with  distance.  At  a  frequency  of  100  MHz  or 
200  MHz,  the  normalized  intensity  stays  fairly  close  to  unity  for 
distances  as  great  as  one  kilometer.  This  means  that  the  random 
component  of  the  intensity  is  relatively  small  for  frequencies  be¬ 
low  200  MHz  and  distances  as  large  as  1  kilometer.  The  curve  for 
300  MHz  increases  rapidly  with  distance,  however,  and  reveals  that 
the  random  component  will  be  important  for  distances  in  excess  of 
200  meters. 

It  appears  that  the  random  component  of  the  intensity  increases 
with  frequency  and  propagation  distance.  Although  the  computed  val¬ 
ues  presented  in  Figure  7-1  are  interesting,  to  compute  the  inten¬ 
sity  for  a  forest  with  larger  trunks,  the  Rayleigh  or  low-frequency 
approximation  can  not  be  employed.  The  correlation  equation  must 
be  used  with  the  transition  kernel  for  resonant  trunks.  This  has 
yet  to  be  done. 
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8.0  Conclusions 

Cyber Com  has  developed  a  stochastic  radiowave  propagation 
model  useful  for  assessing  the  effects  of  forests  and  other  vege¬ 
tation  upon  radio  communication  systems  operating  in  the  200  - 
2000  Megahertz  frequency  band.  The  model  considers  the  forest  to 
be  representable  as  a  time-invariant  ensemble  of  randomly  posi¬ 
tioned  and  oriented,  discrete,  canonical  scatterers.  Tree  trunks 
are  modeled  as  infinitely-long,  circular,  lossy-dielectric  cylin¬ 
ders;  branches  as  finitely-long,  circular,  lossy-dielectric  cylin¬ 
ders;  and  leaves  as  flat,  circular,  lossy-dielectric  discs.  Math¬ 
ematical  expressions  characterizing  the  scattering  properties  of 
these  canonical  scatterers  in  terms  of  their  dyadic  scattering 
amplitudes  have  been  developed  and  integrated  into  the  multiple¬ 
scattering  theory  of  Foldy  [53]  to  develop  equations  describing 
the  behavior  of  a  propagating  radiowave.  In  the  case  of  tree 
trunks  and  leaves,  the  model  is  suitable  even  in  the  difficult  re¬ 
sonant  scattering  region  where  the  physical  dimensions  of  the 
scatterer  are  comparable  to  the  radio  wavelength  within.  However, 
in  the  case  of  branches,  additional  study  would  be  required  to 
extend  the  model  into  the  resonant  regime. 

The  propagating  radiowave  may  be  considered  to  consist  of  two 
components  -  a  mean  (or  coherent)  component  and  a  random  (or  non¬ 
coherent)  component.  CyberCom  has  shown  that,  so  far  as  the  mean 
radiowave  component  is  concerned,  the  ensemble  of  discrete  canon¬ 
ical  scatterers  representing  the  forest  can  be  replaced  by  an  equi¬ 
valent  continuous  medium  characterized  by  an  effective  dyadic  sus¬ 
ceptibility  Simple  mathematical  expressions  have  been  developed 
relating  X  to  such  biophysical  parameters  of  the  forest  as  tree 
trunk,  branch  and  leaf  dimensions,  permittivity,  number  density, 
and  orientation  statistics.  Expressions  relating  the  specific 
attenuation  of  a  radiowave  propagating  within  an  unbounded  forest 
of  tree  trunks,  branches  or  leaves  to  salient  biophysical  para¬ 
meters  of  the  forest  have  been  developed  and  shown  to  agree  favor¬ 
ably  with  experimental  data.  Utilization  of  these  expressions  in 
radio  communication  systems  design  still  requires,  however,  a 


global  statistical  classification  and  quantification  of  these  bio¬ 
physical  forest  parameters. 

The  effective  dyadic  susceptibility  has  been  used  by  CyberCom 
in  a  half-space  representation  of  the  forest  to  extend  an  earlier 
isotropic  model  proposed  by  Sachs  and  Wyatt  [101]  and  studied  by 
Tamir  [107] .  This  newer  anisotropic  model  encompasses  the  lateral- 
wave  contributions  observed  experimentally  on  long-distance  paths 
above  200  Megahertz.  For  generality,  it  should  be  refined  by  con¬ 
sideration  of  the  ground  plane. 

Using  linear  system  theory,  CyberCom  has  extended  the  harmonic 
anisotropic  half-space  model  of  the  forest  to  the  transmission  of 
arbitrarily -modulated  waveforms.  This  extension  has  been  used  to 
show  the  pulse  distortion  arising  in  the  direct-,  reflected-  and 
lateral-wave  components  of  the  mean  scattered  field,  and  also  in 
the  total  mean  field.  This  extension  of  the  model  is  especially 
important  if  the  model  is  to  be  used  in  the  analysis  and  evaluation 
of  digital ,  spread-spectrum  radio  systems . 

CyberCom  has  also  investigated  the  random  (or  non-coherent) 
scattered  field  propagating  within  an  unbounded  forest  of  electri- 
cally-thin,  inf initely-long,  lossy-dielectric,  parallel  tree  trunks. 
For  this  simple  forest  model,  CyberCom  has  shown  that  the  intensity 
of  the  non-coherent  scattered  field  increases  relative  to  that  of 
the  coherent  (mean)  field  with  increasing  distance  from  the  source 
and  with  increasing  carrier  frequency.  This  work  should  be  extended 
so  far  as  possible  towards  fully  realistic  models. 

In  summary,  the  homogenous  isotropic,  refracting  slab  model  of 
a  forest  has  been  replaced  by  CyberCom  with  an  inhomogenous ,  aniso¬ 
tropic,  scattering  ensemble  of  trunks,  branches,  and  leaves.  In 
consequence : 

1.  The  lateral-wave  contribution  has  been  found  even  above 
200  Megahertz. 

2.  Preliminary  validation  of  the  model  has  shown  rough  agree¬ 
ment  with  experiments. 

3.  Results  have  been  obtained  for  narrow  pulse  transmission 
at  600  MHz. 
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4.  A  preliminary  study  (as  well  as  field  reports)  has  empha¬ 
sized  the  importance  of  the  incoherent  component  of  the  transmitted 
field. 

As  the  CyberCom  approach  is  ambitious  in  both  scope  and  depth, 
the  following  remain  to  be  done: 

a.  The  techniques  already  developed  must  be  exercised  to 
determine  the  relative  importance  of  forest  components  and  the 
effects  of  varying  critical  parameters.  For  example,  if  the  con¬ 
tribution  of  branches  is  major,  the  difficult  characterization  of 
the  scattering  properties  should  be  advanced  in  the  resonant  region. 

b.  The  effects  of  antenna  directivity,  the  ground,  and  ter¬ 
minal  movement  should  be  incorporated  into  the  model. 

c.  The  difficult  transport  theory  for  the  important  non¬ 
coherent  scattered  wave  should  be  developed  as  far  as  practical. 

d.  Forest  studies  must  be  pursued  to  quantify  important 
biophysical  parameters  in  areas  of  interest. 
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APPENDIX  A 


Susceptibilities  of  Wood  and  Leaves 

The  complex  susceptibilities  of  green  wood  and  leaves  required 
for  the  modeling  approach  of  CyberCom  have  been  investigated.  Meas¬ 
urements  by  Broadhurst  [1]  of  the  National  Bureau  of  Standards  ap¬ 
pear  to  be  the  best  available.  They  are  summarized  in  Figure  A-l. 

The  spread  of  the  results  for  the  leaves  represents  the  effects  of 
different  percentages  of  moisture  (30-78%)  and  for  the  branches, 
different  orientations  relative  to  the  branch  axis.  The  curves 
show  that  the  real  parts  of  the  susceptibilities  x&'  are  nearly 
constant  with  frequency;  the  imaginary  parts  of  the  susceptibilities 
X^"  are  inversely  proportional  to  frequency  at  the  lower  frequencies 
of  interest,  and  rise  with  frequency  at  the  higher;  the  results  for 
bamboo  leaves  and  for  tulip  tree  leaves  and  branches  overlap. 

CyberCom  has  developed  and  employed  three  simple  models  for 
the  susceptibility  of  green  wood  and  leaves,  summarized  in  Table 
A-l.  Model  I  describes  a  frequency-independent  susceptibility  which 
is  convenient  for  ascertaining  the  effects  of  scatterer  geometry  on 
the  scattered  wave.  Model  II  describes  a  susceptibility  character¬ 
ized  by  a  real  part  which  is  frequency-independent  and  an  imaginary 
part  derived  from  a  constant  conductivity  of  0.1  mho/meter.  Model 
III  best  reflects  the  typical  behavior  of  measured  data  and  includes 
both  a  frequency-independent  real  part  and  an  imaginary  part  related 
to  both  conduction  and  relaxation  losses.  It  is  shown  in  Figure  A-l. 

The  imaginary  parts  of  the  susceptibilities  (the  loss  factors) 
are  shown  in  Figure  A-2  for  the  three  CyberCom  models.  The  recovery 
of  the  Model  III  loss  factor  from  its  minimum  near  1  GHz  is  clearly 
evident,  and  may  be  attributed  to  increased  relaxation  losses  as 
the  resonant  frequency  of  the  water  molecule  (approximately  22  GHz) 
is  approached. 


TABLE  A-l  Susceptibility  of  Wood  and  Leaves  -  CyberCom  Models 


CyberCom 

l+Xj-JXj  -  e£“EA'-je£" 

Model 

KWLUtl 

Xa"*E£"  (Loss  Factor) 

I 

39 

10 

II 

39 

1.8/f_„  (a  «=  0.1  mho/m) 

GHZ 

III 

39 

(1.5/fGBi)  ♦  (2fGH2)/U+(fGHz/20)M 

Appendix  B 


Scattering  Amplitude-Transition  Operator  Relationship 
for  Two-Dimensional  Scatterers 


In  this  Appendix  the  relationship  between  the  scattering  ampli¬ 
tude  and  the  transition  operator  will  be  obtained  for  two-dimensional 
scatterers.  This  will  be  done  for  both  the  vector  case  (oblique 
incidence)  and  the  scalar  case  (normal  incidence) .  The  methodology 
used  will  parallel  that  of  Lang  [57]  who  derived  the  relationship  for 
three-dimensional  scatterers.' 

Assume  a  dyadic  plane  wave  is  incident  on  a  two-dimensional 
scatterer  as  shown  in  Figure  2-4.  The  incident  wave  is  of  the  form 


io(x)  =  (I-ii)e"j  (— t— t+Jczz) 


(B-l) 


where  i  is  in  the  direction  of  incidence,  i.e. 


i  -  <*t  +  kzr>/ko 


with 


+  ' 


kt  =  I  St 


(B-2) 


(B-3) 


The  dyadic  scattered  field  is  given  by 

e  (x)  *  /  dx'  G(x-x')  •  /  dx"  t(x',x") 
“8  *  ™"  ■“  *  **  ■r" 

•  e  (x") 

*o  — 

where  the  dyadic  Green's  function  G  satisfies 

as 

L  •  G(x)  -  16  (x) 


(B--4) 


(B-5) 


with  L  given  by  equation  (2-3-5)  and  t  the  kernel  of  the  transition 
operator . 


B-l 


Since  the  scatterers  are  two-dimensional ,  they  are  invariant  in 


the  z  direction,  i.e.. 


t(x',x")  =  t(x£,x£,z'-z") 

Putting  equation  (B-6)  into  equation  (B-4)  and  introducing  the 


(B-6) 


Fourier  transform  of  the  transition  operator  with  respect  to  z  yields 
gs(x)  *  /  dx^dz'  G(xt~x^, z-z' ) •  /  dx£  t ( x ^ , x^ ,  k2) 

•  (I-ii)e~^  t'— t+ks2'  *  (B-7) 


where 


t^t'St'V  *  ^  dz  I (xt/x^»*)e+:’kzz 


(B —  8 ) 


Now  introducing  the  Fourier  transform  of  the  dyadic  Green's  function 
with  respect  to  z  gives 

fig(x)  -  /  dx^.dx£  G(xt-x^.,kz)  •  t(x£,x£,kz> 

•  (I-ii)e^<it-5i+kzz)  (B'9) 


where 

G(xt,k2)  =  /  dz  G(xt,z)ejkzz  .  (B-10) 

Before  proceeding,  the  transformed  Green's  function  must 
be  found.  This  is  accomplished  by  transforming  the  well-known 
expression  for  the  three-dimensional  dyadic  Green's  function 
[Van  Bladel,  81].  The  three-dimensional  dyadic  Green's  function  is 


G (x)  -  (I  -  ——tt’)  g  ( x)  (B-ll) 

«  =  k  2 
o 

where  g(x)  is  the  scalar  Green's  function; 

9<2)  -  ,  r  -  |x|  (>-12) 

This  scalar  Green's  function  satisfies  the  scalar  wave  equation 

(V2+k2)g(x)  *  -  6  (x)  .  ( B— 1 3 ) 


(B— 1 3 ) 


B-2 


To  find  the  transformed  Green's  function,  first,  the  trans¬ 
formed  scalar  Green's  function  is  found.  This  is  done  by  Fourier- 


transforming  equation  (B-13)  and  then  solving  the  resulting 
two-dimensional  equation.  The  result  is 

g(x)  -  -  i  /  H^21  (kt|xt| dkz 


(B-14) 


where 


Kt  - 
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z  o 


(B-15) 


Now  putting  equation  (B-14)  in  equation  (B-ll) ,  bringing  the  dyadic 
portion  of  the  operator  through  the  integral  and  finally  using 


equation  (B-10)  yields 

A  A 

£<*t'V  -  -  + 

Ko 


(B-16) 


where  V  «  V  -  jk  z°. 

t  a 

To  proceed  with  developing  the  relationship  between  the  scatter¬ 
ing  amplitude  and  the  transition  operator,  an  asymptotic  expression 
for  e  (x)  will  be  needed.  An  examination  of  equation  (B-9)  shows 
that  this  in  turn  requires  a  large  |xt|  approximation  for 

A 

t'kz) *  Assuming  |xt|  >>  |x£|  and  using  the  asymptotic 


expression  for  the  Hankel  function  results  in 


H^1)(ktlxt-Xtl)  ~ 


nktxt 


rj(ktXt-ktOt«X^-TT/4) 


(B-17) 


where  ofc  ■  xt/|xt|.  Now  substituting  this  result  in  equation 
(B-16)  and  keeping  highest-order  terms  yields 

G(xt-x;,k2>  -  -jltggi.  (ktxt-ktot-x;-x/4) 

/filc^xt 


(B-18) 
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*  M 


where  o 


o.  +  oz 
— t  z— 


with  o. 


k2/ko* 


The  filial  form  of  the  asymptotic  expression  for  the  scattered 
field  e  (x)  is  obtained  by  substituting  equation  (B-18)  into 
equation  (B-9) .  The  result  is 

es00  -  <I~00>  ,^kt°t>ktit?kx) '  (ktxt+ks*Hr/4)  (B-19) 

t  t 


where  i_t  is  the  transverse  portion  of  i  (ifc  *  kt/kQ)  and  the 
two-dimensional  transformed  transition  operator  is  given  by 
ifc^ik,)  =  /  dxtdx£  i(xt,x^;k2) e^  (~t*— t”— t*— t)  (B-20) 

The  desired  relationship  can  be  established  by  comparing  the 

definition  of  scattering  amplitude  with  equation  (B-19) .  The  two- 

dimensional  dyadic  scattering  amplitude  is  defined  as 

-j  (k.x.+k  z) 

S  (x)  -  f (o. ,i. ;k  )  2—  z  *  Z  .  (B-21) 

-s  -  X.  — X,  Z  r— 

Comparing  equations  (B-19)  and  (B-21) ,  one  obtains 

!(Ot,it;kz)  -  Y(I-Oq)-t(ktot,ktit;kz)«(I-ii)  (B-22) 

where 


Y 


[27r3/kt] 


1/2 


e-j7F/4 


(B-23) 


Before  concluding  this  appendix  a  similar  relationship  will 
be  derived  for  the  scalar  problem.  Because  the  scalar  problem  is 
a  special  case  of  the  two-dimensional  vector  work  just  presented, 
those  results  can  be  used. 

The  scalar  problem  restricts  attention  to  propagation  normal 
to  the  tree  trunks  and  to  sources  that  only  excite  vertical  waves. 
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The  vector  problem  reduces  to  the  scalar  problem  as  follows: 

>*-°  *  2n<^ >fi 

(B-24) 

g(xt»Jc2)  *£  *  2tt  g(xt)  6(kz) 

If  one  dots  equation  (2-4-4)  by  z°  and  introduces  the  scalar 
specializations  as  given  by  equation  (B-24) ,  one  obtains  the  scalar 
mean-wave  equation  as  given  by  equation  (2-2-4)  with 

t(xt,x^)  *  z°* t (xt,x£;0) »z°  .  (B-25) 

In  addition,  an  examination  of  the  definition  of  both  the  vector 
and  scalar  scattering  amplitudes  as  given  in  equations  (2-2-10) 
and  (2-2-23)  shows  that 

f(o,i)  =  z°*I(o,i,0) .z°  .  (B-26) 

The  desired  relationship  between  the  scalar  scattering  amplitude  f 
and  the  transformed  scalar  transition  operator  t  is  obtained  by 
dotting  equation  (B-20)  from  the  left  and  the  right  by  z°,  and 
then  using  equation  (B-26)  and  the  transform  version  of  equation 
(B-25) .  The  result  is 

f (o,i)  =  6  t  <k0o,k0i> 


where  6  is  given  by  equation  (B-23). 


Appendix  C 
Correlation  Equation 


Ii.  this  Appendix  the  approximate  equation  for  the  correlation 
equation  will  be  derived.  The  methodology  employed  in  the  deri¬ 
vation  will  parallel  the  Foldy  development  for  the  mean  equation. 
The  starting  point  is  Equation  (2-2-29)  for  the  scalar  field 
quantity  .  Multiplying  this  equation  through  from  the  left  by 
L”1  gives  an  expression  for  V  in  terms  of  the  incident  and  scat¬ 
tered  fields.  It  is 


l  L“1T,¥(^) 
j-1  3 


(C-l) 


Here  the  observation  point  for  the  field  ¥  is  x,  i.e.,  Y^'Ffx). 
An  equation  for  the  conjugate  field  - V* (x)  at  point  x  can  be 
obtained  from  Equation  (C-l) .  The  result  is 


f*  *  ¥*  + 


f  l-wi>* 

j-i  j 


<02  > 


In  the  remainder  of  the  appendix,  the  dependence  of  quantities  on 
x  or  x  will  be  denoted  by  whether  the  symbols  are  starred  (con¬ 
jugate)  or  not. 


Now  forming  the  product  and  averaging  yields 
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1  1  1j-l  3  j-1 
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j-1 


N 

l 
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The  double  sum  appearing  in  the  above  equation  can  be  broken  into 
like  terms  (i*j)  and  unlike  terms  (i/j).  Doing  this,  the  follow¬ 
ing  expression  results: 


<¥¥*>  -  ¥*  +  4'.  I  L +  9  L“1<T.y^)>4'J 

1  1  1  j-1  3  j-1  3  1 

+  l  L“1l"1*<T.T}>¥^4'^)* 
j-1  *  j 


N  N 


+  l  l  l“1L~1*<T.T{>4'^4'*1c)  *  (04) 

j-1  k-1  3  K 


To  obtain  an  equation  for  the  correlation  function  a  closure 
approximation  is  employed.  The  approximation  to  be  used  is 

^(j)  *  <'ji> 


y(j)y*(j)s: 

y(j)y*(k)  2  <^><y*>  jjlfc 


(05) 


This  is  a  generalization  of  Foldy's  closure  approximation  used  to 
derive  the  mean  equation.  Using  this  approximation  in  Equation 
(C-4)  gives: 


<¥¥*> 


N  ? 
v  ,  -1. 


y  y*  +  I  l"a<t. 

IX  1  j-l  3 
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Since  the  particles  are  identically  distributed,  the  statistics 
for  each  particle  are  the  same,  i.e.. 


I  <T ■ >  -  N<T> 
j-1  3 

N 

T  <T.T*>  «  N<TT*> 
j-1  3  3 


(07) 
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Employing  these  simplifications  in  Equation  (C-6)  results  in 

-I 

<yy*>  «  +  n^l  A  <T*xf*> 

* 

+  NL"1<TX,r>'i'*  +  NL“1L“1<TT*><Y¥*> 

* 

+  (Na-N)L“1<T><y>L_1<T*x,r*>  (C-9) 

This  is  the  desired  equation  for  the  correlation  function  <¥¥*>. 

It  can  be  simplified  substantially  by  making  use  of  the  me an -wave 
equation. 

The  mean-wave  equation  is  given  in  the  text  by  Equation 
(2-2-48).  It  is 

£<lp>  ■  g  (C-10) 

where 

£  -  L  -  N<T>  (C-ll) 

In  order  to  simplify  Equation  (C-9) ,  the  mean  equation  is  put  in 
a  slightly  different  form.  With  this  purpose  in  mind,  multiply 
Equation  (C-10)  by  L*1  from  the  left.  This  gives 

<ip>  *  +  nl~1<T><'V>  (C-12 ) 

Now  using  this  result.  Equation  (C-9)  becomes 

* 

<yy*>  *  <¥><¥*>  +  (N2-N)L"1<T><'l'>L”1<TXf*>  (C-13) 

Multiplying  from  the  left  by  ££*,  using  the  fact  that  £L_1  and 
neglecting  the  N2  terms  in  Equation  (C-13)  yields 

££*^ij>*>  -  N<TT*Xi|>iJ>*>  ■  gg*  (C-14) 

The  above  steps  can  be  justified  by  scaling  the  problems  with  re¬ 
spect  to  the  characterstic  size  of  the  particle,  then  N  is  replaced 
by  the  fractional  volume  which  is  Bmall. 

The  final  form  of  the  correlation  equation  as  it  appears  in 
Equation  (2-2-44)  is  obtained  by  writing  the  spatial  average  of 
TT*  out  explicitly  and  by  using  Equation  (2-2-38). 
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Dyadic  Notation 

According  to  Lindell  [124],  dyadic  formalism  is  the  most  suit¬ 
able  notation  for  linear  vector  functions  when  applying  the  gibbsian 
vector  notation  to  electromagnetic  fields.  Nevertheless,  although 
conceived  by  J.  Willard  Gibbs  nearly  one-hundred  years  ago  (1884), 
dyadic  formalism  has  been  introduced  into  electromagnetics  only  com- 
paritively  recently,  and  is  not  generally  included  in  most  engineer¬ 
ing  curricula.  Thus,  for  the  convenience  of  the  uninformed  reader, 
those  basic  definitions  and  properties  of  dyadics  employed  in  the 
body  of  the  report  are  reviewed  here. 

D. 1  Basic  Definitions 

A  dyadic  (more  properly  called  a  dyadic  polynomial)  is  an  oper¬ 
ator  representable  as 

£  ~  f*ibi  +  ®2b2  + . +  a  b  (D-l) 

where  a^a^,...*^  and  b  i ,  b2 , . . .  ^  are  vectors  and  which,  with  an  ar¬ 
bitrary  vector  v,  forms  the  scalar  products 

v*D  *  (v*ai)bi  +  (v*a2)b2  +  .  +  (v«a  )b  (D-2) 

—  —  —  — n  — n 

D*v  *  ai(bi'V)  +  a2(b2»v)  + . +  a  (b  *v)  (D-3) 

—  —  ~  ~  —  — n  — n  — 

The  vectors  {ak>  are  called  antecedents;  the  vectors  {b^}  are  called 
consequents.  Each  term  of  the  dyadic,  a^b^  is  called  a  dyad. 

Dyadics  arise  naturally  in  electromagnetics  when  vector  operators 
are  to  be  separated  mathematically  from  the  field  vectors  being  oper¬ 
ated  upon.  For  example,  the  projection  of  the  electric  field  e  along 
a  linear  antenna  having  the  direction  of  a  unit  vector  u  can  be  writ¬ 
ten  as  u(u»e) .  Here,  the  vector  u  performs  the  projection  operation 
upon  the  electric  field  vector  e.  If  these  two  vectors  are  to  be 
separated  mathematically,  it  is  necessary  to  adopt  the  dyadic  repre¬ 
sentation  uu  in  the  expression  (uu) »e. 


D. 2  Basic  Properties 

Because  in  three-dimensional  vector-space  any  vector  v  can  be 
expressed  in  terms  of  any  three  vectors  a,b,c  forming  a  base  (i.e., 
satisfying  the  condition  a  x  b  *  £  ¥  C) ,  any  dyadic  can  be  written 
as  a  polynomial  of  no  more  them  three  dyads.  For  example,  the  dyadic 
polynomial  of  Equation  (D-l)  can  be  recast  in  the  form 

D=ae+bf+cg  (D-4) 

where 


e  = 

bxc 

[a  b  c] 

*  E 

(D-5) 

f  — 

cxa 

•  E  5ibi 

(D— 6) 

x  — 

[a  k  s) 

a  = 

axb 

[a  b  c] 

•  E  ai^i 

(b-7) 

Since  the  consequents  e,f,c[  can  also  be  expressed  in  terms  of 
an  arbitrary  three-vector  base,  it  is  clear  that  the  dyadic  D  is 
representable  in  terms  of  nine  scalar  coefficients. 

For  example,  if  the  same  three-vector  orthonormal  base  is 
employed  for  both  the  antecedents  and  the  consequents,  D  might 
be  written  in  the  form 

fi  -  dll*°*°  +  di2^#Z#  + 

+  d2l£°£#  +  d22X°E°  +  d2  3^°— #  ( D—  8) 

+  d3lE°*0  +  d32*0Z°  +  d33£°2° 

Although  these  nine  scalar  coefficients  can  be  written  as  a  matrix 
[dij],  a  dyadic  is  not  a  matrix.  The  matrix  [dij]  depends  upon  the 
base  employed  and  so  is  coordinate-dependent;  the  dyadic  is  not. 

A  dyadic  serves  as  a  linear  mapping  from  one  vector  to  another. 
Conversely,  any  such  linear  mapping  can  be  expressed  in  terms  of  a 
dyadic.  This  can  be  demonstrated  by  applying  an  orthonormal  base 
{ui)  to  an  arbitrary  linear  function  f(a)  and  noting  that 

f  (a)  •  f  (UjUj  «a)  (D-9) 

=  CZui  *  *  2.  ID-10) 


D-2 


The  dyadic  in  the  square  brackets  corresponds  to  the  linear  function 
f. 

The  identity  mapping  has  the  property  of  mapping  every  vector 
into  itself.  The  corresponding  dyadic,!:,  is  called  the  unit  dyadic 
or  idemf actor.  Thus, 


a  =  a  =  a  •  I 


Expressing  a  in  an  orthonormal  base  {vy. }  shows  that 


X  *  HlEl  +  H2H2  +  H.3H.3 


More  generally. 


(D-ll) 


(D-12) 


X  =  a  a'  +  b  b'+  c  c 1  (D-13) 

where  (a' ,b',c'}  is  a  base  reciprocal  to  the  arbitrary  base 
{a,b,c>.  Thus,  as  for  any  dyadic,  there  exists  an  infinity  of 
representations  for  the  unit  dyadic.  As  a  consequence  of  Equation 
(D-ll)  the  matrix  of  coefficients  Id^j]  associated  with  a  given  base 
iai }  can  be  found  in  the  relation 


£  =  *  EE  <»•  ’£•«-) >*i-j  “££Dij£i£j 


(D-14) 


D-3 
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The  Fast  Fourier  Transform 

The  fast  Fourier  transform  (FFT)  is  a  highly  efficient  algo¬ 
rithm  for  computing  the  discrete  Fourier  transform  (DFT) .  The  FFT 
can  be  used  in  place  of  the  continuous  Fourier  transform  only  to  the 
extent  that  the  DFT  could  before,  but  with  a  substantial  reduction 
in  computer  time.  This  appendix  reviews  the  basic  properties  of  the 
DFT  and  its  relation  to  the  continuous  Fourier  transform,  introduces 
the  FFT  and  identifies  the  major  pitfalls  likely  to  arise  in  its  use, 
and  exemplifies  the  application  of  the  DFT  and  the  FFT  by  considering 
a  band-limited,  high-frequency,  rectangular  pulse. 

C. 1  Discrete  Fourier  Transform 

The  Fourier  transform  pair  for  continuous  signals  can  be  writ¬ 
ten  in  the  form 


X(f)  *  /°°  X(t)exp{-j2trft)dt 

(E-l) 

to 

x(t)  -  /*X(f)exp{+j2Trft}df 

(E-2) 

-40 


The  integrals  in  Equations  (  E-l)  and  (  E-2)  are  actually  finite  in 
any  practical  case  since  only  a  finite  segment  of  an  essentially 
band-limited  signal  will  be  available.  This  situation  is  commonly 
represented  by  setting  x(t)  «  0  for  0  >  t  >  T  when  the  signal  x(t) 
is  available  only  over  the  time  interval  (0,T),  and  X(f)  *  0  for 
-B/2  >  f  >  B/2  when  the  amplitude  spectrum  X(f)  is  essentially  band- 
limited  to  the  frequency  interval  (-B/2,+B/2). 

The  analogous  discrete  Fourier  transform  (DFT)  pair  that  applies 
to  samplod-data  signals  can  be  written  in  the  form 


X(fa)  *  At  2^  x(tk)exp  {-jSflfjntfc} 
k-0 

(E-3) 

*/2-l 

x ( tjc )  -  Af  £  '  ?i(4>«*p{+j2*f*tk} 

(E-4) 

m— U/2 

E-l 


If  the  signal  x(t)  has  been  sampled  uniformly  at  the  Myguist  rate  so 
that 


At  -  T/N  «  1/2B,  Af  -  2B/N  «  1/T 

then  upon  letting 

t*  *  kAt,  fm  =  mAf 


(E-5) 


(E-6) 


Equations  (e-3)  and  (E-4)  become 

N-l 

X(m)  *  At  2  x(k)exp{-j2Trmk/N}  (E-7) 

N-l 

x(k)  *  Af^  x(m)exp{+j2Trmk/N) 

m«0 

The  periodicity  of  X(m)  apparent  in  Equation  (E-7)  has  been 
used  in  re-ordering  (for  ease  of  computation)  the  summation  appearing 
in  Equation  (e-8) •  More  generally,  however,  it  is  apparent  from 
these  equations  that  both  X(m)  and  x(k)  must  be  considered  periodic, 
i .  e . 


x(.mtN) 

*  X(m) 

(E-9) 

xfk+N) 

"  *00 

(E-10) 

This  property  has  been  anticipated  in  writing  Equation  (E-4)  where 
the  missing  end  point  X(fN^2)  is  considered  tc  be  the  first  point  of 
the  next  period  in  a  periodic  extension  of  X(fn).  Mathematically 
speaking,  the  interval  from  which  samples  of  X(f)  [and  x(t)]  are 
taken  must  be  considered  closed  on  the  left  and  open  on  the  right 
(or  vice-versa).  This  property  is  illustrated  in  Figure  E-l.  Fur¬ 
ther  discussion  of  this  topic  can  be  found  in  [157] . 

When  the  sampled-data  signal  x(k)  is  real,  the  real  part  of 
X(m)  [denoted  hereafter  as  X' (m) ]  is  symmetric  about  the  folding 
frequency  B  [m  *  N/2]  and  the  imaginary  part  of  X(n»)  [denoted  here¬ 
after  as  XH (m)  is  antisymmetric.  Since  X(m)  has  been 


E-2 


interpreted  as  being  periodic,  these  symmetries  are  equivalent  to 
saying  that  X* (m)  is  an  even  function  of  m  and  X" (m)  is  an  odd  func¬ 
tion.  This  also  means  that  the  Fourier  coefficients  between  N/2  and 
N-l  can  be  viewed  as  the  "negative  frequency"  harmonics  between  -N/2 
and  -1,  respectively.  Likewise,  the  last  half  of  the  sampled-data 
signal  x(k)  can  be  interpreted  as  negative  time  (that  is,  as  occur¬ 
ring  before  t  *  C) . 

Since  the  FFT  algorithm  provides  an  efficient  transformation  to 
the  frequency  domain,  it  is  interesting  to  consider  the  product 


Y(m)  «  F(m) X(m) 


(E-ll) 


where  Y(m)  is  the  DPT  of  y(k),  F(m)  is  the  DFT  of  f(k),  e.nd  X(m)  is 
the  DFT  of  x(k).  According  to  Equation  (E-7) 


N-l  ,  , 

y  (k)  *  Af£  F(m)X(m)  exp{j2imk/to} 
m-0 1  J 


(E-12) 


According  to  Equations  (E-5)  and  (E-7)  this  can  be  written  as 


A  A 

T(W  -  ^  E 


At  V"  r|jLx(x)e"j2lT“T/N||N  1t(f)e^2inn^ 


m*o 


T-0 


j2mnk/N 


(E-13) 


Since  all  of  the  sums  are  finite  this  can  be  written 
N-l  N-l  r  N-l 


y (k)  -  £ |f  Nf x(t>f(f)  f  fV2^-^ 

w  t*o  t«o  Lm-o 


Ob/n) 


!E-14) 

This  can  be  simplified  through  the  use  of  the  orthogonality  relation 

N,  if  n=m  (modulo  N) 


N-l 

^  ^2-nlm/H  e“j2irJln/N  _ 


A-0 


0 ,  otherwise 


(E-15) 


Thus,  Equation  (e-14)  is  zero  unless  f  =  k-T  and  hence 


N-l 


y (k)  =  At  £x(.xlf(k-T) 


T-0 


(E-1C) 


Equation  (E-16)  may  be  recognized  as  the  discrete  form  of  the  convo¬ 
lution  integral 


y (t)  «  /  x(t) f (t— t) dx 


V 


Equation  (E-17)  can  be  associated  with  the  response  y(t)  a  fixed 
linear  system  to  the  input  x(t) ,  if  that  linear  system  can  be  char¬ 
acterized  by  the  impulse  response  f(t,T)  where 


f(t,T)  -  *<*-*>.«  t»T 

O,  if  t<T 


(E-18) 


C.2  Fast  Fourier  Transform 

The  fast  Fourier  transform  (FFT)  is  a  highly  efficient  algorithm 
for  computing  the  discrete  Fourier  transform  (DFT) .  Taking  advantage 
of  the  fact  that  the  complex  exponential  weights  appearing  in  the  DFT 
can  be  calculated  iteratively,  the  algorithm  employs  a  clever  compu¬ 
tational  technique  for  sequentially  combining  progressively  larger 
partial  sums  of  the  exponentially  weighted  data  samples  to  realize 
the  DFT.  lvo  mathematically  equivalent  versions  of  the  FFT  algorithm 
usually  provide  the  basis  for  implementation:  the  Cooley-Tukey 
(decimation-in- time)  algorithm  and  the  Sande-Tukey  (decimation-in- 
frequency)  algorithm.  The  algorithm  chosen  for  implementation  is 
usually  selected  to  exploit  the  characteristics  of  the  sampled  data 
and/or  the  hardware-software  properties  of  the  computer  being  used. 


The  value  of  the  fast  Fourier  transformation  lies  in  the  reduc¬ 
tion  of  computer  execution  time  in  evaluating  the  discrete  Fourier 
transform.  Whereas  a  direct  N-point  evaluation  of  the  DFT  typically 
requires  a  computer  execution  time  proportional  to  N2 ,  the  FFT  re¬ 
quires  a  computer  execution  time  proportional  to  only  N'Aog^N.  The 
approximate  ratio  of  FFT  to  direct  execution  time  is  given  by 


N-logaN  =  I022N  (E-19) 

N2  N 

For  example,  if  N  =  256  =  28,  the  FFT  requires  only  about  3%  of  the 
time  required  by  direct  computation. 

Four  problems  often  encountered  in  using  the  FFT  are:  aliasing, 
leakage,  picket-fence  effect,  and  round-off  error.  The  term  "ali¬ 
asing"  refers  to  the  fact  that  high-frequency  spectral  components  of 
a  band-limited  signal  can  "impersonate"  low-frequency  spectral  com¬ 
ponents  if  the  sampling  rate  is  too  low  [153] .  This  circumstance  is 
reflected  in  its  amplitude  spectrum  where  the  sidebands  of  frequency 


E-5 
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translated  baseband  replicas  are  observed  to  overlap.  This  problem 
can  be  eliminated  by  sampling  the  signal  at  a  rate  at  least  twice  as 
high  as  the  highest  frequency  present  in  its  amplitude  spectrum. 
"Leakage"  refers  to  distortion  in  the  amplitude  spectrum  of  a  band- 
limited  signal  arising  as  a  result  of  utilizing  only  a  finite  number 
of  data  samples  from  the  signal's  data  record.  This  problem  can  be 
mitigated  (but  not  eliminated)  either  by  increasing  the  record  length 
of  the  signal  or  by  applying  a  data  window  to  the  data  record  [157] . 
The  "picket-fence  effect"  refers  to  enhanced  spectral  responses  at 
the  discrete  frequencies  of  the  DPT  due  to  interstitial ,  unresolved, 
spectral  components  of  the  signal.  This  problem  can  be  mitigated 
(but  not  eliminated)  by  using  an  interpolation  function  between  the 
discrete  frequencies  of  the  DFT  or  by  extending  the  data  record  arti¬ 
ficially  with  zero  samples  in  order  to  increase  the  spectral  resolu¬ 
tion  of  the  DFT.  "Round-off  error"  arises  in  numerical  computation 
as  a  result  of  the  finite  word  length  employed  in  digital  computers 
for  the  representation  of  numbers.  According  to  Kaneko  and  Lin 
[159,160]  the  total  relative  mean-square  error  (MSE)  due  to  round-off 
error  is  bounded  by 

,-2b  „  2'2b 

n  ^ —  <MSE<3n  —  (E-20) 


where. 


MSE 


< 


N-l 

£  |  e(m)  |  2> 


(E-21) 


e(m)  =  ft(m)-X(m)  (E-22) 

n  *  log2N  (E-23) 


3 


and  &(m)  is  the  computer  value  of  X(m)  based  upon  a  b-bit  representa¬ 
tion  for  the  mantissa  (not  including  sign)  of  a  real  variable  using 
binary  floating-point  arithmetic  with  double  precision  accumulation. 
For  example,  using  double-precision  real  variables  of  4  bytes  each 
(3-byte  mantissa,  1-byte  exponent)  so  that  b  «  23  and  a  512-point 
transform  so  that  n  *  8,  the  total  relative  MSE  is  bounded  by 


W 


■  V 
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(E-24) 


3. 79x10“ 14<MSE<1. 14x10' 13 

The  total  relative  root -me an- square  error  (RMSE)  in  the  amplitude 
spectrum  is  thus  bounded  by 

1.95X10'7  <RMSE<3. 37x10“ 7  (E-25) 

The  relative  accuracy  in  the  calculation  of  the  amplitude  spectrum 
is  thus,  roughly  speaking,  about  1  part  in  107. 

In  passing  it  may  be  noted  that  the  FFT  not  only  reduces  the 

computation  time  of  the  DFT,  it  also  substantially  reduces  the 

round-off  errors  associated  with  direct  computation.  In  fact,  both 

computation  time  and  round-off  error  are  reduced  by  the  same  factor 

of  (£og  N)/N  [156].  Further,  although  the  above  discussion  of  the 
2 

problems  encountered  in  using  the  FFT  was  couched  in  terms  relating 
to  the  forward  trans formation  of  variables  from  time  to  frequency 
[Equation  (e-7)  ] »  the  discussion  remains  valid  in  the  inverse  trans¬ 
formation  of  variables  from  frequency  to  time  [Equation  (E-8)  ] . 

C. 3  Example:  Bandpass  Rectangular  Pulse 

The  basic  properties  of  the  discrete  Fourier  transform  (DFT)  and 
its  computation  by  the  fast  Fourier  transform  (FFT)  may  be  exempli¬ 
fied  by  considering  the  bandpass  rectangular  pulse 

x (t )  -  p(t)  cos( 2irf •  t )  (E-26) 

where  f  is  the  r-f  carrier  frequency  and  p(f)  is  the  pulse  envelope 
o 

defined  by 


p(t) 


jl,  I t | <  L/2 
io,  |t|>  L/2 


(E-27) 


where  L  is  the  pulse  length  (duration) .  In  the  numerical  computa¬ 
tions  which  follow,  f  will  be  taken  to  be  600  MHz  and  L  to  be 

O 

5  nanoseconds  (3  r-f  cycles). 

E-7 


The  amplitude  spectrum  of  the  pul be  envelope  can  be  determined 
from  Equation  (E-l)  to  be 


P(f)  ■  L  •inc{tffL} 


(E-28) 


It  should  be  noted  that  although  the  pulse  p(t)  is  undefined  for 
1 1 1  *  L/2 ,  the  inverse  transformation  of  P(f)  using  Equation  (E-2  ) 
converges  to 


P ( | L/2 | )  -  0,5 


(E-29) 


because  of  Gibb 1 s  phenomenon  [  ] .  As  a  consequence ,  the  numerical 

representation  of  p(t)  at  the  point  of  discontinuity  (|t|  *  L/2) 
must  be  that  afforded  by  Equation  (E-29) .  Failure  to  accomodate 
ribb's  phenomenon  at  points  of  discontinuity  in  the  signal  (or  even 
the  spectrum)  will  lead  to  spurious  oscillations  in  the  transformed 
variable . 

The  r-f  carrier  cos(2Trfot)  and  its  amplitude  spectrum  [two  6 
functions  at  ±fQ]  and  the  pulse  envelope  p(t)  fluid  its  amplitude 
spectrum  P{f)  are  shown,  respectively,  in  Figures  E-2  (a)  fluid  (b). 
The  bandpass  rectangular  pulse  x(t)  [refer  to  Equation  (E-26)  ]  and 
its  flunplitude  spectrum. 


X(f)  »  bine  {til (f+f«)  >  -  sine  (irL(f-f*)}] 

are  shown,  respectively,  in  Figure  E-2  (c) .  In  consonance  with 
Equations  (E-ll)  and  (E-16)  and  the  time- frequency  duality  of  the 
Fourier  trcuisform  [refer  to  Equations  (E-l)  and  (E-2)  j.  Equation 
(E-30)  and  Figure  E-2  (c)  reflect  the  fact  that  the  amplitude  spec¬ 
trum  of  the  product  of  two  time  functions  is  equal  to  the  convolution 
of  their  respective  amplitude  spectra. 

The  sampled-data  signal  x(k)  is  obtained  by  multiplying  the 
bandpass  recteuigular  pulse  x(t)  by  an  infinite  train  of  unit  impulses 
(or  Dirac  comb)  defined  mathematically  by 


c(t)  «  5^6(t-kAt) 

k--*» 


(E-31) 
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The  Dirac  comb  and  its  amplitude  spectrum 

CD 

c<0  ■  sr  2 6 <*-  h->  <*-«>  * 

are  shown  in  Figure  e-2  (d) .  The  sampled-data  signal  x(t)  and  its 
amplitude  spectrum  C(f)*X(f)  are  shown,  respectively r  in  Figure  e-2 
(e> .  It  is  apparent  from  this  figure  that,  because  the  amplitude 
spectrum  of  the  pulse  |P ( f ) ]  is  not  strictly  bandlimited,  aliasing 
cannot  be  entirely  avoided  and  some  spectral  overlapping  is  ines¬ 
capable.  However,  Cyber Com  has  found  that  the  aliasing  associated  ^ 

with  non-bandlimited  signals  (which,  strictly  speaking,  arise  only 
in  mathematical  models  and  not  in  real  systems)  can  be  held  within 
acceptable  bounds,  in  most  cases,  by  letting 


At  -  [2Cf.4B99))“X 


(E-33) 


where  is  the  low-pass  equivalent  bandwidth  containing  99%  of  the 
signal  power.  Note  that  this  rule  is  consistent  with  the  low-pass 
Nyquist  sampling  rate  when  the  signal  is  strictly  bandlimited  to  B 
Hertz.  For  P(f)  [refer  to  Equation  (e-28)  ],  the  99%  containment 
bandwidth  is  [153] 


(E-34) 


In  the  numerical  computation  which  follow 


so  that 


At  « 


(E-35)  * 


(E-36) 


which  corresponds  to  a  sampling  rate  of  8  samples  per  r-f  cycle.  It 
should  be  noted  that  the  aliasing  appearing  in  Figure  E-2  (e)  will 
not  lead  to  distortion  in  recovering  the  sampled-data  signal  x(k) 
from  C(f ) *X( f)  via  the  inverse  Fourier  tranformation.  This  is  as¬ 
sured  by  the  isomorphic  (one-to-one)  relation  between  a  Fourier 
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transform  and  its  inverse.  However,  subsequent  signal  processing 
operations  which  modify  the  amplitude  spectrum  C(f)*X(f)  [e.g., 
convolution]  will  be  affected. 

The  continuous  frequency-domain  function  C(f)*X(f)  shown  in 
Figure  e-2  (e)  can  also  be  made  discrete  (sampled)  by  treating  the 
sample -data  signal  x(k)  as  one  period  of  a  periodic  function  with  the 
data  sample  x(N/2)  taken  as  the  first  sample  of  successive  periods. 
This  forces  both  the  time-domain  and  frequency  domain  functions 
[x(k)  and  X(m) ,  pectively]  to  be  infinite  in  extent,  periodic, 
and  discrete  as  shown  in  Figure  E-2  (f).  The  resolution  in  the  fre¬ 
quency-domain  is  determined  by  the  relation 


Af  - 


1 

NA  t 


(E-37) 


where  T  is  the  period  in  the  time  domain.  At  is  the  sampling  rate, 
and  N  is  the  number  of  data  samples.  It  is  apparent  from  Equation 
(E-37)  that  the  resolution  can  be  increased  by  increasing  the  sam¬ 
pling  rate  and/or  the  number  oi  data  points.  The  number  of  data 
points  can  always  be  increased  without  changing  the  sampling  rate 
by  filling  out  the  (truncated)  data  record  x(k)  with  zeroes.  In 
the  numerical  calculations  which  follow  N  was  taken  to  be  a  power  of 
2,  spe.  fically 


N  =  2®  =  256 


(E-38) 


As  a  consequence  of  Equations  (E-36)  and  (E-37) 


(E-39) 
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APPENDIX  F 


Optical  Theorem  for  Scalar  Fields  in  Two-Dimensional  Media 

The  purpose  of  this  appendix  is  to  apply  Poynting's  vector 
theorem  to  two-dimensional  scatterers  that  are  normally  illuminated. 
This  will  lead  naturally  to  definitions  of  the  total,  scattering 
and  absorption  cross-sections  for  the  scalar  case. 

Consider  a  two-dimensional  scatterer  of  cross-section  S  as 

P 

shown  in  Figure  F-l.  Assume  that  the  scatterer  is  oriented  in  such 
a  way  that  its  generating  elements  are  parallel  to  the  z  axis. 

The  scatterer  is  illuminated  by  a  unit  amplitude  plane  wave. 

Ei<Xt)  =  *  e”^*0— '-t  (F-l) 

where  i*z°  -  0. 

Apply  Poynting's  vector  theorem  to  a  circular  cylindrical 
volume  of  unit  height  and  surface  SQ.  The  cylinder  has  radius  r 
which  is  assumed  to  be  large.  If  E  and  H  are  the  total  electric 
and  magnetic  fields  then 

R  f  S  *da  +  i  f  aE*E*dV  =  0  (F-2) 

6  '  — C  —  Z  Vp  —  — 

where  Sc  is  the  complex  Poynting  vector, 

Sc  =  SExH*  (F-3) 

a  is  the  conductivity  of  the  scatterer 's  material  and  VQ  is  the 
volume  of  the  unit  height  cylinder.  Because  the  scatterer  and 
the  incident  field  are  both  independent  of  z,  it  follows  that 

E  =  E(xt)  =  iMxt)z°  (F-4) 

The  scalar  field  <P(xt)  is  composed  of  an  incident  and  scattered 
portion,  thus 

iMxt)  =  4»i(xt)  +  ^s(xt)  (F-5) 

The  scattering  amplitude,  f(o,i),  is  introduced  when  the  scattered 
wave  is  represented  in  the  far  field  as  follows: 


*S  Ut>  ~  f  (Q/i)e'“jkt-,-t  #  r  «  lxj.1  »  l 

Jr 

where  o  =  xt/r. 

Finding  the  incident  and  scattered  magnetic  fields  from 
Maxwell's  equations  for  large  r  and  then  using  them  in  equation 
(F-2)  results  in  the  following  equation: 


/a'|f<o,i)  I'de+tf-*.  /2"[f(o,i)ejlt°r<co‘e-l) 

O  o  - 

-  j Jc0r  (cos0-l) 

+f*(o,i)e  cos0]  d6  «  oa 

where  the  angle  6  is  defined  in  Figure  F-l  and  the  absorption 
cross-section  oa  is  given  by  [55] 

V  if  /V  °®*£dV  (F- 

The  derivation  is  completed  by  asymptotically  evaluating  the 
the  integral  in  equation  (F-7)  for  large  k0r.  The  integral  can  be 
divided  into  two  integrals  each  of  which  can  be  treated  by  the 
method  of  stationary  phase.  A  calculation  shows  that  there  are 
stationary  points  of  9*0,r  and  2tr.  Now  using  this  asymptotic 
result  in  equation  (F-7) ,  one  obtains 

J2  1 2(16  +  Ref  (i,i)e“jir/4  +  -  0  (F- 

°  ^ 

Because  the  scattering  cross-section  a8  is  defined  as 


°s  =  /o  jf  (o,i)  1 2de 


and  the  total  cross-section  is  the  sum  of  the  scattering  and 
absorption  cross-section  i.e.. 


one  obtains 


o.  -  —  R  [f(i,i)e~jlr/4J  I 

*  *  “  “ 

The  expressions  for  a  and  a  in  terms  of  the  scattering 

8  t 

amplitude  f  can  be  put  in  terms  of  the  Fourier  transform  of  the 
transition  operator  t.  By  using  equation  (2-2-24)  in  equations 
(F-10)  and  (F-12) ,  the  following  expressions  are  obtained: 


APPENDIX  G 


Alternative  Representation  for  the  Intensity  of  Scalar  Fields 

in  Two-Dimensional  Media 


In  Chapter  7  an  integral  representation  for  the  normalized 
intensity  [equation  (7-28)3  has  been  obtained.  Unfortunately,  the 
integral  is  very  slowly  convergent  for  large  values  of  tq.  To 
rectify  the  difficulty,  the  integral  will  be  converted  into  an 
alternate  form  which  is  more  rapidly  convergent. 

For  convenience,  the  original  expression  for  T(t)  is  repeated 
here.  It  is  given  by 


I ( t ) =ie 


-Wc 


(G-l) 


The  alternative  representation  for  T(t)  can  be  obtained  by  deforming 
the  integration  path  in  the  q  plane.  To  do  this  an  integral  from 
-®  to  +«*>  is  needed.  This  is  obtained  by  representing  J0(xq)  in 
terms  of  the  Hankel  function. 


J  (Tq)  =  [H^Mtq)  +H*2\iq)]/2 


(G-2) 


and  then  substituting  this  expression  into  equation  (G-l) 
This  gives 


I  (T) 


/ 


qH^>1)  (Tq) 


-  W 


(G-3) 


♦  if /■  aqqH°,21(Tq) 

2  ^!+T  -  W0 


The  last  integral  in  equation  (G-3)  is  now  transformed  by  letting 
q'  *  qe***-6^,  e  >  0.  The  result  when  combined  with  the  first 
integral  of  equation  (G-3)  is  given  by 


I(t) 


<2Holl>(Tq) 

'  -  v* 


(G-4) 


G-l 


where  P  is  the  resulting  integration  path  shown  in  Figure  G-l. 

The  small  parameter  e  has  been  introduced  in  order  to  avoid  the 
branch-cut  of  the  Hankel  function  lying  along  the  negative  real 
q  axis. 

An  examination  of  the  singularities  of  equation  (G-4)  shows 
that  in  addition  to  the  Hankel  function  branch-cut,  there  are 
simple  poles  at  q  *  q±  *  +  i  vl-W|  and  branch  point  singularities 
at  q  «  ±i.  The  branch  cuts  associated  with  these  branch  point 
singularities  are  chosen  alo;  cj  the  positive  and  negative  imaginary 
axes.  The  alternative  form  of  the  integral  is  now  obtained  by 
deforming  the  integration  path  P  into  the  upper  half  of  the  q  plane. 
Because  the  integrand  tends  to  zero  as  |q|  -*■  <*>  in  the  upper  half 
q  plane,  the  inte  ral  is  represented  wholly  by  the  contribution  due 
to  the  pole  q  *  q±  and  the  upper  branch-cut.  The  result  is 


I (t)  =  2tt1tW  exH<15  (it vd-W7) /2 
o  o  o 

(1) 


reT  r  dq  q  H° 

Jc1+c2 


(xq) 


(G-5) 


✓q*+T  -  W, 


where  the  paths  C1  and  C2  are  defined  in  Figure  G-l. 

The  integrals  over  and  C2  can  each  be  transformed  to  the 
real  positive  axis  and  combined.  The  result  is 

I(t)  *  2W  teTK  (t  ^L-w“ ) 


«  K  (x  »l+pl) 
+  —  eV  dp  p2— - — 

71  °  p2+  w2 


(G-6) 


Here  K0(z)  is  the  modified  Hankel  function.  Because  K0(z)  is 
decaying  exponentially  for  large  z,  the  integral  converges  rapidly 
for  large  values  of  p. 
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